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APRESENTAÇÃO 


O ramo das engenharias veem ganhando cada vez mais espaço no decorrer dos 
anos, sendo hoje um dos principais pilares para o setor empresarial. Analisar os campos 
de atuação, bem como pontos de inserção e melhoria dessa desta área é de grande 
importância, buscando desenvolver novos métodos e ferramentas para melhoria continua 
de processos. 

Desta forma estudar temas relacionados a engenharia é de grande importância, 
pois desta maneira pode-se aprimorar os conceitos e aplicar os mesmos de maneira mais 
eficaz. O aumento no interesse se dá principalmente pela escassez de matérias primas, a 
necessidade de novos materiais que possuam melhores características físicas e químicas 
e a necessidade de reaproveitamento dos resíduos em geral. Além disso a busca pela 
otimização no desenvolvimento de projetos, leva cada vez mais a simulação de processos, 
buscando uma redução de custos e de tempo. 

Neste livro são apresentados trabalho teóricos e práticos, relacionados a área de 
engenharia, dando um panorama dos assuntos em pesquisa atualmente. De abordagem 
objetiva, a obra se mostra de grande relevância para graduandos, alunos de pós-graduação, 
docentes e profissionais, apresentando temáticas e metodologias diversificadas, em 
situações reais. Sendo hoje que utilizar dos conhecimentos científicos de uma maneira 
eficaz e eficiente é um dos desafios dos novos engenheiros. 
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ANÁLISE MODAL E REPRESENTAÇÃO MUSICAL 
APLICADAS AO DESIGN DE ESTÚDIOS DE 
PRODUÇÃO FONOGRÁFICA 


Data de aceite: 01/11/2021 
Data de submissão: 20/09/2021 


José Augusto Mannis 

Universidade Estadual de Campinas 
(Unicamp), Instituto de Artes, Depto. de Música, 
Laboratório de Acústica e Artes Sonoras 
(LASom) 

Campinas — SP 
http:/lattes.cnpg.br/6940068973325101 
htitps://orcid.org/0000-0002-9484-5674 
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RESUMO: Este trabalho relata e desenvolve 
uma análise acústica realizada em 2015, 
apresentando a fundamentação teórica e a 
metodologia adotada; discutindo resultados 
obtidos e apontando para novas perspectivas de 
pesquisa. Como métodos inovadores destacam- 
se a análise do comportamento modal em acústica 
de salas pelo Critério de Bonello a partir de 
resultados obtidos por simulação por elementos 
finitos; a inclusão de representação musical para 
dar mais visibilidade a imagens mentais próprias 
a escuta humana. Medições acústicas realizadas 
nas duas salas simuladas, depois de construídas 
e devidamente instaladas de mobiliário e 
equipamentos, permitiram estender a análise a 


outros aspectos de interesse como a disposição 
dos monitores de áudio. Este estudo é uma 
segunda versão revisada de palestra ministrada 
14º Congresso de Engenharia de Áudio realizado 
pela Sociedade Brasileira de Engenharia de 
Audio (AES Brasil) em maio de 2016, em São 
Paulo sendo originalmente publicada nos Anais 
deste evento. 

PALAVRAS-CHAVE: Acústica de salas; Conforto 
ambiental; Engenharia de áudio; Estúdio de 
produção fonográfica; Sala de audição crítica. 


MODAL ANALYSIS AND MUSICAL 
REPRESENTATION APPLIED TO AUDIO 
RECORDING STUDIO DESIGN 


ABSTRACT: This work describes and develops an 
acoustic analysis carried out in 2015, presenting 
a theoretical foundation and the methodology 
adopted; discussing the results obtained and 
pointing to new research perspectives. As 
innovative methods, the analysis of modal 
behaviour in room acoustics using the Bonello 
Criteria ffom results obtained by finite element 
simulation results are highlighted; the inclusion of 
musical representation give us more visibility to 
mental images proper to human listening. Acoustic 
measurements carried out in the two simulated 
rooms, after building and properly installed 
furniture and equipment, allowed the analysis to 
be extended to other aspects of interest such as 
the arrangement of audio monitors. This study is 
a second revised version of a lecture given 14h 
Congress of Audio Engineering conducted by 
the Brazilian Society of Audio Engineering (AES 
Brazil) in May 2016 in Sao Paulo and originally 


1 http://aesbrasil.org.brAwp-content/uploads/2018/anais/Anais AESBR2016.pdf 
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published in the Proceedings of the event. 
KEYWORDS: Room acoustics; Environmental comfort; Audio engineering; Phonographic 
production studio; Critical listening room. 


11 INTRODUÇÃO 


Os comprimentos de onda das baixas frequências que caracterizam as salas são 
comparáveis às suas próprias dimensões. Isso significa que temos muito mais frentes de 
onda que se deslocam, do que raios acústicos. Segundo Gerges, a teoria de geometria ou 
raio acústico é inadequada para estudar o comportamento acústico de salas para baixas 
frequências. A mais adequada é a teoria da análise modal, na qual a sala pode ser tratada 
como um ressonador complexo comportando vários modos acústicos, cada um com uma 
frequência característica de ressonância e amortecimento (GERGES, 2000). Os modos de 
vibração em acústica de salas estão especificamente vinculados às ondas estacionárias 
que se estabelecem em seu volume interior única e exclusivamente por propagação aérea. 
Existem três tipos de modos característicos de uma sala retangular: (1) Modos Axiais, 
referentes a ondas estacionarias resultantes da propagação na direção de uma única linha, 
ou melhor, reta, portanto, entre superfícies paralelas. Se a sala for perfeitamente retangular 
e a origem estiver num de seus vértices, os modos axiais se alinharão aos eixos (x, y 
e z); (2) Modos Tangenciais, referentes a ondas estacionárias que se estabelecem pela 
propagação sonora em duas direções, seguindo não mais uma linha, mas agora um plano 
definido por duas retas, como por exemplo os planos xy, xz e yz; (3) Modos Oblíquos nos 
quais as ondas estacionárias resultam de propagações nas três direções ao mesmo tempo. 


Axial Tangencial Oblíquo 


Figura 1 - Modos de vibração de uma sala produtos do percurso de propagação de ondas estacionárias 
em relação à geometria da sala: (1) Axial — em uma única direção, portanto em uma reta (eixo da 
sala); (2) Tangencial — em um plano: propagação paralela à superfícies de seu contorno; (3) Oblíquo — 
propagação em todas as direções (EVEREST, 1988). 
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Figura 2 - Alturas musicais (C, D, E, F. GA, B— do, ré, mi, fá, sol, lá, si) com as respectivas alterações: 
sustenido (  ) um semitom acima; bemol (>) um semitom abaixo; bequadro (:) cancela cada um dos 
anteriores. 


As alturas musicais inscritas em pentagrama nas representações empregadas neste 
trabalho, não se referem a nenhum sinal gerado por fonte sonora, mas apenas a potências 
de ressonâncias que somente ocorrerão mediante a incidência de ondas sonoras com 
frequências correspondentes. 

Cabe aqui ainda explicitar a distinção a ser feita nas áreas de conhecimento de 
acústica, psicologia e afins entre os conceitos de frequência, grandeza atribuída a uma 
característica da onda sonora enquanto fenômeno físico, e altura, sensação produto de 
um fenômeno da percepção. Cada um desses conceitos remete, portanto, a uma instância 
distinta. O primeiro refere-se a um fenômeno físico externo que se pretende observar, 
enquanto o segundo, a um fenômeno subjetivo, na instância do observador, como uma 
sensação que se apresenta à sua consciência. Portanto atenção para que a associação 
de alturas aos modos de vibração de uma sala não induza ninguém a ter a equivocada 
ideia de que uma sala possa gerar sons de alguma maneira. Salas por si não geram sons, 
elas simplesmente soam e ressoam. Os modos de vibração de uma sala farão com que 
determinadas ondas sonoras permaneçam ressoando por mais tempo do que outras, 
simplesmente porque sua conformação geométrica favorece em alguns casos encontros de 
máximos e mínimos de pressão que se somam, e em outros cancelamentos entre máximos 
e mínimos que tendem a ocorrerem no mesmo lugar. Por essa simples razão é que o 
principal elemento de design acústico de uma sala é antes de tudo a cuidadosa definição 
de suas proporções, ou seja, seu volume em um contorno geral delineado, harmonizando 
altura, largura e profundidade. 

Nomenclatura dos modos: Os modos obtidos na simulação de cada sala estão 
ordenados de M001 a M200. Os modos axiais entre teto e piso levam a letra Hacompanhada 
de um número n referente à sua ordem: H1 (primeiro Modo Axial entre teto e piso). 


21 METODOLOGIA 


A análise modal deste estudo não foi realizada a partir de predição dos modos 
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inferida sobre uma superfície regular com faces paralelas (SRP), um paralelepípedo, com 
dimensões totais médias de largura (L ), profundidade (L) e altura (L ), mas através de um 
aplicativo de elementos finitos simulando o comportamento dos modos de vibração de um 
volume de ar contido no contorno de uma superfície irregular (SIR) com admitância nula, 
ou seja, considerando todas as superfícies como totalmente reflexivas. O método a partir 
das dimensões medidas consiste por exemplo em obter um contorno equivalente como na 
figura a seguir, compensando as irregularidades da superfície de contorno equilibrando os 


pontos além e aquém da linha média, mantendo a área e o volume constantes. 





Figura 3 — Exemplo de duas tentativas de soluções de contornos médios, cada um reduzido a dois 
pares de linhas paralelas, buscando equilibrar as áreas que excedem e as que faltam em relação à 
figura inicial, de maneira a que se chegue a uma área equivalente. 





Figura 4 — Modelo de sala retangular (paralelepípedo) de referência para compreensão dos Modos 
Normais de Vibração de uma sala. (DAVIS, 2006, p. 180) 


Uma vez obtidas as três dimensões médias da sala L,, L e L, os modos de vibração 


no volume de ar no interior do contorno médio é dado, como proposto por Davis (2006), por 





Equação 1 


Sendo, 
c, a velocidade de propagação do som 
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N. N,e N, são números inteiros a partir de 0, 1, 2, 3... 

O diferencial significativo do presente método, baseado em simulação em aplicativo 
de elementos finitos, está em lidar diretamente com modos efetivos que se verificam na 
sala, e não com modos de um modelo aproximado como o que acabou de ser exposto logo 


acima. 


31 PREPARAÇÃO DOS DADOS PARA SIMULAÇÃO 


Realizamos a análise modal de dois projetos de salas de controle de áudio multicanal 
denominados Sala de controle de áudio A e Sala de controle de áudio B, com as seguintes 
características: 





Área: 34,6 m? 
Perímetro: 22,3 m 
Área total: 148,6 m? 
Volume: 123,2 mº 


Tabela 1 — Coordenadas (em metros) dos vértices localizados na estrutura de alvenaria de cada sala 
(piso, paredes e laje) da Sala de controle de áudio A 
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Área: 28,7 m? 
Perímetro: 21,4 m 
Área total: 120,5 m? 
Volume: 84,7 mé 


Tabela 2 - Coordenadas (em metros) dos vértices localizados na estrutura de alvenaria de cada sala 
(piso, paredes e laje) da Sala de controle de áudio B 
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Admitância da superfície de contorno 


Elemento da malha de elementos finitos 
(A 16) 


max 


Frequência máxima a ser considerada 
(t=21ºC c=343,78m/s? À «= 1,32M) 


X 





Frequência Final 3000 Hz 


Tabela 3 - Parâmetros aplicados para a análise modal - simulação em aplicativo de elementos finitos 


41 RESULTADOS, ANÁLISE E DISCUSSÃO 


Os dados resultantes foram (1) uma listagem com os Modos ordenados de 1 a 300 
e suas respectivas frequências; (2) plotagem em perspectiva compreendo três dimensões 
da conformação da variação de pressão no volume de cada sala para cada um dos seus 
modos. Foram desconsiderados os modos com frequência acima da frequência máxima, 
sendo esta delimitada pelo elemento empregado para compor a malha estrutural (22 
cm) que não deve exceder 1/6 do comprimento de onda da onda estacionária simulada. 
Finalmente foram considerados os primeiros 270 modos da Sala A e 196 modos da Sala B. 


2 Velocidade de propagação do som definida em função da temperatura conforme ISO 9613-1:1993 Acoustics — At- 
tenuation of sound during propagation outdoors — Part 1: Calculation of the absorption of sound by the atmosphere 
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SET FREQ (Hz) 
27,71 
31,27 
48,35 
49,40 
50,74 
55,77 
57,63 
60,81 
64,34 
69,23 
70,20 
76,12 
77,84 
79,44 
80,62 
81,82 
90,36 
92,86 
92,97 
93,19 
95,22 
97,06 
101,00 
101,04 
102,09 
104,29 
104,91 
105,02 
109,16 
109,81 
110,17 
112,22 
112,35 
114,94 
115,20 
116,85 
119,53 
119,59 
120,58 
122,57 


FREQ (Hz) 
153,72 
153,90 
155,04 
155,07 
155,11 
155,47 
155,52 
156,74 
157,49 
157,94 
159,36 
160,72 
161,47 
163,14 
164,42 
164,76 
165,10 
165,82 
166,04 
166,28 
167,71 
168,81 
169,95 
170,45 
170,97 
171,31 
173,11 
173,25 
174,59 
174,83 
174,96 
175,82 
176,56 


FREQ (Hz) 
199,85 
200,09 
200,30 
200,45 
200,71 
202,54 
203,22 
208,51 
203,63 
203,71 
204,12 
204,16 
204,96 
205,54 
207,33 
207,79 
208,17 
208,18 
208,28 
209,32 
209,49 
209,62 
209,67 
212,27 
212,54 
213,28 
214,00 
214,10 
214,79 
215,12 
215,20 
215,21 
215,52 

176,70 216,47 

177,00 216,52 

178,05 170 217,18 

179,43 171 217,68 

180,25 172 217,75 

181,37 173 218,16 

182,16 174 218,56 
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127,58 
129,26 
129,28 
131,62 


138,56 
139,55 
140,54 
140,90 
142,34 


183,53 
184,71 
185,12 
185,91 
186,40 


191,16 
191,29 
192,95 
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Tabela 5 - Frequências dos 200 primeiros Modos de Vibração da Sala B obtidas na simulação. 


Everest (1988, p. 58-59) e Davis (1988, p. 343-344) citam os critérios desenvolvidos 
por Oscar Bonello (1979; 1979b; 1981), Bonello Criteria, para avaliação da qualidade acústica 
de uma sala no tocante a seus modos de vibração. Bonello prescreve a observação e análise 
da configuração e do comportamento dos modos de vibração de uma sala. O Critério 1 
verifica como os modos estão distribuídos, mais particularmente quanto à densidade modal 
em função da frequência. Para verificar a progressão da acumulação de modos, monta-se 
um diagrama contendo o número de frequências modais por banda de frequência de um 
terço de oitava e observa-se como o gráfico evolui. Se for progressivamente crescente, de 
forma regular, as frequências modais estarão com uma distribuição aceitável. Entretanto, 


se a linha for irregular, poderá haver problemas de coloração na resposta acústica da sala. 
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Figura 7 — Análise pelo Critério de Bonello das frequências dos Modos de Vibração da Sala A obtidas 
na simulação. 


Na análise através do Critério 1 de Bonello constata-se que a Sala A apresenta uma 
lacuna na densidade modal na banda de 40 Hz e uma forte contenção de progressão de 
densidade modal nas bandas de 63 e 80 Hz. 

Na análise da Sala B pelo Critério de Bonello (Figura 8) temos uma progressão com 
uma falsa lacuna [FAUX GAP] em 32 Hz, portanto uma progressão contínua, havendo uma 
leve contenção na progressão da densidade modal na banda 80 Hz após a qual, de 90 a 
180 Hz, ocorre uma retomada na progressão. A aparente ausência de modo na banda de 
32 Hz corresponde a uma falsa lacuna, pois os modos iniciais estão em constante e regular 
progressão, tendo ocorrido de ter dois dos modos posicionados, um muito próximo, mas 
antecedendo o início da banda de 32 Hz e o seguinte, sucedendo imediatamente o final da 
mesma banda, conforme se observa na Figura 10. 

Nesta análise de comportamento da progressão modal, empregaremos, dentre 
outros procedimentos, um método apresentado por Renato Cipriano na AES de 2010, 
infelizmente em palestra sem publicação, mas cujo processo consiste em distribuir os 
modos em um diagrama logarítmico de maneira a ter uma impressão mais próxima da 
sensação auditiva quanto à distribuição dos modos em função das alturas percebidas pelo 
ouvido humano, como apresentado em Figura 9e Figura 10. 
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Figura 8 - Análise pelo Critério de Bonello das frequências dos Modos de Vibração da Sala B obtidas na 
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Figura 9 - Distribuição das frequências dos Modos de Vibração da Sala A (16 — 128 Hz) obtidas na 
simulação, plotadas em eixo logarítmico. 


A Figura 9 explicita o vazio que se forma na região de 40 Hz [GAP 40] da Figura 
7 bem como a progressão fortemente contida da acumulação de modos nas bandas de 
63 e 80 Hz [FLOAT 63-80], aglutinando vários modos muito próximos um do outro em 
pequenos grupos, cercados por segmentos vazios, proporcionalmente iguais ou maiores 
que o próprio agrupamento. 

Observando na Figura 10 como o ritmo dos modos representados pelos pontos que 
se acumulam progride regularmente, pode-se concluir que a lacuna em 32 Hz no gráfico da 
Figura 8 é unicamente produto de um efeito de fase entre o período das bandas e do passo 
de progressão modal, não reportando a uma irregularidade de progressão. Por efeitos 
como esse é sempre recomendável checar as possibilidades de predição de qualidade com 
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diversos métodos. Quando todos apontam para a mesma tendência é uma confirmação. 
Quando há contradições, aprofunda-se a verificação. 

Discrepâncias como essa em relação à análise pelo Critério de Bonello ou ainda 
efeitos contrários ao expresso no gráfico de progressão de densidade modal encontrados 
na posição do sweet spot podem ser evidenciados somente após o exame das plotagens 
de pressão para cada modo resultantes da análise modal com simulação. Evidências sobre 
essas discrepâncias encontram-se no item 4, mais à frente. 
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Figura 10 - Distribuição das frequências dos Modos de Vibração da Sala B (16 — 128 Hz) obtidas na 
simulação. 


A Sala B, portanto, começa a progressão modal de maneira regular e suave, 
contrastando neste aspecto com o GAP 40 Hz, o FLOAT 63-80 e a retomada da progressão 
nas bandas de 90 a 180 Hz da Sala A, plotada na Figura 11, apresentando truncamentos 
irregulares em pequenos segmentos. Os agrupamentos inconstantes tendem a reforçar a 
intensidade nas frequências dos modos reunidos em detrimento de uma menor intensidade 
nas bandas de frequências vizinhas, proporcionalmente descobertas, sem a presença de 
modos. Razão pela qual a isso denominamos flutuação. 
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Figura 11 - Distribuição das frequências dos Modos de Vibração da Sala B (90 — 180 Hz) obtidas na 
simulação. 


Buscando aproximar ainda mais a apresentação dos dados de maneira semelhante 
a sensação auditiva, as frequências dos modos foram convertidas em alturas musicais 
com precisão de 1/8 de tomº, através da comparação da frequência de cada modo com 
uma tabela de frequências com precisão de 1/16 de tom, sendo cada altura musical de 
referência centrada e compreendida num âmbito com intervalo de +1/16 de tom. 


ELUPJPMHJFPÃÕJÊPMÃ ÇMPÃÕ1IHOP HM NM 
e E E A ES IS 





Bequadro: 1/8tom T/4tom 3/8 tom Sustenido: 5/8tom 3/4 tom 7/8 tom 
sem alteração — acima acima acima 1/2 tom acima acima acima acima 


Figura 12 — Legenda da notação de alturas musicais em intervalos de 1/8 de tom. A nota de referência 
no pentagrama é um lá 3 (A4) 440 Hz. 
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Figura 13 — Representação em notação musical dos 50 primeiros modos da Sala A. Abaixo de cada 
nota musical, a indicação da ordem do modo e sua respectiva frequência. 


Contudo, apesar de nesta última representação ter havido um sensível e 
significativo ganho de precisão em alturas, a percepção das irregularidades de distribuição 
das frequências ficou menos clara do que estava na plotagem sobre eixo de frequência 
logarítmico (Figura 9, Figura 10, Figura 11). Como intuito de torná-la novamente perceptível, 
atribuímos à posição das notas musicais sequenciadas no pentagrama não mais a ordem do 
modo (1, 2, 3...) cuja frequência lhe corresponde, mas uma progressão constante de alturas 
em intervalos de 1/8 de tom. Cada posição à direita de cada nota, está reservada para 
uma nota 1/8 de tom acima, de maneira que percorrendo horizontalmente o pentagrama 
temos as posições para todas as alturas ordenadas em sequência de 1/8 de tom. Porém, só 
estarão inscritas as notas cuja altura corresponder à frequência de um determinado modo 
encontrado na sala, caso contrário o espaço permanecerá vazio. Assim, nas posições em 
que houver uma nota inscrita, terá como significado que há um modo com a frequência 
correspondendo à altura anotada. Quando não houver nenhum modo a ser associado, o 
espaço permanece vazio, ou seja, com a representação musical de pausa ( - ). Os espaços 
vazios no pentagrama corresponderão à progressão de frequências cujas ondas sonoras 
não serão influenciadas por ondas estacionárias. Além disso, acrescentamos à cada nota, 
sempre que possível, a imagem da plotagem de pressão correspondente a seu modo de 
vibração, conforme obtido na simulação. 

É importante salientar que todas as zonas de máximos e mínimos de pressão 
mostradas nas plotagens produto da simulação são potências de respostas da sala, 
podendo ser mais ou menos evidenciadas em função do posicionamento dos monitores de 
áudio, como mostrado no item. Ou seja, todos os modos da sala aqui mapeados, podem 
ser excitados colocando fontes sonoras em várias posições, tendo uma resposta sonora 
em frequência diferente para cada posição. Porém é importante salientar que todas essas 


respostas serão derivadas da mesma estrutura modal da sala, representada pelo produto 
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do mapeamento, cujo processo está aqui apresentado e proposto como metodologia de 


análise. 
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Figura 14 — Representação da progressão modal da Sala A em alturas musicais em 1/8 de tom (escala 
de 1 a 35 oitavos de tom) — Modos MO1 e MO2Z. 
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Figura 15 - Representação da progressão modal da Sala A em alturas musicais em 1/8 de tom (escala 
de 36 a 71 oitavos de tom) — Modos MOS a M12. 
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Julgamos relevante priorizar a análise e o tratamento acústico em relação ao sweet 
spot da sala de controle de áudio, pois esta é a posição de escuta mais relevante nesse 
tipo de sala. Após a identificação de todos os modos interferindo com máximos e mínimos 
de pressão na zona do sweet spot, eles foram reunidos e distribuídos em representação de 
alturas musicais em sequência de passos de 1/8 de tom, evidenciando assim as incidências 
e os vazios de modos, numa representação de sensação linear para a percepção humana e 
bastante clara à interpretação e análise. Essa representação implementa um novo recurso 
permitindo antever a variação de resposta acústica da sala para cada nota (altura) musical 
e, com isso, tornar possível um controle estreito e preciso por parte da equipe de produção, 
integrando o conteúdo musical às técnicas de produção fonográfica. Esse esquema 
também representa um mapeamento da resposta acústica da sala de controle em função 
das alturas musicais. Se a tonalidade principal da música gravada coincidir com frequência 
de um modo, alguns efeitos, como por exemplo o de finalização em cadência podem ser 
potencializados. Se uma eventual influência desse tipo proporcionar uma impressão de 
maior plenitude e amplidão, isso estaria falseando a mixagem, pois é um efeito acústico 
ocorrendo na sala e não no processo de mixagem. Da Sala A separamos na Figura 18 
quatro modos que nos pareceram particularmente críticos: MO8, M14, M22 e M60. Todos 
atuam fortemente no sweet spot e na altura da cabeça humana, em posturas em pé ou 


sentado. 
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Figura 16 - Representação da progressão modal da Sala A em alturas musicais em 1/8 de tom (escala 
de 72 a 107 oitavos de tom) — Modos M13 a M39. 
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Figura 17 — Distribuição em passos de 1/8 de tom das alturas dos modos com impacto sobre o sweet 
spot da Sala A. 
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Na concepção de contorno de sala em forma de polígono simétrico, como é o caso 
da Sala A, os modos axiais laterais sempre tenderão a atuar na direção do sweet spot 
como M0O2, MO8 e M37, neste caso reforçando frequências de modos associadas a alturas 
musicais la, 1/8 de tom acima, si, 1/8 de tom acima, alturas entre sib, /4 de tom abaixo e 
do, 1/8 de tom abaixo, e entre la, e si, (Figura 13), além de envolver percursos lineares de 
propagação das estacionárias na composição de modos tangenciais próximos como M12 
(mi, /4 de tom abaixo), M1S (mib,), M18e M19 (fáf,) (Figura 14, Figura 15, Figura 16). Essa 
característica de resposta é inerente à conformação do contorno e só pode ser alterada, 
mudando sua forma ou implementando bass trap com interferência em % de comprimento 
de onda no eixo da dimensão lateral mais ampla. O dimensionamento do bass trap, 
compreendendo sua efetiva eficiência, pode ser verificado em simulação em aplicativo de 
elementos finitos. Um outro ponto que chama a atenção na Sala A é a grande concentração 
de modos acerca de si, e sib, /4 de tom abaixo e si, e sib, /4 de tom abaixo, pelo fato de ter, 
na raiz dos modos, três frequências fundamentais próximas decorrentes das proporções 
entre as dimensões da sala: Profundidade M1 (la, 1/8 de tom acima), Largura M2 (si 1/8 de 
tom acima) e altura MS (sol, /4 de tom abaixo ou fa, % de tom acima). 





MODOS : 60,81 Hz MOLA: 79,44 Hz 





MOZa : 47,06 Hz MD6O : 146,57 Hz 


Figura 18 — Sala A - Impacto no sweet spot: (1) Modo de Vibração Axial L1: primeiro modo lateral. MO08 
60,81 Hz; (2) Modo de Vibração Tangencial: M014 79,44 Hz; (3) Modo de Vibração Axial H2: MO22 
97,06 Hz; (4) Modo de Vibração Axial H3: M060 146,57 Hz. 
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Figura 19 — Representação musical das Alturas correspondendo as frequências dos primeiros modos 
axiais da Sala A. Pentagramas inferior: (1) altura real e (2) razão proporcional entre as frequências 
(intervalos musicais) sendo a mesma para as dimensões dos comprimentos de onda entre x, y e z. 

Pentagrama superior: (3) representação das alturas na configuração mais compacta possível. (4) 
intervalo entre a nota mais aguda e a mais grave desse agrupamento, em unidades de 1/8 de tom. 


Observa-se na parte superior da Figura 19 três alturas musicais concentradas em 
apenas um intervalo de terça maior mais % de tom, equivalente a 18 oitavos de tom, ou 
seja, 1/3 da oitava. Já a Sala B apresenta uma relação entre dimensões mais espaçada, 
tendo na representação compacta um intervalo de 21 oitavos de tom entre a nota mais 


aguda e a mais grave. 


Sala B 





Figura 20 — Representação musical das Alturas correspondendo às frequências dos primeiros modos 
axiais da Sala B. 


Se houver correlação para a qualidade da resposta modal entre a configuração 
de alturas sonoras e as proporções das dimensões dos comprimentos de onda dos três 
primeiros modos axiais poderíamos investigar, mas em outro trabalho, a hipótese de que 
a constituição desse acorde compacto, quanto mais ela for ampla, mais equilibradamente 
estariam distribuídos os modos, com menos coincidências de alturas sobrepostas. 
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Comparando a progressão modal de Salas Retangulares (SRT), dada pela Equação 1, 
com outras salas com superfícies irregulares (SIR) como as Salas A e B aqui estudadas, 
com pé direito H constante, observa-se que os modos axiais apresentam frequências mais 
elevadas do que as esperadas para as frequências múltiplas inteiras da frequência de seu 
modo axial H7, como se a altura da sala diminuísse virtualmente em relação à altura (pé 
direito) real. No caso da Sala Ao primeiro Modo Axial H1 corresponde ao modo M003 48,35 
Hz cujo comprimento de onda a 21º C é 7,11 m praticamente duas vezes o pé direito da 
sala, 3,56 m, com um desvio de 0,01 m. 
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Alturas dos Modos potencialmente interferindo 
em região acerda do SWEETSPOT - sala B 
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Figura 21 - Distribuição em passos de 1/8 de tom das alturas dos modos com impacto sobre o sweet 
spot da Sala B. 


ANSTS 





MOOS: 48,35Hz 
Figura 22 - Sala A: Modo de Vibração Axial H1: MOOS 48,35 Hz. 


O segundo modo axial vertical H2 é o M022 97,06 Hz (Figura 18 - 3) cujo 
comprimento de onda a 21º C é 3,54 m. Seu H1 virtual seria de 48,53 Hz correspondendo a 
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um comprimento de onda de 7,08 m. 97,06 Hz para o modo axial H2 é ligeiramente superior 
a frequência esperada para o segundo modo axial H2 96,70 Hz (2 x H1). O terceiro modo 
axial vertical H3 146,57 Hz (Figura 18 - 4), correspondendo ao M060, está também acima 
do H3 esperado 3 x H7 = 145,05 Hz, assim como H4 197,18 Hz (M128) (Figura 24) está 
acima do esperado 4 x H1 = 193,40 Hz. Esses desvios encontram-se em progressão linear, 
conforme apresentado a seguir (Figura 23). Observa-se que o modo axial H4 (Figura 24), 
correspondendo ao modo M128 197,18 Hz apresenta na conformação de seus máximos 
e mínimos de pressão a influência de modos tangenciais, pois as camadas em si estão 
delineadas em alturas estáveis e as concentrações de máximos (vermelhos) se dão 
em áreas definidas em um único plano horizontal, intercaladas de camadas estáveis de 
mínimos (azuis). Essa interferência de modos tangenciais pode estar vinculada à resposta 
modal causada pela geometria irregular de contorno das paredes laterais, neste caso mais 
complexa do que numa sala retangular, desestabilizando progressivamente os modos 
axiais. Isso poderia ser verificado em posteriores trabalhos e seria uma influência bastante 
positiva, pois os modos axiais são os mais intensos e potencialmente problemáticos 


podendo induzir ao desequilíbrio na resposta em frequência. 


Intervalo de Aumento da Frequencia de cada 
Modo em relacao a F, 48,35 Hz 


1,02 y = 0,0079x + 0,9954 





URSS) 
MODO AXIALh 2 MODO AXIALh 3 MODO AXIAL h 4 


Figura 23 — Intervalo crescente de aumento da frequência dos Modos Axiais Verticais (H2, H3, H4) em 
relação à frequência do Modo Axial H1 48,35 Hz, modo MOOS da Sala A. 
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NODEL SOLUTION 





M128:197,18Hz 
Figura 24 - Sala A: Modo de Vibração Axial H4: com influência do modo tangencial M128 197,18 Hz. 


Medições realizadas ulteriormente, depois dos estúdios concluídos, apontam para 
discrepâncias em relação à análise do Critério de Bonello, no tocante à progressão da 
densidade modal, mesmo para as frequências produto da simulação. Alguns máximos de 
pressão localizados próximo ao sweet spot podem sugerir tendências contrárias as do 
Critério de Bonello, como ocorre nas medições que resultaram nas Curvas de Respostas 
Operacionais (CRO), abaixo ilustradas, das Sala Ae Sala B. 


51 MEDIÇÕES APÓS A CONCUSÃO DAS OBRAS E A INSTALAÇÃO DOS 
ESTUDIOS 

Depois de concluídas as obras e as instalações de ambos os estúdios, procedeu- 
se a medições acústicas visando a aplicação em produção fonográfica, ou seja, tomando 
como ponto de medição o sweet spot. 

Observando os resultados do Nível de Pressão Sonora em função da frequência 
em bandas de 1/3 de oitava, apresentadas a seguir para cada uma das salas, pergunta- 
se: - Como uma mesma sala pode dar respostas distintas, mesmo algumas opostas em 
relação a outras, ou ainda salientar bandas de frequência que estavam em desvantagem 
no diagrama de Bonello? Esta questão será abordada ao final deste item, após a exposição 
dos resultados. 
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Curva da Resposta Operacional - Sala A 
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Figura 25 - CRO da Sala A: Nível de Pressão Sonora em função da Frequência, medido no sweet spot. 
Monitores em configuração 5.1: L e R, em 120º e 135º, com e sem interferência acústica do monitor de 
vídeo. 


Observa-se que as curvas de respostas operacionais extrapolam, em alguns pontos, 
os limites de tolerância. Conforme a EBU Tech 3276 (EUROPEAN BROADCASTING 
UNION, 1998) obter todos os pontos dentro da margem de tolerância é bastante complexo. 
É comum o uso de equalizadores para pequenas correções em salas de controle. Na 
mesma figura são apresentadas as linhas que correspondem aos limites de tolerância 
indicados na referida recomendação. 

Os resultados das CRO da Sala À para os monitores de áudio frontais esquerdo e 
direito podem ser vistos no gráfico da Figura 26. Ambos apresentam respostas semelhantes 
acima de 160 Hz devido à geometria simétrica da sala. A curva da resposta do monitor 
esquerdo apresenta uma queda de 3 dB em relação à tolerância em 125 Hz, uma banda na 
qual temos M034 114,94 Hz, MOS7 119,53 Hz, MO42 123,23 Hz, M0O45 129,26 Hz e MO47, 
131,62 Hz (Figura 27, Figura 28), todos atuando ao mesmo tempo no sweet spot ou bem 
perto deste. Para o monitor direito a resposta foi um pouco mais plana nessa região, porém 
houve uma queda na CRO abaixo de 63 Hz apesar do TR mais elevado nessa mesma 
banda, coincidindo com M008 60,81 Hz (si % de tom abaixo) atuando no sweet spot. Nas 


frequências acima de 160 Hz o resultado está dentro da faixa recomendada. 
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Curva da Resposta Operacional - Monitores Frontais 
Esquerdo e Direito 
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Figura 26 — CRO da Sala A: Nível de Pressão Sonora em função da Frequência, medido no sweet spot 
para monitores na posição L e R. 





Figura 28 — Sala A: (esq.) MO47 : 131,62 Hz; (dir.) MO51 : 135,18 Hz. 
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Tempo de Reverberação - Sala A 
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Figura 29 — Tempo de reverberação medido na Sala A em função da frequência em 1/3 de oitava. 


Observa-se que o Tempo de Reverberação (TR) está dentro ou muito próximo da 
faixa de tolerância recomendada, com exceção da frequência de 80 Hz na curva do monitor 
de áudio central, posição na qual incidem nessa mesma banda de frequência os modos 
M01276,12 Hz, MO13 77,83 Hz, MO14 79,44Hz (Figura 30). O monitor de áudio esquerdo 
traseiro (L.) apresentou uma elevação no TR, em relação às demais curvas, em 160 Hz. 





Figura 30 — Sala A: (esq.) MO12 76,12 Hz, (centr.) MO13 77,83 Hz, (dir.) MO14 79,44Hz 
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Curva da Resposta Operacional - Sala B 
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Figura 31 - CRO da Sala B: Nível de Pressão Sonora em função da Frequência, medido no sweet spot. 


As CRO da Sala Be as curvas de tolerância recomendadas pela EBU Tech. 3276 
(EUROPEAN BROADCASTING UNION, 1998) são apresentadas na Figura 31. Observa-se 
um valor elevado do Nível de Pressão Sonora em 80 Hz na CRO para o monitor central, 
apesar do Diagrama de Bonello apontar para uma flutuação de contenção nessa banda 
específica [FLOAT 80] (Figura 8, Figura 10). Efetivamente nessa posição e nessa banda 
de frequência temos a incidência do segundo modo axial de profundidade P2, MO9 74,86 
Hz (si, 1/8 de tom acima) e do terceiro modo axial lateral LS, M1S 81,81 Hz (mi,1/8 de tom 
abaixo) (Figura 32). O monitor de vídeo altera significativamente o resultado da CRO, uma 
vez que parte do som direto não chega ao ouvinte. Nesse estúdio, o monitor de vídeo é de 
32 polegadas (-80 cm) e exerce grande impacto no resultado final da CRO para o monitor 


central (91< A... /4 < 121cm) podendo atuar para reflexões alimentando os lóbulos de 


80 Hz 
máximos e mínimos de M13. 
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Figura 32 — Sala B: (1) MO9 74,86 Hz; (2) M13 81,81 Hz. 
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Figura 33 - CRO da Sala B: Nível de Pressão Sonora em função da Frequência, medido no sweet spot 
para monitores na posição L e R. 


Na CRO da Sala B referente aos monitores frontais, esquerdo e direito, apresentada 


na Figura 33, observam-se alguns valores de Nível de Pressão Sonora fora da faixa de 
tolerância, porém na maior parte do âmbito de frequência audível o resultado está dentro 
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da faixa recomendada. Observa-se pouca variação acima de 400 Hz. Os estúdios de 
ambas as salas analisadas neste trabalho possuem um equalizador destinado a pequenos 
ajustes de frequência. Contudo tomamos a precaução de desligá-los durante as medições. 


Tempo de Reverberação - Sala B 
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Figura 34 - TR medido na Sala B em função da frequência em 1/3 de oitava. 


A Sala B possui bastante revestimento reflexivo em madeira na parte frontal, 
o que justifica o tempo de reverberação bem próximo do limite superior, e em algumas 
frequências ultrapassando a faixa de tolerância indicada pela EBU Tech 3276 (EUROPEAN 
BROADCASTING UNION, 1998). A região dos graves, onde se concentram as ondas 
estacionárias, apresentou resultado dentro do tempo de reverberação recomendado, com 
exceção do monitor direito em 63 Hz, próxima ao modo axial lateral L2, MO4 54,85 Hz (la,) 
(Figura 35) cujo lóbulo central incide nessa exata região. Nessa posição, o resultado obtido 
ultrapassou em 0,15 s a tolerância máxima do TR. 





Figura 35 — Sala B: MO4 54,85 Hz 
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Diversas interferências modais podem estar associadas a respostas irregulares ou 
inesperadas como as que destacamos neste item, pois quando uma fonte sonora se posiciona 
em uma área na qual há potencialmente um máximo ou mínimo de pressão, o respectivo 
modo associado a essa pressão se instala e a estacionária se estabelece. Quando houver 
mais de um modo com máximos ou mínimos na mesma região, há que se considerar as 
interferências de fase entre os diferentes modos, o que pode inclusive, conforme o caso, 
causar cancelamentos. Entenda-se assim como um pequeno deslocamento de uma fonte 
sonora pode causar uma mudança, até surpreendente, na resposta em frequência do sinal 
sonoro reproduzido pelos monitores de escuta. Isso explica como uma mesma sala pode 
estimular diferentes respostas para variados posicionamentos das fontes sonoras, neste 


caso, os monitores de áudio. 


61 CONCLUSÃO 


Apresentamos, discutimos e comentamos métodos aplicados para análise modal 
de dois projetos de estúdios de produção fonográfica, implementando novos recursos 
como a inclusão de notação musical para aperfeiçoar a representação de maneira a que 
se aproxime o máximo possível de imagens mentais da percepção sonora humana. As 
análises expostas foram discutidas, por um lado nos levando a algumas conclusões e por 
outro suscitando suposições e formulações de novas hipóteses a serem verificadas em 
trabalhos ulteriores. Destacam-se dentre essas conclusões e desdobramentos: (1) para a 
predição da resposta modal de uma sala, além da tradicional redução a uma sala retangular 
média (Figura 3, Figura 4), pode-se proceder a uma simulação por elementos finitos (item 
2. PROJETOS ANALISADOS E SIMULAÇÃO) centrando a atenção da análise nos modos 
axiais que se revelarem em função da conformação do contorno; (2) contornos irregulares 
podem influenciar a evolução dos modos axiais a partir de uma certa ordem, interferindo em 
sua progressão (Figura 24); (3) as posições dos monitores de áudio influenciam a resposta 
em frequência do sistema formado por eles mesmos integrados aos modos de vibração da 
sala (Figura 26, Figura 29, Figura 31, Figura 33, Figura 34). Para cada posição de monitores 
haverá uma resposta correspondente, porém todas remetendo à estrutura do modelo modal 
geral da sala; (4) Como elemento de design, além das proporções entre as dimensões 
da redução a uma sala retangular média, pode-se proceder ao exame das relações de 
proporção entre as frequências dos primeiros modos axiais que se revelarem na simulação 
(Figura 19, Figura 20); (5) Poderá ser feito um estudo ulterior sobre a hipótese de que 
determinadas relações intervalares de frequências estejam correlacionadas a constituição 
de determinados acordes, orientando a uma definição de proporções com resultados de 


distribuição modal mais equilibrada e homogênea. 
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CAPÍTULO 2 
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RESUMO: O desenvolvimento tecnológico tem 
progredido exponencialmente, dando origem a 
avanços como a inteligência artificial, robótica, 
impressoras 3D, nano e biotecnologias, fintechs, 
drones, e uma boa perspectiva futura para 
telefonia 5G e carros autônomos. Tais tecnologias 
viabilzam uma nova Revolução Industrial, que 
já está em andamento, ainda tímida e pouco 
utilizada pelas empresas do Brasil. Fato é que a 
Indústria 4.0 desponta como um caminho natural 
no ganho expressivo, não só de produtividade, 
como também de competitividade, por intermédio 
da integração das tecnologias. Por meio desta, 
observa-se uma maior produção utilizando menos 
recursos, além dacapacidade das empresas em 
fornecer produtos/serviços de alta qualidade, 
baixo preço, com menor tempo de entrega e 
ainda com uma boa rentabilidade. 
PALAVRAS-CHAVE: Industria 4.0, revolução 
industrial, mineração. 


HOW THE EXPRESSIONS DIGITAL 
TWINS, CYBER-PHYSICAL SYSTEMS, 
INTERNET OF THINGS AND CLOUD 
COMPUTING ARE INSERTED IN THIS 
INDUSTRIAL REVOLUTION? 
ABSTRACT: Technological development has 
progressed exponentially, giving rise to advances 
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such as artificial intelligence, robotics, 3D printers, 
nano and biotechnologies, fintechs, drones, and 
a good future perspective for 5G telephony and 
autonomous cars. Such technologies enable 
a new Industrial Revolution, which is already 
underway, still timid and little used by companies 
in Brazil. The fact is that Industry 4.0 emerges 
as a natural path towards expressive gains, not 
only in productivity, but also in competitiveness, 
through the integration of technologies. Through 
this, there is a greater production using fewer 
resources, in addition to the ability of companies 
to provide products/services of high quality, low 
price, with shorter delivery time and still with good 
profitability. 

KEYWORDS: Industry 4.0, Industrial Revolution, 
mining. 


Mas, afinal, o que vem a ser a Indústria 4.0? 


Como as expressões Gêmeos Digitais, 
Sistemas Ciberfísicos, Internet das Coisas e 
Computação em Nuvem são inseridas nessa 

Revolução Industrial? 


11 INTRODUÇÃO 


O desenvolvimento tecnológico tem 
progredido exponencialmente, dando origem a 
avanços como a inteligência artificial, robótica, 
impressoras 3D, nano e biotecnologias, 
fintechs, drones, e uma boa perspectiva futura 
para telefonia 5G e carros autônomos. Tais 


tecnologias viabilzam uma nova Revolução 


Capítulo 2 


Industrial, que já está em andamento, ainda tímida e pouco utilizada pelas empresas do 
Brasil. Fato é que a Indústria 4.0 desponta como um caminho natural no ganho expressivo, 
não só de produtividade, como também de competitividade, por intermédio da integração das 
tecnologias. Por meio desta, observa-se uma maior produção utilizando menos recursos, 
além dacapacidade das empresas em fornecer produtos/serviços de alta qualidade, baixo 
preço, com menor tempo de entrega e ainda com uma boa rentabilidade. 

A cada revolução industrial houve uma transformação radical no modelo produtivo, 
especialmente se tratando de máquinas cada vez mais protagonistas em relação ao 
trabalho humano. A Primeira Revolução foi caracterizada pela mecanização da produção 
através do uso de água e energia a vapor. Na Segunda Revolução, introduziu-se a produção 
em massa com a utilização da energia elétrica. Logo em seguida, a Terceira Revolução, 
também conhecida como Revolução Digital, consolidou o uso de aparelhos e dispositivos 
eletrônicos em uma maior automação da produção, principalmente com as tecnologias de 
inovação. Na indústria 4.0 tem-se a utilização de sistemas inteligentes, rápidos e precisos, 
com um elevado grau de automação e competência na tomada de decisões autônomas. 

Esta revolução considera a visão e execução de “fábricas inteligentes”, com 
estruturas modulares e sistemas ciberfísicos (CPS), com a união entre os mundos físico 
e virtual, conforme apresentado na Figura 1. Por meio da tecnologia de gêmeos digitais é 
criada uma representação virtual desse mundo físico com o auxílio da internet das coisas e 
seus dispositivos inteligentes, que comunicam e cooperam entre si, inclusive com humanos, 
em tempo real. Os serviços internos e intraorganizacionais são oferecidos e utilizados pelos 


participantes da cadeia de valor através da computação em nuvem. 
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Figura 1 — Constituição do sistema ciberfísico! 


O fato de que o modelo digital é constantemente atualizado com informações da 
parte física, faz com que o digital seja um modelo mais fiel do que está, de fato, acontecendo 
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no físico. Além disso, consegue armazenar todo o histórico de dados e prover otimização 
e predição para ele. 

Anteriormente, a escolha dos locais de instalação das indústrias era baseada em 
facilidade de escoamento de produção, disponibilidade de mão-de-obra e recursos e, 
ainda, disponibilidade do fornecimento de água e energia. Atualmente, mais um item torna- 
se imprescindível: conectividade à internet. 

O termo “Indústria 4.0” é originário de um projeto de alta tecnologia para 
informatização da manufatura do governo alemão, utilizado pela primeira vez na Hannover 
Messe, uma das maiores feiras de Inovação do mundo. O conjunto de recomendações 
para implantação desta nova revolução foi apresentado em outubro de 2012 e, em abril de 
2013, nesta mesma feira, apresentou-se o relatório final do grupo de trabalho dedicado a 


este tema. De acordo com este documento, os princípios desta revolução se baseiam em: 


* | Interoperabilidade: a habilidade dos sistemas ciberfísicos (tal como suporte 
de peças, estações de montagem e produtos), dos humanos e das fábricas 
inteligentes de se conectarem e se comunicarem entre si através da internet e 
da computação em nuvem; 


* Virtualização: uma cópia virtual das fábricas inteligentes é criada por sensores 
de dados interconectados, que monitoram processos físicos, com modelos de 
plantas virtuais e modelos de simulação; 


* Descentralização: a habilidade dos sistemas ciberfísicos das fábricas inteli- 
gentes de tomarem decisões sem intervenção humana; 


* || Acompanhamentos e Decisões em Tempo Real: a capacidade de coletar e 
analisar dados e entregar conhecimento derivado dessas análises imediata- 
mente; 


* Orientação a Serviço: oferta dos serviços dos sistemas ciberfísicos, humanos 
ou das indústrias inteligentes através da computação em nuvem; 


* Modularidade: adaptação flexível das fábricas inteligentes para requisitos mu- 
táveis através da reposição ou expansão de módulos individuais; e 


* || Transparência da informação: transparência de dados e informações para 
tomada de decisão, identificando o problema e áreas carentes de desenvolvi- 
mento, investimento ou tecnologia. 

Bloemet et al? citam três razões para implantação da Indústria 4.0 em uma empresa. 
A primeira é a comunicação machine-to-machine (M2M) onde, como no próprio nome diz, 
a comunicação é feita entre máquinas sem a necessidade de intervenção humana. Essa 
condição gera importantes contribuições no aumento da eficiência e segurança, além de 
permitir a conversão dos dados sem informações significativas, ou seja, os dados não 
apresentam influência humana e a capacidade de erro diminui, aumentando a eficácia na 


seguridade das informações. 
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A segunda razão comentada pelos autores é a possibilidade de realização de 
manutenção preditiva de máquinas e aparelhos, com base nos dados que eles enviam 
ao sistema. Com o uso excessivo, o equipamento pode ser danificado antes do tempo 
previsto e, por meio dos dados recebidos das atividades já executadas, a programação da 
manutenção pode ser feita de forma mais assertiva. A mudança da cultura de manutenção 
corretiva para preditiva evita ociosidade do equipamento, provocada por manutenções não 
planejadas, e permite que seja alcançado o pico de produtividade. 

Por fim, a última razão vai além de maximizar a eficiência operacional, pois tem como 
objetivouma maior interação com o cliente (engagement), já que é possívelcompartilhar 
informações de uso do produto vendido, por exemplo, de formaa criar novos valor e forma 
de prestar serviços. 


21 BRASIL NA INDUSTRIA 4.0 


Apesar da economia do Brasil representar cerca de 1/3 do PIB da América Latina, 
sendo uma das 10 maiores economias do mundo, este participa somente de 1,2% do 
comércio internacional. A participação da indústria de transformação no PIB, que já havia 
atingido mais de 20% em meados da década de 80, reduziu para 11% em 2016, fruto de 
mudança na estrutura produtiva do país. 

Além disso, o Brasil ocupa o 69º lugar no Ranking do Índice Global da Inovação, com 
33,1 pontos, onde são avaliados quesitos de crescimento de produtividade, investimentos 
em pesquisa e desenvolvimento, educação, exportação de produtos de alta tecnologia, 
dentre outros. Em primeiro lugar, tem-se a Suíça (67,7 pontos), seguida de Suécia (63,8), 
Países Baixos (63,3), EUA (61,4) e Reino Unido (60,8). 

O relatório Readiness for the Future of Production Report 2018 mostra o Brasil 
na 41º posição em termos da estrutura de produção e na 47º posição em relação aos 
vetores de produção da indústria. 

Com base nestes dados, conclui-se que a integração das tecnologias físicas e 
digitais, em todas as fases de desenvolvimento de um produto, ainda é pouco utilizado 
pelas empresas. O Brasil tem muitos campos de excelência, seja no agronegócio ou até 
mesmo em determinados ramos da indústria, mas ainda tem práticas pouco assertivas 
e voltadas ao tema da Industria 4.0. Em média, cerca de 43% das empresas brasileiras 
não identificam quais tecnologias tem potencial para alavancar a competitividade do setor 
industrial sendo que, nas pequenas empresas, esse percentual sobe para 57% e, nas 
grandes, recua para 32%). 

Segundo a Confederação Nacional da Industria (CNI), que publica pesquisas 
nacionais sobre adoção de tecnologias digitais relacionadas à era da manufatura avançada, 
a indústria brasileira ainda está se familiarizando com a digitalização e com os impactos 
sobre a competitividade. Acredita-se que o Brasil esteja pouco preparado para a adoção em 
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larga escala da Indústria 4.0, tendo em vista aspectos estruturais, educacionais e culturais. 

Reconhecendo a importância do tema, o Governo Federal lançou através do 
Ministério da Indústria, Comércio Exterior e Serviços (MDIC) e da Associação Brasileira de 
Desenvolvimento Industrial (ABDI), a Agenda Brasil para a Indústria 4.0, a qual possui um 
conjunto de iniciativas que visam promover o desenvolvimento deste tipo de Industria no 
país. Segundo levantamento da ABDI, a estimativa anual de redução de custos industriais 
no Brasil, a partir da migração da indústria para o conceito 4.0, será de, no mínimo, R$ 
73 bilhões/ano. As etapas para implantação da indústria 4.0 irá depender do grau de 
maturidade de cada empresa e, no geral, pode-se listar como um sequenciamento viável 


para a indústria brasileira, a implantação das medidas a seguir. 


Parte A —- Conhecimento — compreensão do conceito e avaliação 


|. Difusão do conteúdo, com conhecimento do que vem a ser a indústria 4.0, apoiada 
nos pilares da tecnologia, brasilidade e indústria. 


Il. Auto avaliação, para que seja identificada a maturidade da indústria em relação 
a jornada para implantação desta revolução. Geralmente, os steps de maturidade 
podem ser subdivididos em: 


a. Step 01 — Processos manuais, ausência de sistemas; 

b. Step 02 — Controle por código de barras, necessita de otimização de layout; 
c. Step 03 — Adoção de CAD, desenhada e parametrizada; e 

d. Step 04 —Sistemas integram simulações, visualização 3D da fábrica. 


III. Infraestrutura de banda larga para a indústria: redes de comunicação confiáveis, 
abrangentes e de alta qualidade. 


Parte B — Protótipo — quem são os parceiros tecnológicos e de negócios 


|. HUB 4.0 — plataforma de serviços para geração de valor, totalmente integrada a 
plataforma de auto avaliação, permitindo conexão aos provedores de tecnologia 
para digitalização e modernização dos processos. Parte do princípio de que 
as empresas, através de um levantamento de demanda, acionam consultores 
especializados que determinam uma agenda de implantação tecnológica, baseada 
em soluções aplicadas. 


||. Brasil mais produtivo (B+P), para diminuir o distanciamento tecnológico entre as 
grandes empresas com as micro e pequenas, ampliando o número de empresas 
que utilizam a digitalização industrial. 


III. Test beds, com a formatação de fábricas do futuro, em que é possível prototipar 
processos de implantação e experimentar as tecnologias, reduzindo as incertezas 
e OS riscos. 


IV. Conexão Startup-Indústria 4.0, com a massificação da conexão com foco nas 
necessidades da indústria nacional e utilização de propostas disruptivas e novos 
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modelos de negócio. O cadastro de startups pode ser feito através de workshops e 
imersões nas indústrias. 


V. Mining Hub, onde mineradoras apoiam o desenvolvimento de soluções que 
inovem o setor no Brasil, principalmente tratando-se de eficiência operacional, fontes 
alternativas de energia, gestão hídrica, gestão de rejeitos e resíduos, segurança e 
desenvolvimento social. 


Parte C — Requisito — definição do time e regras 


|. Talentos, com mapeamento de competências, entendimento das demandas de 
mercado, requalificação de trabalhadores e preparação das novas gerações. 


Il. Segurança, com a implantação de arquiteturas de segurança integradas e 
identificadores únicos, juntamente com os aprimoramentos relevantes para 
treinamento e conteúdo de desenvolvimento profissional contínuo. 


Il. Realização de estudos, para identificação de pontos estratégicos onde a 
tecnologia da Indústria 4.0 pode se iniciar na empresa. 


IV. Regras do jogo, com reformas legais e infra-legais para promoção da aceleração 
da jornada da indústria brasileira. Como exemplo, pode-se citar regulamentações de 
robôs colaborativos (COBOT), privacidade e proteção de dados e polos industriais. 


Parte D - Incentivos - para elevação da eficiência, produtividade e 
competitividade 


|. Financiabilidade, através de parcerias com bancos públicos e privados, para 
garantia de um leque acessível a diferentes empresas e necessidades. 


Parte E — Alianças estratégicas - competição igualitária com mercados 
externos 


||. Conexão global, com propostas de zeroing (zeragem) das alíquotas de importação 
de bens e insumos estratégicos e redução do imposto de importação para novas 
tecnologias. Vale também ressaltar a importância de cooperações e projetos 
bilaterais com demais países. 

Dentre 23 setores industriais do país, 14 estão em situação mais vulnerável, tais 
como setores de impressão e reprodução, farmoquímicos e farmacêuticos, químicos, 
minerais não-metálicos, couro e calçados, vestuário e assessórios, têxteis, máquinas e 
aparelhos elétricos, outros equipamentos de transporte, produtos de metal, máquinas 
e equipamentos, móveis, artigos de borracha e plástico, e produtos diversos. Em geral, 
estes segmentos apresentam produtividade inferior à média nacional e baixa inserção no 
comércio exterior, o que irá exigir um maior esforço para a implantação da Industria 4.0. 

A indústria extrativa (mineração - metálicos e rochas), metalurgia, companhias 
produtoras de energia elétrica, alimentícias, fumo, bebidas, celulose e papel, madeira e 
petrolífero, são exemplos de setores industriais com maior potencial para serem os líderes 
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na adoção de tecnologias, com maior proporção relativa de produtividade e altas taxas de 


exportação. 


31 MINERAÇÃO BRASILEIRA E A INDÚSTRIA 4.0 


Apesar da tendência de que a indústria extrativa tenha maior capacidade para a 
implantação da Indústria 4.0, e necessitam migrar para se manterem competitivos, tem- 
se que a mineração ainda está nos primeiros passos dos processos de modernização. 
Um dos principais motivos é o perfil conservador da indústria, onde também percebe-se a 
pequena integração entre os elos da cadeia produtiva. Segundo Moilanent), este cenário 
é denominado de “Sindrome da Tecnologia Comprovada”, onde as empresas querem 
se beneficiar das tecnologias mas, pelo fato de que ainda são recentes e sem provas 
concretas de sua eficácia em um horizonte temporal maior, causam insegurança. Assim, os 
empreendedores ficam receosos em testá-las e preferem utilizar métodos tradicionais, uma 
vez que a maior parte dos riscos já são conhecidos ou foram experimentados e vivenciados. 

Entretanto, estas novas tecnologias trazem inúmeras oportunidades para agregar 
valor aos produtos oferecidos aos clientes. A automação dos processos otimiza as 
operações de mineração e as mineradoras estão começando a perceber que inovação é 
uma solução estrutural de longo prazo ?. Na Tabela 1, tem-se alguns exemplos de como 
as tecnologias podem ser aplicadas na indústria da mineração, com algumas vantagens 
elencadas. 


USO NA MINERAÇÃO VANTAGENS 


- Segurança dos veículos e dos funcionários 

- Maior tempo de vida útil do equipamento 

- Maior eficiência produtiva 

- Maior velocidade média de transporte e, 
consequentemente, menor tempo médio de entrega 


Sistemas de Transportes Autônomos, como 
automação de trens e caminhões!) 


Sistema que detecta a sonolência e emite 
sons visuais e sonoros ao funcionário e à 
central de monitoramento! 


- Evita acidentes 
- Maior segurança e bem-estar dos funcionários 


Algoritmos traçam tendências e 
identificam quais projetos em P&D 
são as melhores oportunidades de 

investimento, principalmente relacionadas 

a sustentabilidade e energia!'?) 


- Detecção de focos de investimento que causem maior 
impacto positivo 


Detecção de corpos de minério em áreas 
novas/remotas, mesmo sem acesso a 
dados geológicos pré-existentes!! 


- Redução no tempo gasto em pesquisas preliminares 
- Melhor capacidade de planejamento da lavra 


- Monitoramento do processo de beneficiamento do 
Geometalurgia automatizada? minério para verificação se está operando conforme 
planejado 
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USO NA MINERAÇÃO VANTAGENS 


Separação de material de forma 
automatizada como, por exemplo, sistemas 
que usam jato de ar comprimido ao invés 
de água(!3:(14) 


Carregadeiras e Caminhões Articulados 
(LoadAssist e HaulAssist)''? 


Construção de banco de dados robusto e 
confiável para identificação de possíveis 
melhorias e controle da operação!'9 


Avaliação do potencial de jazidas por meio 
da comparação com o desempenho de 
áreas e depósitos mineiras similares, além 
de informações geoquímicas, geológicas e 
atividades sísmicas!'? 


Gêmeo digital interativo da mina que 
mostra o estado atual e o planejamento 
futuro do sítio? 


Sistema que detecta mineralizações para 
localização de novos empreendimentos 
viáveis!” 


Sistema de monitoramento das paredes de 
minas a céu aberto! 


- Menor quantidade de finos transportados para as 
barragens de rejeitos 

- Mais sustentável, em termos ambientais 

- Melhor eficiência energética (redução do consumo 
hídrico) 

- Mina teve uma taxa de recuperação global 20% maior 
do que o período em que não usava 


- Redução da variância da quantidade de material a ser 
transportado 

- Auxílio na gestão da produção e da frota de transporte 
- Otimização da produção 

- Redução do desgaste dos caminhões 


- Identificação de melhores práticas das unidades e 
replicação em outras unidades 

- Ganho de eficiência 

- Maior produtividade 


- Predição do tempo de vida útil da jazida 
- Recomendação das frentes de exploração 


- Auxilia na tomada de decisões sobre o futuro da 
exploração 


- Dados do monitoramento são utilizados para tomar 
decisões sobre evacuação de pessoas e equipamentos 
- Melhora a segurança da mina 





Tabela 1 — Exemplos de aplicação tecnológica na mineração 


Com base nestes exemplos de algumas aplicações já implantadas é possível 


vislumbrar a redução do impacto e dano ambiental nas atividade minerárias, tais como 


retirada de menos material, geração de menos resíduos e redução com custos de 


combustível e de manutenção de veículos. Ainda, maior segurança dos trabalhadores, 


aprimoramento da eficiência da produção, principalmente no que tange a melhoria da 


previsibilidade do produto, além da melhora na capacidade de planejamento da lavral'?. 


A redução de tempos ociosos não previstos na produção também deve ser 


considerada, uma vez que a simulação virtual permite que sejam verificados gargalos e 


falhas de implantação e operação. A produção deve ser planejada para que sejam mantidos, 


ou reduzidos, os custos e as tecnologias sejam especializadas. 


Vale ressaltar que, sem o enfoque adequado, podem ser desperdiçados grandes 
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investimentos para poucos resultados. Um ponto imprescindível é que a estratégia de 
digitalização da mina deve apresentar objetivo definido para a organização como um 
todo, ao invés de inovar em pontos desconectados entre si. O passo inicial é melhorar 
a capacidade de transformar dados em insumos importantes para resolver a gama de 
problemas enfrentados pelo setor, principalmente em relação à integração dos sistemas 
operacionais. 

Acredita-se que com a evolução e o barateamento das tecnologias necessárias 
para colocar o programa Indústria 4.0 em prática, haverá a transição de um maior número 
de indústrias do nível de automação atual para a utilização de sistemas ciberfísicos. 
Entretanto, é importante atentar-se que, para os setores em que o Brasil, hoje, se encontra 
em situação relativamente confortável em função das suas vantagens comparativas, O 
avanço tecnológico em outros países poderá provocar grandes pressões competitivas no 


futuro. 


41 ACELERAÇÃO DA INDÚSTRIA 4.0 NO BRASIL 


A pandemia do coronavírus gerou impactos profundos na economia mundial e no 
comportamento das pessoas. Apesar das grandes perdas, a crise forçou uma aceleração 
tecnológica em todo o mundo e as tendências de implantação da indústria 4.0 no 
Brasil indicam certo progresso. A necessidade de que profissionais continuassem suas 
atividades em regime home-office evidenciou, claramente, como processos precisam ser 
compartilhados, integrados e descentralizados. Coneciar-se é mais do que fundamental. 

No que diz respeito ao consumo, a simulação de cenários para planejamento de 
demanda produtiva foi bastante utilizada. Porém, ocorreu de forma limitada, uma vez 
que as empresas não conseguiram prever o volume de saída de muitos itens, gerando 
carências no público consumidor. No Brasil a pandemia não apresentou grandes cadeias 
de desabastecimento mas, ainda, uma grande redução produtiva, principalmente de micro 
e pequenas empresas. 

Os sistemas tecnológicos, geralmente, se baseiam em dados históricos para que 
sejam feitas previsões. Ao se deparar com um cenário novo e, no caso desta recente 
pandemia, nunca visto nesta escala, os cálculos não se consolidaram. Por isso a tendência 
da indústria em buscar modelos mais inteligentes para fazer simulações, a fim de prever 
cenários menos previsíveis, deve estar fundamentalmente em crescimento. 

Além da pandemia, outra fonte catalisadora da Indústria 4.0 é a rede 5G, já que a 
ethernet ainda predomina na infraestrutura de muitas indústrias, com um custo elevado. 
O avanço da Indústria 4.0 ganhará ainda mais força com a disseminação desta nova 
tecnologia no país, uma vez que trará muitas vantagens operacionais para as fábricas 
inteligentes. Dentre elas, pode-se citar maiores qualidade e velocidade de comunicação 


para os dispositivos sem fio ganharão, o que permitirá que os gestores utilizem o máximo 
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do potencial da Internet das Coisas, automatizando uma infinidade de processos. Apenas 
para se ter uma ideia, a 5G permite downloads cerca de 82 vezes mais ágeis que o 4G, 
prometendo uma base sólida para o desenvolvimento das indústrias do futuro. 


5 | CONSIDERAÇÕES FINAIS 


É importante ressaltar que ainda caminha-se de forma bem conservadora para a 
implantação da Indústria 4.0 no Brasil. É imprescindível rever as bases desta implantação 
e adaptar o que já se tem para o que realmente é relevante. 

A pandemia do coronavírus acelerou a disseminação da Indústria 4.0 em setores 
farmacêuticos, vendas, bens de consumo não duráveis, mineração e alimentação, uma vez 
que estes tiveram que acelerar suas produções e inovar com a ajuda das tecnologias. É 
certo que, se antes as empresas já não estavam repensando suas estratégias operacionais 
e adotando as tecnologias da Indústria 4.0, a possibilidade de crise força a busca pela 
inovação para que se mantenham competitivas e bem-posicionadas no mercado. 

Neste contexto o velocidade da comunicação e o diminuição do preço das tecnologias 
5G, além do alcance das mesmas em varias regiões do país, descentralizará o auto custo. 
Esse passo é importante para que a industria consiga alavancar e manter seus produtos 
competitivos, abrindo um leque de oportunidade, em âmbito global, haja visto que o Brasil 
é um país vasto territorialmente e com variadas riquezas. 

Embora exista esse movimento global de desenvolvimento das fábricas inteligentes 
e o avanço das comunicações digitais, uma nova tendência já tem sido discutida e 
apontada como uma possível Indústria 5.0. Com o propósito de abordar a temática de 
indústrias colaborativas, insere, novamente, o toque humano na produção. Ao contrário 
das revoluções anteriores que já tinham a parcela humana no sistema produtivo, nesta, há 
uma redistribuição da criatividade, necessária ao atendimento das exigências de padrões 
consumidores personalizados. As máquinas são excelentes para fabricação de produtos 
padronizados mas, sozinha, a tecnologia não consegue abordar essa exclusividade 
nos produtos. No contexto desta nova possibilidade, os robôs assumem tarefas árduas 
e perigosas e não se destinam a substituição da força de trabalho. A sustentabilidade, 
abertura e inclusão são os pilares defendidos por esta proposta, com uma tecnologia 
centrada na humanidade e auxiliando a qualidade de vida. 
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RESUMO: Embora as maiores reservas de 
carvão do Brasil estejam situadas no Estado 
do Rio Grande do Sul, a Bacia Carbonífera 
Sul-Catarinense, localizada no Estado de 
Santa Catarina, concentra o maior número 
de minas de carvão em subsolo do País. O 
conhecimento das características mecânicas 
das rochas que compõem a camada de carvão 
e as camadas encaixantes é de fundamental 
importância para o dimensionamento das 
estruturas de sustentação das minas (pilares) e 
suporte do teto das escavações. Este trabalho 
reporta os estudos realizados pela CIENTEC 
- Fundação de Ciência e Tecnologia do Estado 
do RS sobre rochas sedimentares presentes 
nas camadas encaixantes da camada de 
carvão Barro Branco. Os trabalhos de campo 
consistiram na execução de sondagens rotativas 
com diâmetro NX utilizando amostrador duplo- 
móvel. Em laboratório foram realizadas análises 
petrográficas macrosocópicas, ensaios para 
a determinação dos índices físicos e ensaios 
mecânicos para a determinação da resistência 
e deformabilidade dos diferentes litotipos. Os 
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resultados permitiram a caracterização das 
rochas encaixantes da camada de carvão Barro 
Branco através da definição de padrões típicos 
de resistência. Esses padrões são usados para 
a análise da capacidade de carga das rochas 
encaixantes no processo de dimensionamento 
dos pilares de sustentação das escavações 
subterrâneas, quando utilizado o método de lavra 
por câmaras e pilares. 

PALAVRAS-CHAVE: Formação Rio Bonito, 
carvão Barro Branco, ensaios de laboratório, 
caracterização geomecáânica. 


GEOTECHNICAL FEATURES OF 
THE HOST ROCKS OF BARRO 
BRANCO COAL SEAM OF THE SOUTH 
CATARINENSE COALFIELD 


ABSTRACT: Although the largest coal reserves 
in Brazil are located in the state of Rio Grande 
do Sul, the South-Catarinense Coalfield located 
in the state of Santa Catarina concentrates the 
largest number of coal underground mines in the 
country. This paper report studies conducted by 
CIENTEC - Fundação de Ciência e Tecnologia 
of the RS State on sedimentary rocks present 
in the enclosing layers of the Barro Branco coal 
seam. The samples were collected through 
vertical boring, using a double tube core barrel 
(NX diameter). The cylindrical samples were 
submitted to macroscopic petrographic analysis, 
physical and mechanical laboratory testis to 
determine the strength and deformability of 
different rock materials. The results allowed the 
characterization of the different sedimentary rock 
types enclosing the Barro Branco coal seam. 
These resistance patterns are used for the 
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analysis of the carrying capacity of the host rocks and design of the underground excavation 
pillars, when the room and pillar mining method is used. 

KEYWORDS: Rio Bonito Formation, Barro Branco coal seam, laboratory tests, geomecanical 
characterization. 


11 INTRODUÇÃO 


Poucos estudos geomecânicos são conhecidos sobre as rochas sedimentares 
do Sul do Brasil, especialmente sobre as rochas do pacote sedimentar encaixante das 
camadas de carvão. Normalmente os estudos existentes estão associados a materiais 
presentes em projetos de engenharia civil (DOBEREINER, 1987; DUARTE et al., 2010). 
No que tange às rochas encaixantes das camadas de carvão, destaque especial deve ser 
dado aos trabalhos realizados pelo IPT (1987), sobre amostras de arenitos laminados, 
siltitos e siltitos laminados do estado de Santa Catarina, e estudos detalhados realizados 
por ZORZI et al. (1989) com mais de uma centena de amostras de siltitos, arenitos, marga 
e pedra areia das Formações Palermo e Rio Bonito presentes na Bacia Carbonífera Santa 
Terezinha, no Estado do Rio Grande do Sul. 

A partir do início da década de 1990, a CIENTEC - Fundação de Ciência e Tecnologia 
desenvolveu estudos detalhados em laboratório para a caracterização geomecânica dos 
materiais presentes nas camadas encaixantes das principais camadas de carvão no Sul do 
Brasil (ZORZI, 1990; 1996). Na Bacia Carbonífera Sul-Catarinense esse trabalho iniciou com 
as rochas encaixantes da camada Barro Branco, a principal camada de carvão minerada 
em subsolo, seguindo-se depois nas encaixantes das camadas de carvão Irapuá e Bonito, 
nessa mesma bacia. O objetivo dos ensaios em laboratório foi definir o comportamento 
mecânico das rochas encaixantes do teto e piso imediatos da camada de carvão. 

Na mineração de carvão em subsolo no Brasil, a geometria das galerias consiste 
em aberturas com 5,5 a 6m de largura e altura correspondente a espessura minerada 
da camada de carvão, da ordem de 1,8 a 2,3m. Dentro desse contexto, os primeiros 
5m acima e abaixo da camada de carvão constituem a parcela do maciço rochoso que 
sofre a maior redistribuição das tensões induzidas em função da escavação. O principal 
sistema de suporte das galerias é o de parafusos de teto (roof bolting). Nesse sistema, os 
mecanismos mais utilizados são o de efeito supensão (supension effect) e o de efeito viga 
(friction or beam-building effect) (PENG, 1986). Conhecer as características mecânicas 
e a distribuição espacial dos diferentes litotipos presentes na encaixante do teto é de 
fundamental importância para o projeto do sistema de suporte das galerias. 

De outra parte, quando tratamos do dimensionamento dos pilares de carvão, 
não deve ser esquecido o papel desempenhado pelas rochas encaixantes, como 
fundação responsável pelo suporte das cargas transmitidas pelos pilares. O método de 
dimensionamento de pilares de carvão recomendado pelo DNPM - Departamento Nacional 
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da Produção Mineral atual ANM - Agência Nacional de Mineração, contempla, na sua 
concepção, a análise da capacidade de carga das rochas encaixantes (ZORZI et al. 1991). 
Portanto, conhecer esse atributo das rochas encaixantes, em especial, do piso é condição 
básica para o projeto dos pilares das minas de carvão brasileiras. 

Neste capítulo, são apresentados os procedimentos adotados para a caracterização 
geomecânica das rochas encaixantes da camada de carvão Barro Branco, minerada na 
bacia Carbonifera Sul-Catarinense, no Sul do Brasil. Os procedimentos adotados in situ 
para a coleta de amostras, as técnicas de preparação das amostras e os procedimentos 
utilizados para a execução dos ensaios de caracterização física e mecânica, em laboratório, 
são detalhadas. Os resultados obtidos foram tratados estatisticamente e permitiram a 
tipificação das principais litologias presentes na Bacia Carbonífera. Essas informações 
servem de base para o projeto do sistema de suporte do teto e cálculo da capacidade de 
carga das rochas encaixantes do piso da camada de carvão Barro Branco. 


2 | GEOLOGIA DA ÁREA DE ESTUDO 


Geograficamente, todos os locais investigados estão localizados dentro da Bacia 
Carbonífera Sul-Catarinense, a qual é parte da porção Leste da bacia do Paraná, conforme 
Figura 1, e se estende por vários Estados do Brasil. Os depósitos sedimentares mais 
antigos pertencem ao período Devoniano. A evolução da bacia possui características 
intracratônicas. Dentro da sequência estratigráfica, as litologias carbonosas estão 
associadas ao Grupo Itararé (basal) e Guatá, com idades do Carbonífero Superior e 
Permiano Superior. Porém, no Grupo Guatá estão localizadas as camadas de carvão mais 
espessas, de melhor qualidade e em explotação. Dentre as sete principais camadas, três 
possuem maior expressão econômica: Barro Branco, Irapuá e Bonito distribuídas na coluna 
estratigráfica como visto na Figura 2. 
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Bacia Carbonifera Sul-Catarinense 


==] Sedimentos do Quaternário 
[++] Derrames de Basalto 
[o] Sedimentos Gonduânicos 
Sedimentos Devonianos 
Escudo Pré-Devoniano 


Figura 1 - Mapa geológico do Sul do Brasil. Números de | a XIll indicam bacias carboníferas 
conhecidas no Brasil. IX = Bacia Carbonífera Sul-Catarinense (Adaptado de ZORZI, 1990). 


O condicionamento estrutural encontrado na bacia do Paraná é representado por 
falhas normais, com planos subverticais. Falhas de empurrão são raras e com pequenos 


rejeitos. A presença de slikensides causados pela compactação diferencial é frequente. 


Sistemas de fraturas podem acompanhar as falhas, formando sistemas conjugados, 


preenchidos ou não por diques de diabásio. 
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Figura 2 - Perfil estratigráfico típico da Formação Rio Bonito, com destaque para a camada Barro 
Branco (adaptado de COSTA, 2000). 


O piso imediato da camada de carvão Barro Branco é formado, na maior parte, por 
arenitos de granulação fina à grossa e composição quarizosa. Estes arenitos apresentam, 
em alguns casos, uma laminação ondulada à lenticular, dada pela presença de finas 
lâminas de material síltico-argiloso. Por vezes podem aparecer camadas de material síltico 
somente, às vezes físsil, intercalado ao pacote de arenitos. Nas proximidades do piso 
da camada Barro Branco, localmente, pode ocorrer camada de siltito preto, compacto, 
fortemente alterável na presença de água. 

Já no teto da camada Barro Branco podem aparecer somente arenitos, somente 
siltitos ou uma sucessão de camadas de arenitos e siltitos. O siltito, que compõe o teto 
imediato em alguns dos locais amostrados, apresenta variações desde siltito cinza escuro 
ou preto, maciço e consistente, até uma rocha físsil e mais micácea (folhelho síltico), ou 
então laminado, com lâminas formadas por arenito fino. 

As demais camadas de siltito que compõem a porção mais superior do teto imediato 
são normalmente arenosas, com presença de micas e com lentes de arenito fino. A presença 
de pirita é constante, com proporções variáveis e ocorrendo sob a forma de plaquetas 
paralelas à fissibilidade e laminação, ou preenchendo fissuras, com ou sem carbonato 
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associado. Os arenitos, por sua vez, são rochas de granulometria fina, as vezes grosseira, 
quartzosos, pouco micáceos e coerentes. Apresentam contatos nítidos ou gradacionais 
com os siltitos de base e topo. Nas camadas mais distantes do teto imediato os arenitos se 
mostram mais micáceos e com lâminas síltico-argilosas. 

A seleção dos locais amostrados dentro da Bacia Carbonífera Sul-Catarinense levou 
em consideração estudos preliminares da geologia e a distribuição dos diferentes litotipos 
das encaixantes das camadas de carvão na jazida. Adicionalmente, a distribuição dos locais 
de amostragem contemplou as minas em atividade mais importantes e representativas da 
bacia carbonífera. As variações litológicas dentro das minas definiram o número de furos 
realizados em cada mina ou dentro de uma mesma mina. A amostragem do teto e do piso 
da camada de carvão Barro Branco abrangeu as principais minas em operação no final da 
década de 1980. Nesse contexto espacial, foram coletadas amostras de 4 minas importantes 
a época da execução dos trabalhos de amostragem (1988), quais sejam: Verdinho (Site 1), 
A-Sangão (Site 2), Santana (Site 3) e 3G (Site 4), conforme sites designados na Figura 3. 
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Figura 3 - Localização geográfica dos pontos de amostragem (Sites 1, 2, 3 e 4) das rochas encaixantes 


das camadas de carvão Barro Branco na Bacia Carbonífera Sul-Catarinense (adaptado de ZORZI, 


1990). 


31 AMOSTRAGEM 


As sondagens foram realizadas em subsolo a partir do acesso direto em galerias 


nas frentes de lavra. Esse procedimento procurou buscar amostras numa condição mais 


preservada dos efeitos de desmonte e tensões induzidas in situ. 


No total foram realizados 20 furos verticais, amostrando-se em cada local 5m no 


piso e 5m no teto, distribuídos em 4 minas localizadas nos 4 sites indicados na Figura 3, 


totalizando 100m de testemunhos coletados. 


A coleta das amostras foi realizada com o emprego de sonda rotativa, utilizando 


amostrador duplo móvel com diâmetro NX (54mm). A qualidade da sondagem proporcionou 


recuperação acima de 95%. Imediatamente após a realização das sondagens as amostras 
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foram embaladas em papel alumínio e/ou filme plástico, como forma de preservar a 
qualidade do material frentes aos agentes intempéricos, especialmente processos de 
oxidação devido à presença de pirita em algumas litologias. 

Nas Fotos 1 e 2 podem ser vistos testemunhos típicos das litologias presentes no 
site 1, região em que são observadas as principais litologias presentes nas demais minas 
na Bacia Carbonífera Sul-Catarinense, enquanto na Figura 4 são apresentados os perfis 
litológicos de todos os furos de sondagem realizados nos diferentes locais de amostragem. 


4 | ENSAIOS DE LABORATÓRIO 


Previamente à seleção das amostras para ensaio em laboratório as diferentes 
litologias foram descritas macroscopicamente e determinado o índice de qualidade da rocha, 
RQD - Rock Quality Designation. O Quadro 1 apresenta uma síntese das características 
petrográficas macroscópicas das diferentes classes geomecânicas identificadas nas 
amostras estudadas. 

A preparação dos corpos de prova obedeceu aos procedimentos específicos 
sugeridos, para cada modalidade de ensaio, pela ISRM (1978a; 1978b; 1979a; 1979b; 
1983). Previamente aos ensaios mecânicos, foram determinados a massa específica 
aparente e velocidade sônica (onda P), esta última determinada utilizando um equipamento 
“PUNDIT” (Portable Ultrasonic Non-Destructive Digital Indicating Tester) com transdutores 
de 54 KHz, 1 MHz and 2 MHz. Com o intuito de determinar a envoltória de ruptura das 
diferentes litologias, foram realizados ensaios de compressão diametral, compressão 
uniaxial e compressão triaxial. 

No total, 702 corpos de prova foram ensaiados, conforme quantitativos vistos 
no Quadro 2. Os quantitativos de massa específica e de velocidade sônica podem ser 
considerados também a expressão aproximada da proporção de cada classe geomecânica 
na jazida. 

Os ensaios de compressão diametral foram realizados em uma prensa Losenhausen 
com capacidade de carga de 100kN. A tensão de tração indireta foi determinada na direção 
paralela aos planos de estratificação da rocha. 

Os ensaios de compressão uniaxial e compressão triaxial foram realizados utilizando 
dois sistemas de carga: - um composto por um conjunto pórtico de reação-macaco 
hidráulico, com capacidade de carga de 2MN e; - outro consistindo de uma prensa EMIC, 
com capacidade de carga de 600kN. Em ambas as modalidades de ensaio a aplicação da 
tensão principal maior foi na direção perpendicular aos planos de estratificação da rocha. 
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Foto 1 - Testemunhos de sondagem diâmetro NX (54 mm) obtidos nas sondagens das rochas 
encaixantes do teto e do piso da camada de carvão Barro Branco - Local 1 na Mina Verdinho - Site 1. 
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Foto 2 - Testemunhos de sondagem diâmetro NX (54 mm) obtidos nas sondagens das rochas 
encaixantes do teto e do piso da camada de carvão Barro Branco - Local 3 na Mina Verdinho - Site 1. 
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Figura 4 - Perfis geológicos típicos das rochas encaixantes da camada de carvão Barro Branco das 
diferentes minas amostradas na Bacia Carbonífera Sul-Catarinense (adaptado de ZORZI, 1990). 


Classe 


Geomecânica 


Arenito HS 
(Alta 
Resistência) 


Arenito MS 
(Média 
Resistência) 


Arenito LS 
(Baixa 
Resistência) 


Descrição Macroscópica 


Arenito orto-quartzítico, de granulação fina a média, cor cinza claro com cimento 
carbonático. Zonas com manchas claras são resultado de forte cimentação 
carbonática. Possui estratificação plano-paralela ou cruzada de baixo ângulo, 
apresenta-se compacto com frequentes intercalações de filmes de siltito, matéria 
orgânica, micas e concentrações de pirita. A rocha é muito coerente com RQD 
de bom a excelente (>75%) nas zonas maciças e pobre (<25%) nas zonas com 
intercalações. 


Arenito de granulação fina a média, cor cinza claro, intercalado com finas lâminas 
ou leitos de siltito escuro, micáceo. Cimentação quartzítica. Fraturas frequentes de 
contato. RQD de excelente (>/5%) a muito pobre (<25%). 


Arenito de granulação grossa, cor cinza claro com manchas escuras. Intercalações 
de filmes de matéria orgânica dão origem a fraturas horizontais de acamamento 
frequentes. A cimentação é carbonática, embora a rocha possua aspecto friável. 
RQD pobre (25-50%). 
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Siltito cinza escuro, com finas lâminas de arenito intercalado. Fraturas muito 
frequentes nos planos de acamamento. RQD pobre (<50%). 
Siltito HS 
(Alta 
Resistência) 


Itito carbonoso intercalado com arenito fino formando pequenas lentes. Apresenta 


estratificação plano-paralela, extremamente fraturado nos planos de acamamento. 
RQD muito pobre (<25%). 
Siltito LS 


(Baixa 
Resistência) 





Quadro 1 - Características petrográficas macroscópicas típicas das diferentes classes geomecânicas 
das rochas encaixantes da camada de carvão Barro Branco na Bacia Carbonífera Sul-Catarinense. 


Nos ensaios triaxiais a tensão confinante foi aplicada com o auxílio de uma câmara 
triaxial de Hoek. A ruptura seguiu os procedimentos de ensaio individual e/ou multiestágio. 
Todas as litologias foram ensaiadas na faixa de tensões confinantes entre 10 e 30MPa. 


Modalidade de Ensaios 


Classe 
Geomecânica 


M Módulo de Young/ 

o Velocidade a J 
Específica Ago Coeficiente de 
Sônica 
Aparente Poisson 


Arenito LS 
Siltito LS 
Total de Ensaios 


Quadro 2 - Quantitativos de ensaios realizados com materiais das diferentes classes geomecânicas. 





Nos ensaios de compressão uniaxial com determinação das constantes elásticas, as 
deformações foram medidas com extensômetros elétricos tipo strain gages e um sistema 
eletrônico Kyowa. O cálculo do módulo de Young foi feito pelo método médio, enquanto o 
coeficiente de Poisson baseou-se no método secante 50%. 

Todos os ensaios mecânicos foram realizados baseados nos métodos sugeridos 
pela ISRM com as amostras na condição de umidade natural entre 0,5% e 3%. 

Para melhor visualização da dispersão dos resultados das diferentes modalidades 
de ensaios foram utilizados diagramas de caixas (Box and Whisker), conforme detalhado 
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no Anexo 1, com o auxílio do Software Statgraphics Plus, version 4.1. Os resultados das 
diferentes propriedades físicas e mecânicas são apresentados graficamente nas Figuras 
5a a 5f. Após o tratamento estatístico dos dados foi possível a definição de cinco classes 
geomecânicas (três classes de arenitos e duas de siltitos). O processo de classificação 
geomecânica dos diferentes tipos litológicos estudados foi feito com base na análise 
petrográfica macroscópica, seguida da análise estatística dos resultados da velocidade 
sônica, resistência à tração e resistência à compressão uniaxial. Finalmente, utilizaram- 
se os resultados de todas as modalidades de ensaios mecânicos (compressão diametral, 
uniaxial e triaxial), individualizados por classe geomecânica, para o estabelecimento da 
envoltória de ruptura média de cada classe. 

A envoltória de ruptura de cada uma das 5 classes foi definida com base no critério de 
ruptura de Hoek-Brown (HOEK et al., 2002). Os cálculos dos parâmetros representativos da 
envoltória da rocha intacta foram feitos com o auxílio do programa RocLab, version 1.038. 
A envoltória média de cada classe geomecânica levou em consideração todos os ensaios 
das diferentes modalidades de ensaio. Valores espúrios e/ou resultados de amostras pouco 
representativas da classe foram desconsiderados, adotando-se uma postura conservadora 
na definição da envoltória de ruptura. 
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Velocidade Sônica (km/s) 
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Módulo de Young (GPa) 


Arenito HS | Arenito MS Arenito LS siltito HS ailtito LS 
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Figura 5 - Características físicas - (a) massa específica aparente, (b) velocidade sônica e - 
propriedades mecânicas - (c) resistência a tração, (d) resistência à compressão uniaxial, (e) módulo de 
Young médio e (f) coeficiente de Poisson secante 50% - das rochas encaixantes da camada de carvão 

Barro Branco. 


As Figuras 6a a 10a apresentam os resultados dos ensaios das amostras ensaiadas 
em cada uma das 5 classes, juntamente com a respectiva envoltória de ruptura. Já 
nas Figuras 6b a 10b podem ser vistas as envoltórias das amostras juntamente com a 
envoltória média da respectiva classe geomecânica. Nesse conjunto de resultados estão 
representadas as 5 classes geomecânicas: 


*  3classes de Arenitos HS, MS e LS; 
*  2classes de Siltitos HS e LS. 


As envoltórias médias das 5 classes geomecânicas são apresentadas no Quadro 3 
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e graficamente representadas na Figura 11. 


5 | ANÁLISE DOS RESULTADOS 


Os estudos realizados permitiram a caracterização geomecânica das principais 
litologias presentes nas encaixantes da camada de carvão Barro Branco, dentro da Bacia 
Carbonífera Sul-Catarinense. Litologicamente, são identificadas duas litologias bem 
distintas: arenitos e siltitos. Sob o ponto de vista das propriedades físicas, podem ser 
observadas 5 classes, especialmente quando observamos a velocidade sônica, conforme 
Figura 5b. É notável a diferença da velocidade (mediana) entre as 3 classes de arenitos, 
assim como entre as duas classes de siltitos. Da mesma maneira, mecanicamente, tanto 
a resistência à tração como a de compressão evidenciam a presença de 5 classes de 
resistência, embora as classes geomecânicas Arenito LS e Siltito LS guardem boa 
semelhança, como visto nas Figuras 5c e 5d. Esse mesmo comportamento é observado 
no que se refere ao módulo de Young (Figura 5e), embora o Siltito LS apresente uma 
variabilidade muito maior nos seus resultados. 

Do ponto de vista geomecânico, os arenitos podem ser divididos em 3 classes: - 
Arenito HS, de grão fino, o mais resistente; - Arenito MS, de grão médio a fino, praticamente 
impossível de ser separado apenas pela descrição macroscópica; - Arenito LS, de grão 
médio a grosso, o menos resistente. Os siltitos, por seu lado, são divididos em duas classes: 
- Siltito HS e o Siltito LS. Macroscopicamente esses siltitos se distinguem apenas pelo 
grau de fissilidade, uma vez que pela textura e cor são muito semelhantes. Do ponto de 
vista prático, tendo como base as características petrográficas, físicas e as propriedades 
mecânicas das rochas estudadas, foram definidas 5 classes geomecânicas, conforme visto 
no Quadro 3. 

Os materiais pertencentes às classes geomecânicas Arenito HS and Siltito HS 
representam, aproximadamente, 50% e 20%, respectivamente, das litologias amostradas. 
Cada uma das 3 outras classes representa cerca de 10%. 

Do ponto de vista da resistência ao cisalhamento, por sua vez, a envoltória de 
ruptura das diferentes classes geomecânicas sugere a existência de apenas 4 classes com 
características distintas, sendo duas de arenitos e duas de siltitos. 
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Figura 6a - Amostras das encaixantes da camada de carvão Barro Branco classificadas 
geomecanicamente como Arenito HS. Resultados de ensaios De cada amostra e envoltória de ruptura 
da rocha intacta, de acordo com o critério de Hoek-Brown. 





—a— Amostra 1 
—- Amostra 2 
—— Amostra 3 
—e- Amostra 4 
—*— Amostra 5 
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-e- Amostra 7 
—e Amostra 8 
—e- Amostra 9 
—+- Amostra 10 


—- Amostra 11 


se Amostra 12 
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Tensão Principal Menor, o; (MPa) 


Parâmetros da envoltória de ruptura da rocha intacta (média de todas as amostras): 
Arenito HS: oc = 102,6MPa; m;=10,2;s=1;a = 0,5 c = 20,6MPa; é = 43,2º 





Figura 6b - Envoltórias de ruptura, de acordo com o critério de Hoek-Brown, das amostras de rocha 
intacta classificadas geomecanicamente como Arenito HS. 
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Figura 7a - Amostras das encaixantes da camada de carvão Barro Branco classificadas 
geomecanicamente como Arenito MS. Resultados de ensaios de cada amostra e envoltória de ruptura 
da rocha intacta, de acordo com o critério de Hoek-Brown. 
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Figura 7b - Envoltórias de ruptura, de acordo com o critério de Hoek-Brown, das amostras de rocha 
intacta classificadas geomecanicamente como Arenito MS. 
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Figura 8a - Amostras das encaixantes da camada de carvão Barro Branco classificadas 


da rocha intacta, de acordo com o critério de Hoek-Brown. 
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os resultados da amostra 1 são pouco representativos dentro dessa classe geomecânica, por isso a pequena 
influência na definição da envoltória média de resistência. 


Parâmetros da envoltória de ruptura da rocha intacta (média de todas as amostras): 


Arenito LS: o; = 58,1MPa; m,=7,8;s=1;a=0,5c' = 12,4MPa; dé = 40,2º 





Figura 8b - Envoltórias de ruptura, de acordo com o critério de Hoek-Brown, das amostras de rocha 


intacta classificadas geomecanicamente como Arenito LS 
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Figura 9a - Amostras das encaixantes da camada de carvão Barro Branco classificadas 
geomecanicamente como Siltito HS. Resultados de ensaios de cada amostra e envoltória de ruptura da 
rocha intacta, de acordo com o critério de Hoek-Brown. 
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Figura 9b - Envoltórias de ruptura, de acordo com o critério de Hoek-Brown, das amostras de rocha 
intacta classificadas geomecanicamente como Siltito HS. 
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Figura 10a - Amostras das encaixantes da camada de carvão Barro Branco classificadas 
geomecanicamente como Siltito LS. Resultados de ensaios de cada amostra e envoltória de ruptura da 
rocha intacta, de acordo com o critério de Hoek-Brown. 
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Parâmetros da envoltória de ruptura da rocha intacta (média de todas as amostras): 
Siltito LS: oa = 32,7MPa;m;=4,1;s=1;a=0,5c' = 8,4MPa; 6 = 32,5º 





Figura 10b - Envoltórias de ruptura, de acordo com o critério de Hoek-Brown, das amostras de rocha 
intacta classificadas geomecanicamente como Siltito LS. 
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Propriedades Típicas da Rocha Intacta a Partir de Ensaios de Laboratório 


Parâmetros do Critério de 
Ruptura da Rocha Intacta 
Classe ess e Br 
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* Valores entre parênteses = desvio padrão. 


r=massa específica aparente natural; Vp=velocidade da onda P; o =resistência à tração indireta; 
o =resistência à compressão uniaxial; E =módulo de Young médio; n....o=Coeficiente de Poisson 
método secante 50%; s =resistência à compressão uniaxial, m, s, a=parâmetros do critério de ruptura 


de Hoek-Brown; c'=coesão, f'=ângulo de atrito do critério de ruptura de Mohr-Coulomb. 


Quadro 3 - Parâmetros médios das classes geomecânicas definidas para as rochas encaixantes da 
camada de carvão Barro Branco na Bacia Carbonífera Sul-Catarinense. 
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Figure 11 - Envoltórias de ruptura média, de acordo com o critério de ruptura de Hoek-Brown, da rocha 
intacta representativas das 5 diferentes classes geomecânicas das encaixantes da camada de carvão 
Barro Branco na Bacia Carbonífera Sul-Catarinense. 


Embora na região de tração as classes geomecânicas Arenito LS e Siltito LS sejam 
semelhantes mecanicamente, na condição confinada a classe geomecânica Arenito LS 
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se assemelha ao Siltito HS, conforme observado na Figura 11. A explicação para esse 
comportamento na região de compressão da envoltória de ruptura pode estar associada 
a textura do Arenito LS (grãos maiores, mais angulosos e rugosos), quando comparado 
com as classes de siltitos. Esses aspectos texturais do Arenito LS tendem a produzir um 
aumento na dilatância durante a ruptura, com o consequente aumento da resistência ao 
cisalhamento. 

De acordo com a proposição de ZORZI et al (1991), conforme Figura 12, as 
classes geomecânicas Arenitos HS e MS e Siltito HS podem ser enquadrados no modelo 
geológico de rochas competentes. Já as classes geomecânicas Arenito LS e Siltito LS são 
enquadradas no modelo de rochas fracas. Essas duas classes representam cerca de 20% 
das rochas encaixantes presentes na bacia carbonífera e requerem análise da capacidade 
de carga no processo de dimensionamento dos pilares quando diretamente em contato 
com a camada de carvão. Também são as classes que exigem o sistema de suporte do 


teto mais robusto. 


Modelo a): não requer 


análise da capacidade 


ROCHA COMPETENTE carga. 
0. = 50 MPa 


o, = 5 MPa 
Em = 20 GPa 


CARVÃO 


ROCHA COMPETENTE 


ROCHA FRACA 
caca Modelos bje c): 


necessidade de análise da 
capacidade de carga 
ROCHA FRACA depende da posição da 


encaixante no teto, 
ROCHA COMPETENTE espessura de cada litologia 


CARVÃO e interferência do sistema 
de suporte adotado. 


Modelo d): requer análise 
da capacidade carga 


0. = 25 MPa 
o: = 2,5 MPa 
Em =6 GPa 





Figura 12 - a) Modelos básicos propostos para a análise da capacidade de carga das rochas 
encaixantes no dimensionamento de pilares das minas de carvão do Sul do Brasil. A figura representa 
modelos para o teto imediato (cerca de 5m). O mesmo modelo é válido para as condições do piso. b) 

Características geomecânicas médias de cada camada (ZORZI et al., 1991). 
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61 COMPORTAMENTO DAS ROCHAS DO TETO NAS MINAS DE CARVÃO 


As litologias classificadas geomecanicamente como arenitos HS and MS e siltito HS 
constituem os maciços rochosos de melhor qualidade quanto à estabilidade do teto das 
galerias. Quando no teto imediato, essas rochas são quase autosustentáveis (arenito) ou 
exigem sistema de suporte leve (siltitos). Quando acima do teto imediato, constituem os 
pontos de ancoragem do sistema de tirantes para a suspensão das camadas inferiores. 

Os maiores problemas de estabilidade do teto nas minas de carvão na Bacia Sul- 
Catarinense estão associados à presença do arenito LS ou do sitito LS. Esses materiais, de 
qualidade mecânica inferior, normalmente exigem intervenções com sistemas de suporte 
de elevada capacidade de carga, ou seja, parafusos mais longos e malhas mais apertadas 
(maior densidade de parafusos por área). Nesse tipo de maciço, os sistemas de tirantes não 
são eficientes quando o mecanismo de suporte é baseado exclusivamente no efeito viga. 
Nessas condições de teto, necessariamente é preciso recorrer a sistemas de tirantes mais 
compridos e que atuem por efeito suspensão. Esses tirantes são ancorados nas camadas 
mais resistentes acima do teto imediato, normalmente de arenitos com características 
mecânicas das classes arenitos HS e MS ou siltito HS. 

Os sistemas de suporte e configurações são diversos nas diferentes minas em 
operação. Cada mina adota o sistema que melhor se adapta a distribuição das diferentes 
litologias no teto, espessuras das camadas de rocha das diferentes classes geomecânicas, 
dimensões da escavação, capaz de proporcionar a estabilidade requerida para as diferentes 


àreas da mina, seja nos eixos principais, seja nos painéis de lavra. 


7 | COMPORTAMENTO DAS ROCHAS DO PISO NAS MINAS DE CARVÃO 


Praticamente todas as litologias estudadas podem estar presentes no piso da 
camada de carvão Barro Branco, embora predominem as das classes geomecânicas 
Arenito HS, Arenito MS e Siltito HS. Analogamente ao que ocorre no teto, o Arenito LS 
e, especialmente, o Siltito LS são as litologias sobre as quais se concentram a maior 
quantidade de problemas de alteração e colapso do piso. 

Um dos problemas está associado à baixa capacidade de carga para suportar as 
cargas imposta pelos pilares. Em casos extremos pode ocorrer a cravação do pilar (pillar 
punching) e colapso do piso da galeria. Outro aspecto importante está associado à alteração 
do siltito, na presença de água, e desagregação, devido ao intenso tráfego de máquinas 
pesadas. Esse tipo de problema ocorre predominantemente em áreas onde a camada de 
carvão apresenta-se inclinada. Nessas condições, os esforços de tração transferidos pelos 
pneus das máquinas são mais elevados, levando a ruptura da rocha do piso (peeling). 
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81 CONCLUSÕES 


Os estudos realizados permitiram a caracterização geomecânica das principais 
litologias presentes nas encaixantes da camada de carvão Barro Branco, da Bacia 
Carbonífera Sul-Catarinense. Pelo menos 5 classes geoemecânicas foram definidas em 
função das diferentes propriedades físicas e mecânicas determinadas em laboratório. 

As classes geomecânicas Arenito HS, Arenito MS e sSiltito HS representam 
aproximadamente 80% das rochas presentes nas encaixantes do teto e piso da camada de 
carvão Barro Branco. Essas classes também se constituem nas rochas de melhor qualidade 
mecânica. De outra parte, as classes Arenito LS e Siltito LS representam cerca de 20% das 
rochas encaixantes. De acordo com procedimentos sugeridos por ZORZI et al. (1991) para 
o dimensionamento dos pilares, quando em contato com a camada de carvão, as rochas 


encaixantes devem ser analisadas quanto à capacidade de carga. 
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ANEXO - SIGNIFICADO DA REPRESENTAÇÃO ESTATÍSTICA DE DADOS NA 
FORMA DE DIAGRAMAS DE CAIXA (BOX AND WHISKER) 

Essa técnica de representação de dados é extremamente útil para a observação da 
variabilidade ou a dispersão dos resultados. Esse tipo de representação é extremamente 
Útil para a identificação criteriosa de dados fora de contexto em uma amostra. A técnica 
utiliza 5 descritores básicos: mínimo, 1º quartil ou quartil inferior, mediana, 3º quartil ou 
quartil superior, e máximo. A média, opcionalmente, também pode ser representada. 

A Figura A1 explicita como são tratados os dados de uma amostra. A caixa retangular 
engloba a metade central dos valores (50% dos dados). No interior da caixa, a linha horizontal 
representa a mediana da amostra enquanto o sinal de (+) representa a média considerando 
todos os dados. A linha vertical, que parte da extremidade inferior da caixa, se estende do 
primeiro quartil até o menor (MIN) valor situado dentro de 1,5 vezes o intervalo interquartil 
(IQR=Q3-Q1). Analogamente, a linha vertical, que parte da extremidade superior da caixa, 
se estende do terceiro quartil até o maior (MAX) valor situado dentro de 1,5*(IQR). Valores 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 3 


que se situam entre 1,5*(IQR) (Lower fence ou Upper fence) e 3*(IQR) (Far Lower fence ou 
Far Upper fence), a partir dos respectivos quartis (Q1) e (Q3), são plotados como pontos 
individuais tratados como expúrios ou fora de contexto (OUTLIER), simbolizados por um 
quadrado ([1). Eventuais pontos que se situem além de 3*(IQR), abaixo do primeiro quartil 
(Q1) ou acima do terceiro quartil (Q3), constituem-se em valores extremos (FAROUTLIER). 

Ainda, para melhor entender o critério utilizado para a definição dos valores ouftliers 
ou extremos, a Figura A2 mostra a relação do diagrama de caixa e a função densidade de 
probabilidade para uma amostra com distribuição normal. Nota-se que os dados designados 
como outliers ou extremos, quando a distribuição é normal, correspondem a apenas 0,/%. 
Os demais 99,7% dos dados estão representados entre os quartis e 1,5*(IQR). 
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<< MIN 
E | IR) 
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ai io comi óuo n oto oê a 


151QR) below 25th percentile 


Figura A1 - Ilustração da representação estatística de dados através de diagramas de caixa (Box and 
Whisker) (adaptado de GALARNYK, 2018). 
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Figura AZ - Ilustração da relação entre o diagrama de caixa e a função densidade de probabilidade de 
uma distribuição de dados normal (adaptado de GALARNYK, 2018). 
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RESUMO: Sabe-sequeamineração de agregados 
para a construção civil, mais especificamente a 
brita, é bastante presente e ativa praticamente 
em todos os estados brasileiros. Porém é muito 
pouco o conhecimento desenvolvido acerca do 
potencial de utilização de seus sub-produtos. 
Estudos recentes apontam que esse material 
pode ser de grande utilidade na agricultura, 
devido às suas propriedades físico-químicas, 
reduzindo assim o consumo de fertilizantes 
industriais e fomentando o desenvolvimento da 
mineração sustentável. Este trabalho avaliou, à 
partir de técnicas de caracterização mineralógica 
e tecnológica, o potencial para utilização do pó 
de brita na agricultura como remineralizador 
de solo. Os resultados indicaram que este 
material apresenta grande potencial para este 
uso, segundo as exigências do Ministério da 
Agricultura, Pecuária e Abastecimento. 
PALAVRAS-CHAVE: Pó de rocha; mineração 
sustentável; agrominerais; pedreiras. 
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USE OF QUARRY ROCK DUST AS SOIL 
REMINERALIZER 


ABSTRACT: It is known that mining aggregates 
for civil construction, more specifically gravel, 
is very present and active in practically all 
Brazilian states. However, there is very little 
knowledge about the potential use of its by- 
products. Recent studies point out that this 
material can be of great use in agriculture, 
due to its physical-chemical properties, thus 
reducing the consumption of industrial fertilizers 
and promoting the development of sustainable 
mining. This work evaluated, using mineralogical 
and technological characterization techniques, 
the potential for using rock powder in agriculture 
as a soil remineralizer. The results indicated 
that this material has great potential for this use, 
according to the requirements of the Ministry of 
Agriculture. 

KEYWORDS: Rock Powder; Sustainable Mining; 
Agrominerals, Quarries. 


11 INTRODUÇÃO 


Sabe-se que no Brasil é gerada uma 
grande quantidade de resíduos provenientes 
do processo de beneficiamento das pedreiras. 
Esses produtos são comumente armazenados 
nas propriedades das mineradoras e/ou 
vendidos para um mercado específico (LUZ, 
2008). Entretanto, segundo Machado (2014), 
esses resíduos apresentam potencial uso na 
agricultura, setor de grande importância para a 


economia brasileira, como melhoradores e/ou 
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condicionadores de solos. 

Isso ocorre porque os solos brasileiros são geralmente ácidos, empobrecidos pela 
ação do intemperismo químico e muitas vezes carentes nos teores de fósforo e potássio. 
Esse fato gera uma necessidade de investimento alto em fertilizantes industrializados 
e corretivos de solo, que são, em sua maioria, importados. Sendo assim, a busca por 
alternativas que possam substituir total ou integralmente a dependência desses insumos 
é de grande importância. Com isso, a rochagem apresenta-se como uma prática agrícola 
viável de incorporação do pó de rocha ao solo empobrecido ou degradado, para promover 
sua remineralização (BENEDUZZI, 2011). 

Considerando que o pó de rocha tem como propósito principal fornecer nutrientes 
ao solo (IBRAM, 2016), o mesmo é também aceito como um fertilizante inteligente, visto 
que uma de suas características é o fornecimento gradual de nutrientes, de acordo com a 
necessidade da planta (BENEDUZZI, 2011). 

Portanto, é possível afirmar que a utilização de resíduos de pedreiras como 
remineralizadores de solo pode promover o aproveitamento de um material antes 
considerado um passivo ambiental, e até mesmo impulsionar o surgimento de um novo 
mercado consumidor regional para este material. Nos últimos anos, este tema tem sido 
muito estudado e debatido por instituições de pesquisa e Órgãos governamentais, de 
forma a subsidiar o conhecimento de diferentes tipos de pós de rocha, sua aplicação e 
também a sua regulamentação. O Ministério da Agricultura, Pecuária e Abastacimento 
(MAPA), através da IN nº5/2016, estabeleceu as especificações e garantias mínimas dos 
remineralizadores de solo, sejam eles sub-produtos de pedreiras ou rejeitos dos processos 
de beneficiamento. 

Dessa forma, este trabalho tem por objetivo avaliar o potencial do pó de migmatito 
de uma pedreira de brita na região nordeste do estado de SP como remineralizador de 
solo à partir da caracterização da rocha quanto à sua composição química, mineralógica e 
distribuição granulométrica. 


21 METODOLOGIA 


Foram coletados cerca de 30kg de produto “pó de pedra” (produto da segunda etapa 
de britagem, equivalente à brita zero), diretamente do britador secundário, em uma pedreira 
de brita para construção civil. Em campo foram também coletadas amostras de rocha para 
confecção de lâminas petrográficas. Em laboratório, esse material foi homogeneizado, 
quarteado e separado em três alíquotas de 300g cada para os ensaios subsequentes. 
Para caracterizar o material, foram realizadas análises mineralógica, petrográfica, química 
e granulométrica, conforme descrito a seguir. 

Os procedimentos metodológicos aqui discorridos basearam-se na análise potencial 
do material segundo os critérios estabelecidos na IN nº05/2016 do MAPA, onde ficam 
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estabelecidas as principais especificações e garantias mínimas dos remineralizadores, a 
saber: 
| — em relação à soma de bases (CaO, MgO, K,0): deve ser igual ou superior a 9% 
(nove por cento) em peso/peso; 


|| — em relação ao teor de óxido de potássio (K,O0): deve ser igual ou superior a 1% 
(um por cento) em peso/peso; 


|| — em relação ao SiO, livre presente no produto: não poderão conter teor superior 
a 25% (vinte e cinco por cento) em volume/volume; 


IV — em relação aos elementos potencialmente tóxicos: não poderão conter teores 
superiores a 15 ppm de Arsênio (As), 10 ppm de Cádmio (Cd), 0,1 ppm de Mercúrio 
(Hg) e 200 ppm de Chumbo (Pb). 


Além disso, quanto à natureza física dos remineralizadores, a mesma IN refere que 
o pó de rocha deve ter como garantia mínima: 


V - 100% das partículas passantes em peso/peso na peneira 2,0 mm (ABNTnº 10); 


VI - no mínimo 70% das partículas passantes em peso/peso na peneira 0,84 mm 
(ABNTnº 20); 


VII — no mínimo 50% das partículas passantes em peso/peso na peneira 0,3 mm 
(ABNTnº 50). 


Análise química 
Uma amostra de 300 g do pó de rocha foi destinada para o laboratório SGS-Geosol 


para a realização da técnica de espectrometria óptica de emissão atômica por plasma de 
argônio (ICP-OES), de elementos maiores e traços. 


Análise mineralógica 


Uma alíquota de 300 g foi destinada ao laboratório da Mineração Curimbaba 
Ltda. para realização da difração de raio X (DRX), que permite a identificação das fases 
minerais presentes na amostra através da caracterização de sua estrutura cristalina para 
quantidades superiores a 5%. 


Análise petrográfica 
Foi realizada à partir da análise de uma lâmina petrográfica da rocha sã em 


microscópio petrográfico. Para esta análise foi utilizado o microscópio Zeiss modelo Scope 
A1, do laboratório Multiusuário Ill da UNIFAL-MG, através de uma lente objetiva de 2,5x. 


Análise granulométrica 


Uma amostra de 300 g foi peneirada com o auxílio de um peneirador suspenso 
durante 10 minutos, seguindo as normas da ABNT. Para o procedimento utilizaram-se as 
seguintes aberturas de peneiras: 6 mesh, 8 mesh, 10 mesh, 14 mesh, 20 mesh, 28 mesh, 
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40 mesh e 60 mesh. 


31 DISCUSSÃO DE RESULTADOS 


Os resultados ensaio de peneiramento dispostos na Tabela 1 e Figura 1 mostraram 
que o material apresenta aproximadamente 65,09 % de partículas passantes na peneira 
de 2 mm. 

Essas informações elucidam que o material não cumpre as garantias mínimas para 
as especificações físicas exigidas no Anexo | Instrução Normativa nº 5/2016 do MAPA, 
precisando portanto de ajustes físicos para o uso em questão. 


AMO STRA (300, T 


Retido Passante Retido Passante 
Peneiras eg e ah) Peneiras | | 3) BC a ] acu E ) acumulado(*%) | acumu em Me] 


E RT 
pi [ni e o io 
Ba “RR SST DO SR RS E E 7 GS 


Lo Mm CT não 1] ai] ias | Ja JI ag II com i|] 
EEE E EEE 
EE A O, E TR 2 = 
32,44 es [24 [RR 
E BS TO RT A RS E PE 7 





Tabela 1: Dados da quantidade de material passante e retido no ensaio de peneiramento. 
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Figura 1: Curva granulométrica da amostra de pó de rocha. 


A análise química total realizada por ICP-OES na amostra (Tabela 2) registrou os 
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seguintes elementos maiores: CaO, Cr20,, Fe,O,, K,O, MgO, Na,0, P,O,, SiO,, TIO,, SIO,, 
ALO,, MnO, U,O,, ZrO,. Já como elementos traço foram registrados os elementos Ba, Nb, 
Sr, Zn, Zr, V, Ce, Co, Cs, Cu, Dy, Er, Eu, Ga, Gd, Gf, Ho, La, Lu, Mo, Nb, Nd, Ni, Pr, Rb, 
Sm, Sn, Ta, Tb, Th, TI, Tm, U, W, Y e Yb, segundo a tabela 2a — d. Segundo a Instrução 
Normativa nº 5/2016 do MAPA, a soma das bases (CaO, MgO, K,O) deve ser superior a 
9%. Os resultados aqui apresentados indicam que no caso em questão a soma das bases 
é de 8,6%, ou seja, num valor muito próximo limite mínimo solicitado. 

No cumprimento integral da IN à partir da análise química, tem-se que o material 
contempla o teor mínimo exigido de pelo menos 1% para K,O, e que também não registrou 
teores dos elementos potencialmente tóxicos de Ar, Cd, Hg e Pb à partir da técnica utilizada. 

Os resultados das análises mineralógica (realizada através da técnica de Difração 
de raio X) e petrográfica indicaram a presença de quarizo, plagioclásio, feldspato potássico, 
biotita, flogopita e hornblenda como minerais principais e, como minerais acessórios, 
zircão e opacos. A análise petrográfica foi realizada através da contagem volumétrica, 
cujos resultados da proporção volumétrica dos minerais identificados, estão dispostos na 
Tabela 3. 


Blank SJ-01 SJ-01 (REP) 
AIZO3 (%) 0,03 13,65 14.26 
Ba (ppm) <10 800 875 





Na20 (%) <0,01 3 3.15 
Nb m NA. NA. N.A. 
P205 (%) <0,01 0,13 0,15 
SiO2 (%) >90 69,93 70,76 
Sr (ppm) <10 215 235 
TIO2 (%) <0,01 0,48 0,52 
Zn (ppm) 10 62 73 
Zr (ppm) <10 260 261 
V (ppm) 6 51 56 
LOI (%) 0,24 0,35 n.d 
Ce (ppm) 04 101,4 n.d 
Co (ppm) 0,7 9,3 n.d 
Cs (ppm) 0,2 0,81 n.d 
Cu (ppm 8 21 n.d 
Dy (ppm) <0,05 4,25 n.d 
Er (ppm) <0,05 2 n.d. 
Eu (ppm) <0,05 1,32 n.d 
Ga (ppm) 0.2 16.4 n.d 
Gd (ppm) 0,1 5.72 n.d 
Hf (ppm) <0,05 6,49 n.d 
Ho (ppm) <0,05 0,74 n.d 
La (ppm) 2,1 51,7 n.d 
Lu (ppm) <0,05 0,26 n.d 
Mo (ppm) 2 3 n.d 
Nb (ppm) 1,33 11,06 n.d. 
Nd m 0,2 39,2 n.d 
Ni (ppm) 16 41 n.d 
Pr (ppm) 0,06 11,76 n.d 
Rb (ppm) 0,5 146 n.d 
Sm (ppm) <0,1 7,3 n.d 
Sn (ppm) 04 2,3 nd 
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Tabela 2: Resultados da análise química por ICP-OES. 


Mineral Fórmula quimica Porcentagem 
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Tabela 3: Proporção volumétrica da mineralogia básica e acessória identificada na análise petrográfica. 


41 CONCLUSÃO 


Os resultados alcançados através das análises químicas estão dentro dos padrões 
aceitados pela Instrução Normativa do MAPA, conforme descrito nos incisos Ile Ill do Art 4º. 
Dessa forma, entende-se que a química e a mineralogia do material de estudo apresentam 
um quadro favorável para a prática da rochagem. Ainda com relação aos resultados das 
análises mineralógica e petrográfica, observou-se que foi cumprido o critério estipulado no 
Art. 4º, 8 2º, inciso | da Instrução Normativa nº 5 do MAPA. 

Entretanto, observa-se que a granulometria das partículas não está de acordo com 
o que foi proposto pela norma. Mais especificamente, o material possui uma distribuição 
granulométrica imprópria para utilizá-lo como remineralizador. Para melhorar essa 
distribuição poderiam ser realizadas cominuições utilizando outros equipamentos, tais 
como britador de rolos de alta pressão e ou até mesmo peneiras de alta frequência, para 
melhor classificação do material. 

A aplicação deste material como remineralizador apresenta grande potencial, mas 
entretanto precisa de ajustes físicos em seu processo de beneficiamento para cumprimento 
da norma. Se utilizado também para esse fim, haverá grande contribuição para o não 
acúmulo de passivos ambientais, além de um novo mercado que auxilie na sustentabilidade 


das atividades da empresa. 
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RESUMO: POME (Palm Oil Mill Effluent) é um 
resíduo do processo de aceitação de extração 
de palma; Para 2018, foram geradas 4.887.000 
toneladas destes resíduos a nível nacional, O 
seu elevado teor de matéria orgânica (DQO 
médio de 60.000 mg 02 /L) torna-o um problema 
ambiental que deve ser resolvido com urgência. 
Para o POME de Puerto Salgar (Colômbia), foi 
encontrado um COD de 69650 mg O2/L e 31543 
mg SV/L, o que o torna uma matéria-prima com 
alto potencial metanogênico. No processo de 
biodigestão, foi alcançada uma produção média 
de 25,3 LBiogás / LPOME com uma composição 
de 58,4% CH4, 40,4% CO2, 1,1% 02 e 1156 ppm 
de H2S, apresentando uma redução na carga 
orgânica de 93%. Os resultados obtidos, além da 
alta disponibilidade desse resíduo, tornam-no um 
substrato adequado para processos de digestão 
anaeróbia voltados para a geração de energia a 
partir de fontes não convencionais. 
PALAVRAS-CHAVE: POME, 
anaerobia, biogás, DQO, UASEB. 


digestión 


BIODIGESTION OF POME AS AN 
ENERGY AND ENVIRONMENTAL 
ALTERNATIVE IN PALM OIL 
EXTRACTORS 


ABSTRACT: POME (Palm Oil Mill Effluent) is a 
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residue from the palm oil extraction process; by 2018 4'887.000 tons of this residue were 
generated nationally, its high organic content (COD average 60000 mg O2/L) makes it an 
environmental problem that must be addressed urgently. For POME from Puerto Salgar 
(Colombia), a COD of 69650 mg 02/L and 31543 MO, jaos volatie!L Was found, which makes 
it a raw material with high methanogenic potential. In the biodigestion process, an average 
production of 25.3 Lagas ! LPOME was achieved with a composition of 58.4% CH4, 40.4% CO2, 
1.1% O2 and 1156 ppm H2S, with an organic load reduction of 93%. The results obtained, 
together with the high availability of this residue, make this a suitable substrate for anaerobic 
digestion processes oriented to the generation of energy from unconventional sources. 
KEYWORDS: POME, anaerobic digestion, biogas, COD, UASEB. 


INTRODUCCIÓN 


En el ao 2018 en Colombia se registró una producción de 1.7 millones de toneladas 
de aceite de palma lo que generó un residuo de 4.9 millones de toneladas de POME; es 
importante resaltar el incremento de la producción de aceite de palma en los últimos aos, 
lo cual ha Ilevado a un aumento de 46.6% de generación de POME desde el aho 2014 [2]. 

El aumento de la producción ha convertido el POME en un problema ambiental 
importante, debido a que presenta altos niveles de contaminación, pues es depositado 
en lagunas de oxidación las cuales requieren de grandes extensiones de tierra y largos 
tiempos de degradación en contacto directo a la atmósfera arrojando dióxido de carbono y 
metano, los cuales representa una gran parte de los gases de efecto invernadero [3]. En el 
presente trabajo se expone la caracterización y el aprovechamiento del potencial energético 
de POME (Palm Oil Mill Effluent), el cual es un residuo proveniente principalmente de la 
etapa de esterilización y clarificado en el proceso de extracción de aceite de palma, del 
cual se estima que por cada tonelada de aceite crudo producido se obtienen en promedio 
3 toneladas de POME [1]. 

Sus características fisicoquímicas y su alta carga orgânica, representada 
principalmente en residuos lipídicos, hacen del POME un sustrato de gran interés para 
ser sometido a un proceso de digestión anaerobia como una forma de aprovechamiento 
energético y así revertir el impacto negativo de estos residuos en las comunidades 
aledahas a las plantas de procesamiento del fruto y los dafos causados al ecosistema 
acuático disminuyendo la carga orgânica de estos efluentes y permitiendo encontrar una 
aplicación agroindustrial de las corrientes obtenidas, generando así un valor agregado para 
el proceso. 

Muchos autores presentan la digestión anaerobia como una alternativa promisoria 
para el tratamiento de los residuos de las plantas extractoras de aceite de palma. Pogaku, 
et al., (2015) alcanzaron una producción de biogás de 0.7825 litros por kilogramo de POME 
con una composición de 62.5% CH,, 37.5% CO,, y trazas de H,S [4]. Poh and Chong (2014) 
en un reactor tipo UASB-HCPB obtienen un biogás con un 52% de metano y una remoción 


de 88% de DQO [7]. En el 2011 Fang, et al. alcanzó rendimiento de metano a partir de 
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POME de 438 CH4/g SV con una remoción de 90% en DQO [8]. 

En el presente trabajo se realiza la caracterización y aprovechamiento energético 
del POME, por medio de un proceso de digestión anaerobia para la producción de biogás 
que pueda suplir demandas energéticas en las plantas de extracción de aceite de palma. 


METODOLOGIA 


Material 


El POME es proveniente de una extractora ubicada en el departamento de 
Cundinamarca (Colombia). 
Caracterización 


Se caracterizó el POME siguiendo los lineamientos del Standard Methods [5], 
realizando pruebas por triplicado de los factores más relevantes como solidos totales (ST), 
solidos volátiles (SV) y solidos fijos (SF), DQO, carbono orgánico, nitrógeno orgánico, pH, 
DBO. 


Evaluación del POME como sustrato 


El proceso se realiza en un reacior de flujo ascendente continuo (UASB) de segunda 
generación acoplado a un reactor secundario; el reactor secundario cumple la función de 
finalizar la degradación de materia orgânica que no se degradó en el UASB, el cual trabaja 
con un flujo diario de aproximadamente 3.96 L distribuido a lo largo de las 24 horas, este 
flujo contiene en promedio 2.03 L de POME. En la Figura 1, se presenta un esquema del 
sistema que se utiliza para la degradación anaerobia del POME. 


1. Reactor UASB 

2. Corrente para alimentación del 
UASB 

3. Corrientes de biogás producido en 
UASB 

4. Brazo del UASB, hecho para evitar 
el paso de sólidos alreactor secundario. 
5. Corrente de alimentación de 
efluente del UASB al reactor 
secundario. 

6. Reactor secundario (LL) 

?. Corrente de salida de biogás del 
reactor secundario 

8. Corriente de salida del efluente del 
proceso. 

9. corrente para toma de muestra del 
efluente UASE. 

10. Corrientes para perfil de lodos en el 
UASB. 

11. Corrientes para perfil de lodos en el 
UASB. 

12. Corrientes para perfildelodos en el 
UASB. 

13. Bolsa colectora de biogás. 





Figura 1. Diagrama del sistema de digestión anaerobia UASB. 
Se toman muestras semanalmente, para realizar un control del sistema por medio 
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de la cuantificación de los parâmetros que se consideran pueden afectan la producción de 
biogás de manera directa, esta cuantificación se realiza principalmente en las corrientes 
(2), (8), (9) y (13) que se citan en la Figura 1. La cuantificación de estas propiedades se 
rige bajo los estândares establecidos por el Standard Methods [5], siendo el método (2540) 
utilizado para la determinación de los sólidos (totales, volátiles y fijos), el método 8000 
para determinación de DOO y mediciones de otras variables pH y relación ácidos grasos 
volátiles y alcalinidad por el método de Kapp [6]. 

El análisis se lleva a cabo tomando 600 mL de muestra en la corriente 2, 8 y 9 que 
corresponde a la alimentación, lecho de lodos y UASB, respectivamente. De cada muestra 
se disponen 150 mL para el análisis de sólidos totales y volátiles, cada análisis se realiza 
por triplicado; cada muestra se seca en estufa a 105ºC por 24 horas para cuantificar sólidos 
totales y posteriormente en una mufla a 550ºC durante 1 hora para cuantificar sólidos 
volátiles. 

Para DQO se toman 20 mL de muestra y se realiza el análisis por triplicado, después 
de preparar las muestras se agregan a los viales de DQO con micropipeta un volumen de 
20 yuL y se ponen a reaccionar en un termoreactor por 2 horas a 150ºC, luego de cumplirse 
el tiempo se dejan reposar hasta llegar a temperatura ambiente y se mide la concentración 
en un colorímetro de la misma marca del termoreactor' el cual tiene como blanco un vial 
de DQO sin muestra. 

Para la medición de pH, cuantificación de ácidos grasos volátiles y alcalinidad, 
se necesitan alrededor de 125 mL de cada muestra, se centrifugan por 10 minutos y el 
sobrenadante se titula con ácido sulfúrico 0.1N, este proceso realiza por triplicado. Para el 
cálculo de ácidos grasos volátiles y alcalinidad se hace uso del método Kapp [6]. 

Con los datos anteriores se lleva el control de la producción de biogás la cual puede 
variar debido a inhibiciones que pueda presentar el sistema como aumento en la relación 
ácidos grasos - alcalinidad, al tener el reactor en funcionamiento se mide la producción 
específica de biogás cuantificando esta por medio de un medidor de gas? y cuantificación de 
la composición del mismo por medio de un equipo”, lo cual nos permite saber la proporción 
de metano, dióxido de carbono, oxígeno y sulfuro de hidrógeno que contiene el biogás. 


RESULTADOS 


La caracterización del POME se muestra en la Tabla 1. Los resultados del monitoreo 
del reactor se presentan en las figuras 2, 3, 4, 5 y 6, donde se muestra principalmente 
el comportamiento de las diferentes variables como solidos volátiles (SV), DQO, ácidos 
grasos volátiles (AGV) y alcalinidad (ALK) y la producción de biogás durante el tiempo de 
monitoreo. 


1 Marca HACH 
2 Marca Premex 
3 GIOTECH BIOGÁS 5000 
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; Desviación 
ST(mg/L) 39333 + 4481 


SV (mg/L) 31236 + 3249 


SF (mg/L) + 1360 
DQO (mgO;/L) 69650 


DBO (mgO,/L) 8500 


Nitrógeno orgânico 
(g/L) 
Relación C/N 





-——4— Alimentación 


20 q “ 
| | —eUASB 


SV (mg/L) x 10º 


| ' —&— Lecho de 


Figura 2. Comportamiento de Solidos volátiles en el reactor. 


En la Figura 2, se observa un pico en el día 55, que representa un incremento en 
los SV en el UASB, esto como resultado de una sobrecarga en el reactor; sin embargo, 
la cantidad de SV en el lecho de lodos (reactor secundario) tiene valores que van con el 
comportamiento normal del reactor, lo que permite decir que esta etapa garantiza obtener 
una depuración promedio en los sólidos volátiles (SV) en todo el proceso al controlar las 


sobrecargas. 
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Figura 3. Relación de solidos volátiles y solidos totales durante tiempo de monitoreo del proceso. 


Figura 4. Demanda química de oxígeno del proceso durante tiempo de monitoreo del proceso. 


De acuerdo a la Figura 4, se puede apreciar que la mayor depuración de materia 
orgânica se realiza en el UASB y que el reactor segundario es óptimo para controlar 
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sobrecargas sin necesidad de detener el proceso cuando se presentan. 
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Figura 5. Ácidos grasos volátiles (AGV) y alcalinidad (ALK) durante el monitoreo del proceso tanto en el 
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reactor (UASB) como en el lecho de lodos (LL). 
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Se mantiene un control de los niveles de AGV y Alcalinidad para evitar sufrir 


inhibiciones en el sistema por el aumento de esta relación, se puede apreciar que pese a la 


sobrecarga sufrida el UASB no se ve descompensado y tiende a estabilizarse. 
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Figura 6. Producción de biogás; A: biogás / volumen de POME, B: producción especifica (volumen 


biogás / (volumen de trabajo — día) 


La producción de biogás se ve afectada por el sistema que se empleo para su 


colecta, siendo modificado (día 40), y dando como resultado una mayor producción debido 


a que no se contaba con un deposito adecuado que permitiera la adecuada cuantificación 


de este. 


CONCLUSIONES 
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A pesar de ser un experimento exploratorio, la biodigestión del POME produce 
un biogás con una cantidad de metano apreciable (58.4% CH,, 40.4% CO,, 
1.2% O,). 


Se logra una remoción de materia orgânica promedio de 93%, indicando que el 
proceso de digestión anaerobia es un buen proceso para disminuir su impacto 
ambiental. 


El uso de un reactor secundario garantiza el control de sobrecargas que pue- 
da sufrir el sistema, en donde la cantidad de materia degradada en el UASB 
disminuye, pero al final del proceso se observa que la remoción se mantiene 
promedio a la que presenta el sistema cuando está estable. 


Los procesos de digestión anaerobia dan a residuos como el POME un valor 
energético agregado, contribuyendo así a la valorización de residuos orgánicos 
y revierten el impacto ambiental que muchos causan tras su disposición final 
al reducir su carga orgânica que se ve reflejada en la disminución de DQO y 
sólidos volátiles. 
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RESUMO: Este estudo visa incrementar o 
interesse de estudantes do ensino médio e reduzir 
a evasão do curso de engenharia civil através 
da realização de aula prática de montagem de 
unidade sanitária de plástico reciclado pelos 
discentes da disciplina de Tecnologia das 
Construções contando com a participação dos 
estudantes de uma escola de ensino médio. 
Uma parceria entre uma Universidade pública, o 
fabricante do banheiro de plástico e uma escola 
pública do semiárido do nordeste brasileiro foi 
estabelecida através de um convênio com o 
CNPq e a Companhia Vale do Rio Doce — Forma 
Engenharia - com o fito de realizar análises do 
processo de montagem e desmontagem do 
banheiro de plástico reciclado. A primeira aula 
prática foi realizada pelos técnicos da empresa 
doadora da unidade sanitária. As observações 
dos participantes durante a referida aula serviram 
de base para elaboração do procedimento 
para realização da aula prática de montagem 
e manual do banheiro de plástico reciclado. 
Duas aulas práticas foram realizadas com apoio 
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EXPERIÊNCIA 


dos técnicos da Universidade que realizou a 
pesquisa, e, duas foram realizadas somente 
pelos estudantes do curso de Engenharia civil 
da referida Universidade. Em todas elas ficou 
evidente a dificuldade de encaixe da viga superior 
em virtude da altura e da necessidade de ajuste 
entre os quatro pilares. A partir desse estudo a 
empresa fornecedora realizou testes no protótipo 
e modificou o material dos pilares de ferro para 
alumínio. 

PALAVRAS-CHAVE: Engenharia civil. Educação 
e Aula prática. 


CLASSROOM PRACTICE ABOUT HEALTH 
UNIT ASSEMBLY (BATHROOM) FROM 
RECYCLED PLASTIC — EXPERIENCE 

REPORT 


ABSTRACT: The bass interest of high school 
students in engineering courses is growing 
in Brazil. Moreover, a survey carried out at 
the Polytechnic School of the Universidad de 
São Paulo (USP) with students from the Civil 
Engineering course, about evasion in that course 
showed that students feel that there is excess 
basic subjects and less focused on practice they 
perceive as civil Engineering. In this light, it is 
necessary to take action to increase the interest 
of students in high school and beyond reduce 
evasion in engineering course. As well, this study 
reports the actions taken in the civil engineering 
course, through a practical class assembly of a 
sanitary unit made of recycled plastic by students 
of the construction technology discipline with 
the participation of students from a high school. 
lt has been the model studied a partnership 
between a public university, the plastic bathroom 
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manufacturer and a public school in the semiarid region of northeastern Brazil. The above 
was established in an agreement with the National Research Center - CNPq and Vale do Rio 
Doce company with the aim of carrying out analysis of the assembly process and disassembly 
of recycled plastic bathroom. This study used the technique of participant observation. The 
observations were the basis for preparing the assembly procedure and bathroom manual 
recycled plastic. It concludes by highlighting the difficulties at the time of bathroom assembly. 
From this study the company bathroom idealizing conducted tests on the prototype and 
modified the material of iron pillars for aluminum. 

KEYWORDS: Civil engeener, education, bathroom from recycled plastic. 


11 INTRODUÇÃO 


A formação do engenheiro civil é considerada fundamental para o sucesso e a 
remuneração do profissional no mercado, porém a taxa de evasão em engenharia civil é 
alta e muitos desistem antes do fim do curso. 

Uma pesquisa realizada na Escola Politécnica da USP com estudantes do curso 
de Engenharia civil sobre evasão no referido curso evidenciou que os alunos consideram 
que há excesso de matérias básicas e pouco voltadas à prática do que entendem como 
engenharia civil. Outro problema brasileiro se concentra no baixo interesse dos estudantes 
de ensino médio pelos cursos de engenharia. 

Diante do exposto este estudo visa incrementar o interesse de estudantes do ensino 
médio e reduzir a evasão do curso de engenharia civil através da realização de aula prática 
de montagem de banheiro sanitário de plástico reciclado pelos estudantes da disciplina de 
Tecnologia das Construções do Curso de Engenharia Civil contando com a participação 


dos estudantes de uma escola de ensino médio. 


21 NOVAS TECNOLOGIAS E ENGENHARIA CIVIL 


Sustentabilidade implica em sistemas construtivos que promovam integração com o 
meio ambiente, adaptando-os para as necessidades de uso, produção e consumo humano, 
sem esgotar os recursos naturais, preservando-os para as gerações futuras; além da 
adoção de soluções que propiciem edificações econômicas e o bem-estar social (CÂMARA 
DA INDÚSTRIA DA CONSTRUÇÃO, 2008). 

Há alguns anos reforçou-se uma preocupação em termos de Meio Ambiente, 
poluição, resíduos e desenvolvimento sustentável; fato que gerou, em grande parte da 
população e em vários setores, a necessidade de inovar, procurar soluções renováveis, 
que provocassem menores impactos ou que contribuíssem de certa forma para minimizar 
os impactos gerados no meio. 

O setor da Construção Civil é um dos que produzem maior quantidade de resíduos 
sólidos que ficam dispostos em lixões ou aterros sanitários sem que lhes seja dado 
reaproveitamento, e para os quais ocorre reaproveitamento, constitui-se métodos difíceis 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 6 EEN 


de serem executados e por vezes caro, tornando essas práticas inacessíveis à população. 

A fabricação de diversos produtos com a utilização de material reciclado vem 
ganhando cada vez mais importância e reconhecimento diante de uma realidade de um 
planeta cada vez mais sedento por sustentabilidade. Há algumas técnicas construtivas 
nesse sentido que têm se inserido e ganhado espaço no setor da construção civil, como 
a técnica dos módulos habitacionais, com os quais se pode construir: canteiro de obras, 
banheiro, residência, escritórios, os quais são alguns produtos utilizados em obras 
nacionais, todos fazendo uso do plástico reciclado como matéria prima. 

O método dos módulos habitacionais tem se aperfeiçoado, além de viável e inovador, 
é bastante simples de ser implantado, ou seja, qualquer pessoa seguindo o manual pode 
montar os produtos; além de ser uma medida sustentável, uma forma alternativa de 
reaproveitar materiais anteriormente sem destinação precisa, tornando-os materiais de 


grande utilidade, reduzindo os impactos ambientais. 


31 EDUCAÇÃO FORMAL, NÃO FORMAL E INFORMAL 


Aprender traz consigo a possibilidade de algo novo, incorporado ao conjunto 
de elementos que formam a vida do indivíduo, relacionando-se com a mudança de 
conhecimentos que ele já possui (OLIVEIRA, 2011). 

A ocorrência de educação formal, não formal e informal é independente dos espaços 
onde ocorrem. Assim as três modalidades de educação podem se dar em espaços formais 
e em não formais como equivalentes a espaços escolares e espaços não formais como 
qualquer outro espaço externo à escola (OLIVEIRA e GASTAL, 2009). 

A integração entre a prática e a teoria é necessária para integrar os conteúdos que 
permeiam o curso, além de propiciar a compreensão dos mesmos por meio de vivências do 
cotidiano da área de estudo. (AGUIAR, et.al., 2016). 

A educação não formal é uma atividade aberta, que vem sendo moldada por 
educadores das mais diversas áreas. Ela não é estática, e isso permite o desenvolvimento 
da modalidade, bem como da forma como se processa o ensino em momentos em que os 


ambientes formais não são suficientes para que o processo de aprendizagem ocorra. 


41 UNIDADE SANITÁRIA EM PLÁSTICO RECICLADO 


A unidade sanitária de plástico reciclado, em seu modelo convencional, conta 
com montantes de alumínio, geralmente com 2,70 m de altura; Juntas de encaixe que 
também são constituídas de material metálico formando perfil H revestido com pvc, que 
ficam entre uma placa de plástico e outra, e entre as cintas de alumínio e as placas de 
plástico, facilitando os encaixes, a mesma possui 84 cm de comprimento. Para a base é 


colocado uma cinta de metalon galvanizado, que já apresenta os encaixes necessários 
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para os montantes de alumínio; o piso é constituído por placas de plasterit encaixadas. E, a 
parte superior é constituída por uma cinta de amarração também de metalon galvanizado. 
Nas partes dos nós entre a ligação das cintas é colocado parafuso autobrocante utilizando 
furadeira. 

A caixa d'água fica localizada no interior da unidade sanitária suspendida na parte 
superior por meio de uma estrutura metálica apoiada na cinta superior, a caixa d'água é de 
200 litros, se tratando de um uso básico, convencional; sendo que sua capacidade máxima 
dependerá da necessidade do projeto, onde poderá ser feito reforço do metalon de suporte 
de acordo com o esforço demandado. A parte de instalações elétricas é feita normalmente, 
assim como a hidráulica contendo bacia sanitária, chuveiro, pia e tudo o que se utiliza 
comumente em banheiros convencionais; porta conforme a norma de edificações da norma 
de acessibilidade 9050, de 80 cm de largura e saída de ar. 

A obra finalizada da unidade sanitária apresenta 1,83 de comprimento por 1,53 de 
largura. Figura 1. A estrutura pode ser tanto comprada como alugada para obras e os 


módulos podem ser montados e desmontados quantas vezes se fizerem necessárias. 


FRODUCED Br AN AUTODESK EDUCATIONAL PRODUCT 


dale E 5 II 
ã 
E] 
a 
b 
E 
e ” 
- N ! 
» 
+ | 
» 
tejo 


TI H E 


Figura 1 — Módulo do banheiro de plástico reciclado 


Fonte: Impacto Protensão (2016). 


A instalação elétrica pode ser feita por cima e pode ficar embutida nos perfis de 
alumínio que já possuem furações predeterminadas para as instalações de tomada baixa 
ou alta. As instalações hidros-sanitárias ficam expostas no interior ou exterior do contêiner 
plástico. A vasão de descarte pode ser feita ligada à rede de esgoto ou através de fossas 
sépticas. As instalações hidráulicas são realizadas junto com a colocação das bases. 

Apesar do banheiro de plástico ser totalmente modularizado, há possibilidade de se 


fazer algumas adaptações para instalações elétricas e hidráulicas. No piso com modulação 
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plástica pode haver necessidade de corte do plasterit. Para passagem das instalações 
hidráulicas através do piso ou das paredes utiliza-se uma serra-copo para perfurações. 
Caso as especificações do projeto sejam passadas com antecedência os ajustes de cortes 
são realizados na própria empresa que fabrica o material, de forma a aumentar a rapidez 
da montagem e evitar a necessidade de mão de obra especializada. 


5 | ESTUDO DE CASO 


A disciplina de Tecnologia das Edificações do oitavo semestre do curso de Engenharia 
Civil visa proporcionar aos discentes conhecimentos das etapas de execução dos serviços 
de uma obra de construção civil com visão dos processos, linguagem técnica e ferramentas 
básicas para gestão e produção de edificações. 

O conteúdo programático é dividido em duas partes denominadas de obra bruta e 
fina. 

As principais estratégias de ensino utilizadas são: Aulas expositivas com discussão 
de conceitos e estudos de caso; Leitura e interpretação de textos; Seminários dos alunos 
para apresentação de artigo científico e de estudos de caso e/ou trabalhos práticos 
realizados por eles e Visitas técnicas. 

Na aula prática de montagem e desmontagem do banheiro de plástico reciclado, que 
é o foco deste estudo, procura-se observar na prática os processos construtivos estudados 
em sala de aula e também inovações tecnológicas utilizadas nas obras. 

A primeira aula prática foi realizada pelos técnicos da empresa doadora do referido 
banheiro. As observações dos participantes durante a referida aula serviram de base 
para elaboração do procedimento para realização da aula prática de montagem e manual 
do banheiro de plástico reciclado. Duas aulas práticas foram realizadas com apoio dos 
técnicos da Universidade que realizou a pesquisa, e, duas foram realizadas somente pelos 
estudantes do curso de Engenharia civil da referida Universidade. Figura 2. 





Figura 2 — Montagem do módulo do banheiro de plástico reciclado. 
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Figura 3 — Colocação da viga metálica superior — andaimes sem tábuas. 


A figura 3 mostra a finalização da montagem do módulo do banheiro de plástico 
reciclado evidenciando a colocação da viga metálica e a ausência de tábua para maior 
segurança e apoio no andaime. 

É necessário difundir a abordagem dessas novas tecnologias para construtoras e 
profissionais da área de engenharia, porém também é bastante interessante que essas 
práticas comecem a ser abordadas e incentivadas desde a graduação, para os futuros 
profissionais, de forma que incentive a utilização futura. 

Tendo em vista a importância do incentivo da aplicação de novas tecnologias, tem 
sido realizada uma prática da apresentação da unidade sanitária de plástico reciclado em 
uma Universidade Federal da cidade de Mossoró, em parceria com uma instituição pública 
de ensino médio 

O banheiro de plástico foi selecionado dentre os diversos produtos que são 
produzidos através da reciclagem do plástico por ser de porte menor que o de um módulo 
habitacional e por envolver parte de instalações, uma vez que constitui uma abordagem 
importante para os alunos. 

O projeto visa promover a inclusão de alunos nos cursos de engenharia apresentando 
o curso como inovador e dinâmico. Busca despertar maior interesse dos alunos que já 
estão na Universidade, de forma a minimizar as frequentes evasões. 

O planejamento da aula prática inicia-se pela organização das atividades a serem 
desenvolvidas, através da elaboração do cronograma de atividades (Tabela 1). 
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Cronograma de Atividades 
1 Planejamento das atividades e parcerias 


Tabela 1 — Cronograma das Atividades 





Fonte: Relatório de pesquisa— Edital 005/ 20012 CNPq /Vale — Forma Engenharia 


Para realização da aula prática foi necessário estudar o processo de montagem do 
banheiro de plástico reciclado e elaborar o procedimento para realização da aula prática de 
montagem e o manual do banheiro de plástico reciclado. Também foi realizada a divulgação 
da referida aula para os alunos do ensino médio, da escola correspondente. As aulas 
práticas ocorrem semestralmente, com a presença dos estudantes de Engenharia Civil e 
os alunos de ensino médio da escola parceira. 

Com base nas orientações de montagem da empresa fornecedora da umidade 
sanitária de plástico reciclado foi elaborado Procedimento de Montagem do Banheiro de 
Plástico Reciclado, descrito a seguir: 


1. Posicionar a cinta metálica de amarração inferior e a tubulação sanitária; 
2. Concretar (Montar) a estrutura do piso; 


3. Posicionar os perfis verticais de PVC (Cloreto de Polivinila) reforçados internamente 
com tubo industrial metálico quadrado para estruturação das paredes; 


4. Encaixar as placas plásticas 0,61cm x 0,61cm em polipropileno EP200K reciclado 
preenchidas com poliestireno expandido (EPS) com isolante térmico acústico 
nos perfis metálicos do piso e nos perfis verticais de PVC (Cloreto de Polivinila) 
reforçados internamente com tubo industrial metálico quadrado para estruturação 
das paredes; 


5. Distribuir as esquadrias (portas e janelas) de acordo com o projeto de montagem; 
6. Posicionar a cinta metálica de amarração superior; 

7. Posicionar a caixa d'água; 

8. Fixar a telha metálica; 

9. Distribuir as instalações, hidros-sanitárias, elétricas e as louças e metais sanitários; 


10. Fixar através de parafusos as cintas metálicas inferiores e superiores nos perfis 
verticais de PVC reforçados internamente com tubo industrial metálico quadrado 
para estruturação das paredes; 
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11. Fixar a coberta com telhas trapezoidais em aço galvanizado; 


O processo de desmontagem foi realizado em agosto de 2013 pela equipe do projeto 
com auxílio de técnicos do laboratório de tecnologia das construções e de instaladores da 
Universidade correspondente. O tempo e as dificuldades nos procedimentos de montagem 
foram similares aos de montagem. 

Realizou-se o acompanhamento da desmontagem e montagem do banheiro de 
plástico reciclado. Esta tarefa foi executada por duas turmas da disciplina de Tecnologia das 
Edificações de uma Universidade Federal da cidade de Mossoró, com o acompanhamento 
de alunos do ensino médio de uma instituição publica de ensino médio. Esse processo 
constitui uma parceria que compõe um projeto, cujo objetivo é atrair alunos para o curso de 
engenharia civil na instituição e minimizar as evasões que ocorrem frequentemente. 

A primeira estratégia para obter bom desempenho na tarefa é observar bem a forma 
como se dá a montagem, seja através do manual ou do vídeo que explicam a montagem e 
realizar um planejamento rápido com base nas ferramentas que se pode utilizar. O segundo 
passo deve ser organizar todo o material ao redor da obra, para facilitar o manuseio sem 
desperdícios de tempo. Aplicando essas estratégias e juntamente obtendo empenho e 
dedicação da parte da equipe de obras, é possível que a realização da montagem ocorra 
de forma fácil, eficaz e com economia de tempo. 

O primeiro processo de montagem ocorreu com a turma de 20183, a qual realizou o 
referido processo com base em um manual de montagem do próprio produto reciclado. No 
procedimento observou-se que houve cooperação mutua e trabalho em equipe, porém não 
foi dada tanta importância a questão da organização do canteiro, o que acabou sendo um 
fator relevante, pois como inicialmente não pesou-se tanto essa questão, a turma acabou 
executando a tarefa com o dobro de tempo, se comparada à turma de 2016. 

A segunda turma foi acompanhada durante o processo de montagem em 2016. 
Dessa vez os alunos não tiveram acesso ao manual de montagem, mas assistiram um 
vídeo do processo de montagem da unidade sanitária de plástico reciclado. Esta turma 
por sua vez apresentou melhor desempenho na realização da tarefa, pois foram bastante 
reforçados os conceitos comumente utilizados ao longo da graduação, quanto o trabalho 
em equipe, a organização do material o mais próximo possível do local onde está sendo 
realizada a obra, OS prazos, o incentivo mútuo da equipe de obra. Dessa vez o tempo de 
desmontagem foi de 15 minutos. 

A montagem e desmontagem do banheiro são procedimentos simples que com 
poucas ferramentas é possível realizar. Das etapas de montagem a que a equipe apresentou 
maiores dificuldades foi com o encaixe da cinta metálica superior, pois para sua realização 
todas as peças deveriam estar bem alinhadas. 

Outra dificuldade com relação à desmontagem e montagem da unidade sanitária, 


foi sobre à ausência de equipamentos de proteção individuais, como: capacete, bota, 
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óculos e cinto. Porém, a maior dificuldade foi em decorrência da ausência de tábuas para 
deslocamento nos andaimes, diante dessa situação, os alunos que estavam realizando a 
montagem tinham que se equilibrar para realizar as atividades em altura, se submetendo 
a risco de acidente. 

A realização da tarefa mostrou-se bastante dinâmica, foi visível a interação entre 
os participantes da aula prática. A montagem da unidade sanitária se mostrou uma 
tarefa além de possível de ser realizada bastante simples, visto que alunos totalmente 
sem experiência na montagem, pela primeira vez que se depararam com o desafio e 
conseguiram desempenhar com êxito. 

Para as próximas turmas que irão realizar a tarefa, pretende-se implementar um 
novo método pré-montagem da unidade sanitária, no qual além de poderem assistir o vídeo 
de montagem, os alunos também tenham acesso ao manual para que dessa forma se 
verifique o tempo em que eles conseguem concluir a montagem. O intuito é promover 
uma montagem com o mínimo de tempo possível, realizada por alunos, ou seja, pessoas 
comuns que não constituem uma mão de obra especializada para a realização da tarefa, 
de forma a mostrar que a execução pode ser desenvolvida por qualquer pessoa e atestar a 
praticidade que o método oferece. 


61 CONSIDERAÇÕES FINAIS 


O acompanhamento da aula prática possibilitou a comprovação da simplicidade da 
montagem do material. O projeto, que envolve alunos da disciplina de Tecnologias das 
Edificações de uma Universidade Federal de Mossoró, em parceria com uma instituição 
pública de ensino médio, mostrou-se bastante eficiente em atender a finalidade para a 
qual foi designado, pois, visa apresentar as técnicas de construção na prática, envolvendo 
os alunos na montagem para que estes despertem um maior interesse para o curso de 
engenharia civil. 

Com a realização da tarefa de montagem, foi possível aplicar vários conhecimentos 
adquiridos ao longo da graduação como por exemplo planejamento, gestão e controle 
de obras, a importância de identificar que sem tais conceitos na realização das diversas 
atribuições às quais o engenheiro civil está sujeito, uma tarefa que poderia ser facilmente 
executada pode gerar tantos transtornos, atraso e retrabalhos em uma obra. 

Duas aulas práticas foram realizadas com apoio dos técnicos da Universidade 
que realizou a pesquisa, e, duas foram realizadas somente pelos estudantes do curso 
de Engenharia civil da referida Universidade. Em todas elas ficou evidente a dificuldade 
de encaixe da viga superior em virtude da altura e da necessidade de ajuste entre os 
quatro pilares. A partir desse estudo a empresa fornecedora realizou testes no protótipo e 
modificou o material dos pilares de ferro para alumínio. 
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RESUMO: A qualidade da água tem uma 
relação intrínseca com a gestão da bacia 
hidrográfica, onde as formas de uso do solo 
serão os parâmetros para definir o grau de 
impacto na área da bacia hidrográfica. Diante 
dessas questões o objetivo geral deste estudo foi 
realizar um estudo acerca na Lagoa do Colosso 
e no decorrer do estudo procurou-se alcançar 
os seguintes objetivos específicos: conceituar 
recursos hídricos; verificar que instrumentos 
legais definem as normas da gestão ambiental 
no rio; e por fim conhecer os instrumentos legais 
que definem as normas da gestão ambiental. 
PALAVRAS-CHAVE: Água. Recursos hídricos. 
Lagoa do Colosso. 


ENVIRONMENTAL DIAGNOSIS OF 
COLOSSO LAGOON, FORTALEZA 
MUNICIPALITY, CEARÁ: SCIENCE 
AND TECHNOLOGY OF THE FEDERAL 
INSTITUTE OF CEARÁ — IFCE 


ABSTRACT: Water quality has an intrinsic 
relationship with river basin management, where 
land use forms will be the parameters for defining 
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the degree of impact in the catchment area. In 
view of these questions the general objective of 
this study was to carry out a study about Colosso 
Lagoon and during the study the following specific 
objectives were sought: conceptualizing water 
resources; verify that legal instruments define the 
norms of environmental management in the river; 
and finally to know the legal instruments that 
define the norms of environmental management. 
KEYWORDS: Water. Water resources. Colosso 
Lagoon. 


11 INTRODUÇÃO 


A questão da gestão de recursos 
hídricos mobiliza a cada ano mais cientistas, 
organizações não governamentais e a 


sociedade civil organizada, explicitando a 
preocupação com a manutenção de um dos 
recursos ambientais que, apesar de renováveis, 
tem sofrido diminuições sucessivas de sua 
potabilidade/qualidade diante do crescimento 
populacional e das atividades econômicas. 
Constitui, portanto, matéria de interesse 
internacional, transcendendo as barreiras entre 
as diversas regiões. 

Sendo a água um dos recursos naturais 
mais intensamente utilizados torna-se nítida a 
necessidade de administrá-la, garantindo seus 
requisitos de qualidade e uma oferta que atenda 
a demanda dos polos industriais, dos grandes 
centros urbanos, zonas de irrigação, bem como 
as necessidades metabólicas do homem e de 


outros seres vivos. 
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A qualidade da água tem uma relação intrínseca com a gestão da bacia hidrográfica, 
onde as formas de uso do solo serão os parâmetros para definir o grau de impacto na área 
da bacia hidrográfica. 

Vale relembrar que os estudos relacionados à análise ambiental têm sempre como 
referencial uma determinada comunidade e sua relação com o contexto espacial em que 
vive e desenvolve suas atividades, podendo ser um país, estado, município, lugarejo, 
ou microbacia hidrográfica etc. Desse modo os estudos ambientais apresentam uma 
preocupação holística de relacionar a complexidade social, cultural e econômica de uma 
sociedade com o processo de apropriação e ocupação do meio ambiente em que vivem 
(ROSS, 2006). 

De acordo com Gonçalves et al. (2013) a Lagoa do Colosso fica localizada no bairro 
Edson Queiroz, na cidade de Fortaleza. A localização geográfica da lagoa é de 0558707W 
latitude e 95822645 de longitude. No entorno da lagoa é possível perceber inúmeras 
interferências humanas, dentre essas intervenções estão: construções civis, vias de acesso 
e também a utilização da lagoa como um meio de entretenimento. 

A Lagoa do Colosso é um local interessante para prática do turismo, quanto para 
moradores de Fortaleza, sobretudo no tocante a turismo ecológico. Nesse espaço da 
capital cearense é possível fugir da agitação que é característica dos dias em Fortaleza. 

A Lagoa do Colosso é um local bastante utilizado na capital para diversos esportes, 
graças as suas calmas águas. Dois desses esportes náuticos são o Wakeboard e o Stand 
Up Paddle e também para balneabilidade e lazer em geral. 

Além dos programas aquáticos para aqueles que gostam de praticar esportes, há 
também outras opções fora d'água, nos quais é possível aproveitar um momento familiar 
ou mesmo passar o tempo com os amigos, visto que muitos já adotaram o local como ponto 
de encontro. 

Diante dessas questões o objetivo geral deste estudo foi realizar uma análise na 
Lagoa do Colosso e no decorrer do estudo procurou-se alcançar os seguintes objetivos 
específicos: conceituar recursos hídricos e sua área de influência hidrográfica; verificar que 
instrumentos legais definem as normas da gestão ambiental no rio; e por fim conhecer os 
instrumentos legais que definem as normas da gestão ambiental. 


2 | MATERIAS E MÉTODOS 


Este trabalho teve como método a realização de pesquisa exploratória na bacia 
hidrográfica, e criação de um referencial teórico, no qual fez-se uma avaliação geral sobre 
os autores que discutem a temática, com o objetivo de entender a real situação dos recursos 
hídricos e os impactos causados pelo o homem, em seguida com base nos instrumentos 
legais que definem as normas da gestão ambiental na bacia, principalmente utilizando 
a resolução vigente, a CONAMA 357/2005, apresentar uma análise descritiva dos dados 
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coletados na Lagoa do Colosso. 


31 RESULTADOS E DISCUSSÃO 


De acordo com um estudo realizado por Gonçalves et al. (2013) a contaminação da 
àgua é um problema que tem ocorrido de maneira frequente no momento atual, pois ela 
pode ser realizada de várias maneiras, de acordo com os autores essas contaminações 
podem ocorrer das seguintes maneiras: 


* | Contato de produtos químicos com a água; 
* | Despejo inadequado de esgoto domésticos; 
* | Despejo de uma quantidade muito grande de nutrientes. 


Segundo Gonçalves et al. (2013) através dessas situações o processo de 
eutrofização das águas pode ocorrer, e dessa maneira a quantidade de água própria para 
banho vai sofrendo escassez, assim como aquela água que deveria ser utilizada para 
consumo da população. 

Baseado na citada Resolução do CONAMA após a análise de salinidade da Lagoa 
do Colosso, a mesma foi classificada como: Seção |, das águas Doces; Artigo 4º, III — 
Classe 2. (GONÇALVES et al., 2013). 

Gonçalves et al. (2013) afirma que com o passar do tempo os ecossistemas aquáticos 
terminam por incorporar algumas substâncias que tem origem natural, e essas substâncias 
têm a capacidade de trazer alterações no comportamento químico da água. O que também 
vem a ocorrer é que a ação humana decorrente das substâncias que são lançadas no copo 
d'água pode vir a trazer problemas bastante significativos no tocante a qualidade da água. 

A seguir ressalta-se uma tabela apresentando os parâmetros dos estudos relativos 
a água na Lagoa do Colosso por Gonçalves et al. (2013). 


E Ma JT A LM RI] had = ao BEN — q 


(PH) 
Oxigênio Dissolvido (OD) | 4 a e 3,6 
Salinidado 
Condutividado Elótrica | Ss! Í 
Temperatura 831 l à, =: 


Turbidez 13,9 NTU 
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Figura 4 - Parâmetros analisados 


Fonte: Gonçalves et al. (2013) 
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As temperaturas encontradas no fundo, no meio e na superfície foram 
respectivamente de 28,31ºC, 29,62ºC e 29,58ºC. De acordo com o CONAMA 357/2005, a 
Turbidez encontrada nas águas está dentro do padrão, sendo que nesse padrão se vai até 
40 NTU. 

Gonçalves et al. (2013) relata que os valores encontrados no fundo, meio e 
superfície são respectivamente de 14,2 NTU e 13,5 NTU. O que indica que há pouco 
material particulado em suspensão quando foi realizada a leitura. 

O pH da água deve estar entre 6,0 e 9,0, de acordo com o CONAMA 357/2005. 
Gonçalves et al. (2013) apresentou o resultado de no fundo de 7,23, no meio de 7,23 e no 
fundo de 7,15. O que coloca a lagoa como dentro da Resolução do CONAMA. 

A amostra da água para clorofila “a” foi realizada às 9:40 da manhã e teve adição 
de carbonato de magnésio na sua amostra. Quando analisada a clorofila três medições 
diferentes de onda foram realizadas (663, 645 e 630 nm), o que resulta em 6,487 mg/mº*. 
Tais valores estão dentro do padrão aceito pelo CONAMA. 


De acordo com os resultados obtidos por Gonçalves et al. (2013): 


* || Alagoa do Colosso possui um valor abaixo do aceitável no tocante ao oxigênio 
dissolvido, podendo isso significar que há uma grande quantidade de micror- 
ganismos aeróbios, que são os responsáveis pela decomposição da matéria 
orgânica; 


*- A lagoa conta com a presença de coliformes termotolerantes e E.coli. Com 
isso O uso da água é impossibilitado para a utilização em contato primário, por 
exemplo: agricultura, recreação e a pesca; 


* Quanto a turbidez, esse apresenta-se de acordo com o que é proposto pela 
Resolução do CONAMA 357/2005, estando dessa maneira satisfatório; 

O local onde foi realizada a amostragem apresentou uma grande interferência 
antrópica, de acordo com Gonçalves et al. (2013) isso implica que a população local tem 
realizado a utilização da água do rio para atividades recreativas, plantações e que tendo 
em base essa referência é possível que existam até mesmo o despejo de esgoto ao redor 
do local, podendo esse ser proveniente da população local ou por alguma indústria. 

Assim Gonçalves et al. (2013) concluíram que as análises realizadas não 
apresentaram números satisfatórios baseado na Resolução CONAMA 357/2005 e, que 
dessa maneira com a poluição detectada no local, as águas da lagoa deveriam ser utilizadas 
em atividades que fossem menos exigentes. 


4 | CONCLUSÕES 


O Estado do Ceará possui uma variedade de ambientes e de atrações que fazem 
com que esse seja um dos estados mais visitados do Brasil. Fortaleza conta com 161 
lagoas, rios e lagos, essas são muitas vezes opções de lazer para os finais de semana 
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dos fortalezenses, as atividades que normalmente são realizadas dizem respeito a estar 
com a família ou amigos, aproveitando o ambiente para a realização de piqueniques e de 
pescarias. 

O órgão responsável pela fiscalização das 161 lagoas, rios e lagos de Fortaleza é a 
Secretária Municipal de Urbanismo e Meio Ambiente (Seuma). A limpeza dos reservatórios é 
realizada pela Seuma em parceria com outros órgãos, dentre esses órgãos estão a Defesa 
Civil, Secretarias Regionais e Autarquia de Urbanismo e Paisagismo de Fortaleza (Urbfor). 

A limpeza de lagoas como a do Colosso, Parangaba e do Opaia ocorrem de acordo 
com a Seuma a cada três meses. O Programa Águas da Cidade é gerido pela Seuma,e 
esse programa tem como o seu objetivo a melhora nos índices de balneabilidade dos 
recursos hídricos de Fortaleza. 

A Lagoa do Colosso é uma das que foi considerada própria para o banho, situado no 
bairro Água Fria, assim como o riacho Cachoeirinha e o açude Santo Anastácio, que são 
localizadas no bairro do Pici em Fortaleza. 

Os esgotos são um dos principais problemas relacionados para as lagoas, rios e 
açudes, e eles podem fazer com que o banho venha a se tornar impróprio nesses locais. 

As fiscalizações acerca dos esgotos ocorrem de maneira diária pela Seuma. Essa 
fiscalização pode se dar tanto por denúncias quanto por uma busca ativa do órgão a 
prédios que não possuam ligação à rede pública de esgoto e de ligações clandestinas. A 
balneabilidade é afetada diretamente graças aos resíduos que são despejados no local de 
banho. 
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ABSTRACT: The months of October 2019 until 
the end of January 2020 were very rainy in São 
Jose dos Campos (23º12'S, 45º52'W) tropical 
Brazilian region. In the period from 10/14/2019 
to 01/27/2020, the counting of gamma rays, 
neutrons and rainfall intensity in a tower at 25 
meters high was monitored at every one-minute 
interval. This tower is located in a free area 
without electromagnetic interference from man 
on the site. The average rate of gamma radiation 
count between (0.2-10.0) MeV was 39000 counts 
/ min. The average neutron count observed was 2 
neutrons / min, between 25 eV to 10.0 MevYV. The 
amount of net rainfall in the period was 461 (mm) 
with variations in dry weather, fine, moderate and 
heavy rainfall throughout the monitored interval 
time. It can be seen in these measurements 
that the intense rains correlate very well with 
variations in the gamma rays and delayed of 4 
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to 5 days with the measured neutron intensities. 
This work gives possible explanations about this 
correlation based in in same site of rainfalls, 
gamma and neutron observations. 

KEYWORDS: Neutron, gamma rays and rainfall 
measurements. 


MEDIÇÕES DE GAMA, NÉUTRONS, 
CHUVA E AS POSSÍVEIS CORRELAÇÕES 
NAS REGIÕES TROPICAIS DO BRASIL 


RESUMO: Os meses de outubro de 2019 até o 
final de janeiro de 2020 foram muito chuvosos 
na região tropical brasileira de São José dos 
Campos (23012'S, 45052 W). No período de 
14/10/2019 a 27/01/2020, a contagem de raios 
gama, nêutrons e intensidade das chuvas em uma 
torre a 25 metros de altura foi monitorada a cada 
intervalo de um minuto. Esta torre está localizada 
em área livre, sem interferência eletromagnética 
do homem no local. A taxa média de contagem 
de radiação gama entre (0,2-10,0) MeV foi de 
39.000 contagens / min. A contagem média de 
nêutrons observada foi de 2 nêutrons / min, entre 
25 eV a 10,0 MevV. A quantidade de chuva líquida 
no período foi de 461 (mm) com variações de 
tempo seco, chuva fina, moderada e forte ao 
longo do intervalo de tempo monitorado. Pode- 
se observar nessas medições que as chuvas 
intensas se correlacionam muito bem com as 
variações dos raios gama e retardadas de 4 a 5 
dias com as intensidades de nêutrons medidas. 
Este trabalho dá possíveis explicações sobre 
esta correlação com base em um mesmo local 
de chuvas, observações gama e de nêutrons. 

PALAVRAS-CHAVE: Medições de nêutrons, 
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raios gama e precipitação. 


11 INTRODUCTION 


At the ground level interface of the Earth's surface, ionizing radiation it is composed 
mainly of gamma ray, soil telluric radiation, primary and secondary cosmic ray radiation [1]. 
However, it is difficult to separate over time the intensity of the ionizing radiation emanating 
from each component as the energies overlap. The telluric radiation is given by 2*U, 2ºU, 
*9K and 22Th disintegration's series that are constant for each region. The gamma ray and 
neutrons coming from radon gas arriving through the 2*U in Earth's crust disintegration to 
*eMa and *22Rn reaching the stables isotopes 2!!Pb, “!“Po and ?'*Bi. Radioactive elements 
such as Uranium, Thorium and Potassium are found in almost all types of rocks, sands, 
soils and water [2]. The Radium 2Ra and its decay products are responsible for a major 
fraction of the dose of internal emissions received by humans. *2ºRa has a half-life of 
1,600 years, and decays to Radon ?2Rn, which has a half-life of 3.82 days. The decay 
of 22Rn is followed by successive disintegration of short halí-life alpha, beta and gamma 
ray emitters. After decay stages, the radioactive chain ends with stable lead *ºPb. The 
alpha particles coming from radon gas in interaction with elements of ground level and 
can produce neutrons. With regard to soils and rocks, the *ºRa is present in virtually all 
soils and rocks in varying amounis. Areas with high levels of background radiation found in 
some soils are due to geological conditions and geochemical effects and cause increased 
terrestrial ionizing radiation. Researches in the world, and specifically in Brazil, snow these 
conditions. Several studies report variations throughout the day of radon concentrations. 
Maximum concentrations are observed in the first hours of the day and the lowest values 
are found late near afternoon, when concentrations are about one third of morning values 
[3]. The same profile is observed with the gamma ray intensity variation in the tropics region. 
However, it is likely that variations in concentrations in localities of gamma ray intensity 
are dependent on local meteorological parameters (rain, wind, pressure, temperature and 
cloudiness) in the gamma ray detector site [4]. Electrical discharges in low atmosphere of 
the region also can contribute with production of low energy gamma ray and neutrons near 


ground level according lightning see Figure1 and Figure 2. 
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Fig. 1: Geographical location of Sao Jose dos Campos. 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 8 0 





Fig. 2: Frequency of lightning strikes per km? per year. The frequency of lightning strikes at São José 
dos Campos is about 20-40 strikes /km*year. 


21 METHOD & MATERIAL 


To monitor the gamma radiation in energy interval 200 keV to 10.0 MevV, it has been 
used a portable system detector composed of Sodium lodide scintillator activated with 
Thallium Nal(TI). This crystal (3” x 3”) inches (diameter and height) placed in a thin cylinder 
of aluminum foil and coupled with a PM (photomultiplier) with source power circuit settled in 
1500 VDC and with data acquisition system provided by the company (Aware Electronics- 
Inc., USA) [5]. Detector and associated electronics of gamma ray were previously calibrated 
in ITA (Technological Institute of Aeronautics) laboratory using radioactive sources (Cs- 
137) and (Co-60) in terms of energy from emitted photons 662 keV and 1,17 MeV, 1,33 
MeV respectively [6]. The data acquisition in terms of gamma radiation and intensity of 
rainfall was performed using 1-minute time interval between each measurement. This detail 
contributes to verify possible correlations between variation of rain intensity, and local 
ionizing gamma radiation. 

The Ludlum Model 25311 tube Neutron Detector is designed for detection of thermal 
and fast neutrons (0.025 eV to 10 MeV) [7]. The neutrons are detected, not directly, but 
through nuclear reactions, which result in energetically charged particles such as alpha 
particles. In many instances, intense fields of gamma rays are also found with neutrons. 
Therefore, it is important to choose a method of neutron detection with the ability to 
discriminate against these gamma rays in the detection process. The most common reaction 
used in neutron detection today is: (n+3He>3H+1H + 0.764 MeV) where both the proton H 
and the ,H are detected by gas-filled “H, the Ludlum Model 25311 tube. The helium-3(3He), 
which fills the gas proportional tube of the detector, with 1500 VDC and sensitivity: 100 com/ 
mrem/hr using the neutron source of (241 AmBe fast neutrons) [8]. 

The rainfall intensity in (mm) was measured with a pluviometer (bascule/bucket) 
rain gauge and data logger acquisition developed in ITA according to the international 
recommendations. The data acquisition in terms of ionizing radiation and intensity of rainfall 
was performed using 1-minute time interval between each measurement [9]. Using these 
three tools properly calibrated for the measurements of gamma rays, neutrons and rainfall 
intensity, he was placed in the tower at 25 meters high for simultaneous monitoring see 
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Figure 3. 





Fig. 3 - Aerial and ground view of the tower ACA and his environmental field region in São José dos 
Campos, SP, Brazil (23º 12'45” S, 45º 52'00” W) 


31 RESULTS AND DISCUSSIONS 


During the period from 10/14/2019 to 01/27/2020, 8 peaks of gamma radiation are 
observed, as shown in Figure 4. However, in Figure 5, the measurements of rainfall intensity 
show 5 peaks indicating heavy rains. Therefore, it is admitted that not only heavy rains 
cause an increase in the observed gamma radiation. Figure 6 shows the neutron monitoring 
observed every minute. The intensity during each day varies very related to weather. After 4 
to 5 days of intense rain there is a significant increase in neutrons in the region. This fact can 
be explained by the increase of hydrogen atoms of water in the earth close to the detector 
that through the reaction (protons / alpha) produce neutrons. That kind of cosmogenic 
neutrons produced from cosmic radiation in the Earth's atmosphere or surface is very well 
studied todays [10]. 

With the soil on the wet surface conditions, the back scattering neutrons increase 
at the measurement site. In times of heavy rain with very humid soil, there was always an 
increase in neutron counts caused by interactions of cosmic rays. Also in times of dry and 
very hot weather, the measurements obtained show a variation and increase in neutrons, 
but now originated by the exhalation of radon gas from the earth's surface. Radon gas 
produces alpha particles that interact with the soil and produce neutrons in this energy 
range. For this reason, the neutron monitoring graphic is very variable in this period studied 
here see Figure 5. 
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Fig. 4 - Measures every minute of gamma radiation intensity between 10/14/2019 and 01/27/2020 


showing 8 peaks of radiation increase. 
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Fig. 5 — Measures of rainfalls each minutes during 10/14/2019 to 01/27/2020. 
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Fig. 6 - Detection and monitoring of thermal and fast neutrons (0.025 eV to 10 MeV) during the period of 


10/14/2019 to 01/27/2020. Green line represents 1-day variation. 


Figures 4, 5, 6 show a good correlation between rainfall intensities and intensities of 
gamma rays during monitoring time. In the case of neutron and rainfall measurements, the 
correlation is good but there is also a delay of a few days as clearly seen at the end of the 


monitoring curve in the Figure 6. 


41 CONCLUSION 


Variations in the intensities of gamma rays, neutrons and rains were monitored 
between 10/14/2019 to 1/27/2020 in tropical regions with only 1 minute intervals between 
each measurement. These measures show a good correlation between rain and gamma 
radiation due to the presence of radon gas in heavy rains. The correlation between rain 
and neutrons is due to the wet soil coming from rain and the phenomenon of neutron back 
scattering via cosmic rays. Also in the case of neutrons in dry times is due to exhalation of 
radon gas that increases the production of neutrons via interaction of alpha particles with 
elements of the Earth surface. 
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RESUMO: No presente trabalho, foi desenvolvida 
a Simulação referente a um Reator Contínuo 
de Tanque Agitado (CSTR), ao qual foram 
aplicadas perturbações de diferentes magnitudes 
nas variáveis de entrada. Observou-se o 
comportamento transiente do sistema e efetuou- 
se a sintonia de um Controlador Pl. Em seguida, 
aplicaram-se ao processo elementos dinâmicos 
denominados desacopladores, os quais visam à 
atenuação da interferência entre as malhas em 
um sistema de controle automático. Executado o 
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controle multivariável através do desacoplamento 
entre as malhas, das variáveis temperatura e 
concentração de reagente, foram verificadas 
as respostas do processo e a obtenção de 
resultados de precisão significativa utilizando um 
software de simulação gratuito. Concluiu-se que, 
em geral, todos os controladores apresentaram 
desempenho satisfatório, excetuando aquele 
com as malhas totalmente desacopladas. Logo, 
para a operação analisada, há dependência 
das duas variáveis manipuladas, corrente de 
alimentação e calor retirado, para o ajuste preciso 
das variáveis controladas. 

PALAVRAS-CHAVE: Simulação, Reator CSTR, 
Controle, Desacoplamento. 


APPLICATION OF DECOUPLING 
BETWEEN CONTROL MEDS IN 
CONTINUOUS STIRRED TANK REACTOR 
(CSTR) 


ABSTRACT: In the present work, a simulation 
of a Continuous Agitated Tank Reactor (CSTR) 
was developed, to which perturbations of different 
magnitudes were applied to the input variables. 
Thetransientbehavior ofthe system was observed 
and a PI Controller was tuned. Then, dynamic 
elements called decouplers were applied to the 
process, which aim to attenuate the interference 
between the meshes in an automatic control 
system. After performing the multivariable control 
through the decoupling between the meshes, the 
temperature and reagent concentration variables, 
the responses of the process and the obtaining of 
results of significant precision were verified using 
a free simulation software. It was concluded that, 
in general, all controllers presented satisfactory 
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performance, except for the one with fully uncoupled meshes. Therefore, for the analyzed 
operation, there is dependence of the two manipulated variables, supply current and removed 
heat, for the precise adjustment of the controlled variables. 

KEYWORDS: Simulation, CSTR reactor, Control, Decoupling. 


11 INTRODUÇÃO 


Na sociedade atual com o aumento populacional e a incessante demanda por 
produtos de alta qualidade, é essencial a inserção de profissionais com a qualificação 
necessária para a implementação desses requisitos. Desta forma, com a visão global e 
acurada dos processos químicos pelos profissionais, se desenvolve uma análise crítica 
sobre a utilização de insumos, eficiência na integração energética, mitigação de impactos 
ambientais e certificação de produtos finais dentro das normas do mais elevado rigor. 

No contexto abordado acima, a criação de metodologias de melhoria contínua e 
otimização de processos é fundamental. Um dos métodos associados ao estabelecimento 
de práticas eficientes é o controle de processos. Este visa garantir que uma propriedade do 
sistema (nível, temperatura ou concentração) fique em um valor previamente estabelecido. 
Assim, pode-se garantir que o processo apresentará menores desperdícios de recursos e 
aplicar diretrizes à otimização do lucro, redução de custeios operacionais e de tempo, além 
da menor geração de rejeitos. 

Nesse viés, a simulação do processo ganha contornos fundamentais. A partir 
do artifício computacional, pode-se visualizar o comportamento de plantas industriais 
completas ou operações específicas. É inteligente, e até prudente, a verificação preliminar 
dos possíveis atenuantes à robustez dos diferentes processos químicos, e também o 
estabelecimento de alternativas na escolha das variáveis para o controle automático de 
processos. 

No presente trabalho é efetuado o controle multivariável com desacoplamento 
entre as malhas do processo referente ao Reator Contínuo de Tanque Agitado (CSTR). 
Ademais, é apresentada uma estrutura completa através da modelagem matemática e 
desenvolvimento das estruturas em malha aberta, malha fechada e com desacoplamento 
do referido sistema, o que possibilita uma visão amplificada dos processos intrínsecos ao 
reator. 

Objetivou-se a aferição da precisão e robustez do sistema de controle, frente as 
perturbações nas variáveis de entrada, através de indicadores de desempenho e o esforço 
de controle dos equipamentos eletromecânicos, visando à condição para a eficaz execução 


da operação. 


21 REFERENCIAL TEÓRICO 


Sistemas SISO (single-input, single-output) são aqueles que apresentam uma 
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variável manipulada (entrada), também chamada de variável de controle, e uma variável 
controlada (saída), podendo ser denominada, também, de variável de processo. Já os 
sistemas MIMO (multiple-input, multiple-output) são constituídos por duas ou mais variáveis 
manipuladas as quais visam ajustar múltiplas variáveis controladas. 

Os processos químicos, em sua maioria, são sistemas multivariáveis. 
Consequentemente, é necessário implementar várias malhas de controle em um mesmo 
equipamento. Devido à interferência de uma variável manipulada em mais de uma variável 
controlada, as malhas de controle interagem entre si, resultando em um controle de baixo 
desempenho (KALID, 2004; STEPHANOPOULOS, 1983). Logo, o uso de desacopladores 
se torna essencial ao pleno desenvolvimento dos processos mencionados. 

A Figura 1 A ilustra um sistema com duas variáveis de processo e duas variáveis 


manipuladas, sem a aplicação do controle automático (conhecido como processo em malha 
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Figura 1- A-Diagrama de blocos para um processo em malha aberta; B-Diagrama de blocos para um 
processo em malha fechada 


Fonte: adaptado de Stephanopoulos (1983) 


em que: H(s) representam as quatro funções de transferência intrínsecas ao 
processo; m(s) denota a variável manipulada e Y(s) a variável de processo. 

A Figura 1 B apresenta o mesmo processo da ilustração anterior, no entanto, há a 
presença de dois controladores (Gc, e Gc,), caracterizando um sistema em malha fechada. 
O sistema de controle aplica uma ação corretiva na variável controlada, de modo a esta 
apresentar um valor equivalente ou próximo à definição do set-point (Y,sp)- 

Nessa linha, considerando que, inicialmente, ambas as saídas estejam nos valores 
de set-point, ao se efetuar uma perturbação neste causada pela malha de controle 2, tal 
alteração acarretará uma mudança no valor de m,. Não obstante, também será gerada uma 
perturbação indesejada na malha 1. Desta forma, é crucial a alteração no valor de m, de 
modo a cancelar o efeito da interação de m,. 

Deve-se, então, introduzir um elemento dinâmico com uma função de transferência 


que gere uma ação compensativa para atenuar os efeitos de acoplamento entre as 
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variáveis trabalhadas. Os seguintes equacionamentos abordam de forma mais clara o que 
fora enunciado previamente. 
—Hiz2 


D'(s) = H 
11 





(1) 


Esse elemento dinâmico é chamado desacoplador (D) e, quando instalado no 
sistema de controle minimiza qualquer efeito que a malha 2 possa ter sob a malha 1. 

Para a atenuação da interação da malha 1 na 2, segue-se o mesmo raciocínio 
acima, chegando à função de transferência mostrada na Equação 11: 


D'“(s) = E (2) 





Em um processo com duas entradas e duas saídas, são requisitados dois 
desacopladores para gerar malhas não interativas. Assim, têm-se ações compensativas 
para minimizar ou até mesmo cancelar por completo as interferências entre as malhas 
de controle. A seguir é apresentado um sistema contendo dois elementos para o 
desacoplamento entre as malhas. 


Y1 


Y2,sp TZ 





Figura 2- Sistema Contendo dois elementos 


Fonte: adaptado de Stephanopoulos (1983) 


31 METODOLOGIA 


A metodologia se pautou na modelagem e simulação da reação química em 
Reator CSTR, juntamente com a aplicação de controladores com as malhas acopladas, 
desacopladas e parcialmente desacopladas. Também serão feitos os comentários 
pertinentes às respostas, com base nas perturbações aplicadas às variáveis de entrada, 


dos controladores citados. 
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3.1 Modelagem Matemática 


O reator contínuo de tanque agitado (CSTR) é comumente usado em processos 
industriais, principalmente, para reações em fase líquida. Considera-se que os constituintes 
do meio reacional estão perfeitamente misturados, ou seja, temperatura e concentração 
são idênticas em qualquer ponto do reator (FOGLER, 2009; LUYBEN, 2007). 

Realizou-se a modelagem matemática baseada em um problema prático abordado 
em Fogler (2009). A reação química considerada foi a isomerização em fase líquida no 
n-butano, em um Reator CSTR, para geração do isobutano, o qual é um produto valioso 
utilizado na fabricação de aditivos para a gasolina (FOGLER, 2009). 

Considerou-se que o volume do meio reacional (V) se manteve constante ao longo 
do processo, portanto, a vazão volumétrica de entrada (F.) é igual à vazão de saída (F). 
Além disso, é destacável que a reação química abordada era exotérmica, de primeira ordem 
e elementar, do tipo: A — B. A partir dessas considerações desenvolveu-se o balanço para 
a concentração do reagente A e o balanço energético do sistema. 


3.2 Funções de transferência e parâmetros do controlador 


Com base nas formulações matemáticas, implementou-se na plataforma x-cos os 
diagramas de blocos para a resolução das equações diferenciais. 

O reator em questão apresenta duas variáveis manipuladas, vazão volumétrica e 
quantidade de calor retirada por um sistema de resfriamento. Foram exercidas perturbações 
na interface de simulação, sendo estas perturbações do tipo degrau em +20% no valor das 
variáveis citadas. 

O objetivo dessas flutuações é encontrar a função de transferência do processo. 
Para tal, utilizou-se a ferramenta Solver do Excel, de modo a minimizar o erro da diferença 
entre a curva do processo obtida e a curva de identificação, encontrando, assim, os 
parâmetros da função de transferência. 

Logo após, utilizou-se do método de Astrom- Hagglund para realizar as perturbações 
inerentes ao processo e obter os parâmetros do controlador. O método é utilizado em malha 
fechada, no qual se utiliza um relé que provoca o aparecimento de oscilações e gera-se um 
gráfico para avaliação do Kc crítico. 


3.3 Perturbações e controle com e sem desacoplamento 


Aplicaram-se perturbações, de diferentes magnitudes, às variáveis de entrada 
(Ca, e T,), a fim de se observar como o processo em malha fechada se comportava. As 
flutuações em questão foram de 10% em Ca, e 1% em T,, embasadas em uma ordem de 
grandeza coerente. Realizou-se a verificação do desempenho dos controladores, com as 
malhas acopladas e desacopladas, através dos gráficos gerados. 

Por fim, foi aferido o Esforço de Controle (Esf) e o indicador de desempenho ITAE 
(Integral do Erro Absoluto Ponderado pelo Tempo), a partir das simulações desenvolvidas. 
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As Equações 3 e 4 representam de forma elucidativa o cálculo do ITAE e Esf, 


respectivamente. 


ITAE = S, tle(t)ldt (3) 


Esf = J, Audt (4) 


em que t representa o tempo de simulação ou de operação e u é a variável 


manipulada. 


41 RESULTADOS 


4.1 Modelagem Matemática 


Mediante as considerações desenvolvidas a respeito da manutenção do volume do 
meio constante, a Equação final representativa da variação da concentração ao longo do 
tempo é mostrada pela expressão matemática a seguir: 

ds 


F -Ea 
de - 7 (Cao — Ca)— ko.€ RT.Ca (3) 


em que: F representa a vazão volumétrica do processo, V o volume do meio 


reacional, C 


ao à concentração inicial do reagente, C| a concentração do reagente em 


determinado tempo, k, denota o fator pré-exponencial relacionado à frequência de colisões, 
Ea é a energia de ativação da reação, T representa a temperatura, a qual ocorre a reação, 
e Ra constante geral dos gases. 

A Equação descritiva da variação, ao longo do tempo, da temperatura (T) é mostrada 

abaixo: 
| fui ESA 
Epic Ep (6) 

em que: T, é a temperatura inicial, Q é a quantidade de calor retirada do sistema, p 
é a densidade do reagente, Cp é o calor específico da espécie química envolvida e AHr 
representa a variação de entalpia da reação. 

Logo após a modelagem matemática, desenvolveram-se no x-cos os diagramas de 
blocos referentes às equações diferenciais (as quais denotam o comportamento transiente 
do sistema) apresentadas pelas Equações 5 e 6. 

As Figuras 3 Ae 3 B ilustram as resoluções das equações diferenciais inerente à 
variação de Concentração do reagente analisado e Temperatura, ao longo do tempo, pelos 


diagramas de blocos. 
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Figura 3- A- Variação de concentração do reagente com o tempo; B- Variação da temperatura com o 
tempo 


Fonte: acervo do autor 


O uso da plataforma x-cos é respaldado pela facilidade na disposição dos elementos 
essenciais à averiguação do perfil dinâmico do processo e ser bastante intuitivo. Os valores 
dos parâmetros e demais variáveis estão alocados na tabela seguinte. É importante salientar 
que alguns valores de variáveis de processo como: vazão volumétrica, calor removido e 


das condições iniciais de temperatura e concentração foram estabelecidos pelo autor. 


V (mº) Fo (m*/s) Q (W) Ca (mol/m*) T, (K) 
1,0 0,0044 30000 1000 330 
Ea (J/mol) p(Kg/m?) Cp (J/Kg.K) AHr (J/mol) Ko 
65700 5783 2600 -6900 3,50.107 


Tabela 1- Valores das propriedades referentes às equações diferenciais 


Fonte: adaptado de Fogler (2009) 


Além dos dispostos acima, utilizou-se o valor da constante dos gases, R, equivalente 
a 8,314 J/mol.K (CENGEL; CIMBALA, 2012). 


4.2 Funções de transferência e parâmetros do controlador 


As Figuras 4 A e 4B representam as curvas de resposta frente à perturbação 
aplicada. Vale destacar que a representação é feita da concentração do reagente A (Ca) 
e da Temperatura (T ), em variável desvio (diferença do valor atual da propriedade com 
relação ao estado estacionário) em função do tempo (t). Empregam-se as variáveis de 
desvio, com o objetivo da eliminação das condições iniciais dos modelos em função de 
transferência (FT), cuja vantagem é a nulidade adquirida pelo termo de condição inicial 
(BOJORGE, 2021). É interessante abordar que os valores encontrados para o estado 
estacionário envolvendo a concentração de A e temperatura foram, respectivamente, 3851 
mol/mº e 354 K. 
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Figura 4- A-Perturbação aplicada à vazão e verificação do comportamento da concentração; B- 
Perturbação aplicada à vazão e verificação do comportamento da temperatura 


Fonte: acervo do autor 


Os parâmetros intrínsecos ao sistema estão dispostos na tabela seguinte. 


Vazão (m/s) Calor (W) 
Kp 1,360.10º 5455 
T 256,2 269,0 


Tabela 2- Valores de Tt e Kp para a identificação do sistema mediante alteração da vazão 


Fonte: acervo do autor 


Logo, a função de transferência global para a concentração do reagente A é: 


1,360.10º. F" (7) 
256,25 + 1 


No tocante à modificação do calor retirado, as figuras abaixo ilustram o comportamento 
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Figura 5- A-Perturbação aplicada ao calor e verificação do comportamento da concentração; B- 
Perturbação aplicada ao calor e verificação do comportamento da temperatura 


Fonte: acervo do autor 
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Ao se analisar as figuras acima, observa-se que o comportamento de ambas é 
inerente ao sistema de segunda ordem, sendo assim, podem ser aproximadas ao modelo 
de primeira ordem com tempo morto. Efetuando os ajustes, chega-se aos seguintes valores 
dispostos na Tabela 3. 


Vazão (m/s) Calor (W) 

Kp 0,080 5,00.10 
T 482,4 452,0 
td 117,7 39,4 


Tabela 3- Valores de T, Kpe td para a identificação do sistema mediante alteração do calor 


Fonte: acervo do autor 


A FT global para a temperatura referente à alteração da variável manipulada calor é: 


500.40" O (8) 


T(S) = 
(s) 452,05 + 1 


Conseguidas as funções de transferência, foram calculados os parâmetros do 
Controlador Pl a partir do método de Astrom-Hagglund. Uma observação pertinente sobre 
a utilização deste método foi pelo fato de não apresentar restrições. Resumidamente, na 
Tabela 4 estão alocados os valores correspondentes aos parâmetros dos controladores Pl. 


Controladores 


Concentração Temperatura 
Kc TI Kc TI 
5,04.107 1041,7 1720 1125 


Tabela 4- Valores dos parâmetros dos controladores Pl 


Fonte: acervo do autor 


4.3 Perturbações e controle com e sem desacoplamento 


Com o auxílio fornecido pela sintonia dos controladores e os equacionamentos 
matemáticos representativos das variações de Ca e T, construiu-se o sistema completo 
envolvendo diagramas de blocos. A Figura 6 representa a síntese do trabalho desenvolvido. 
Nela se encontram os dois desacopladores e todo o aporte necessário para a geração 
dos gráficos, os quais serão ferramentas importantes para a comparação da performance 


envolvendo os controladores. 
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Figura 6- Sistema com desacoplamento aplicado no x-cos 


Fonte: acervo do autor 


Os gráficos representativos do comportamento do sistema com os controladores 
acoplados e desacoplados estão alocados a seguir, para cada uma das perturbações 
aplicadas nas variáveis de entrada, concentração (Ca,) e temperatura (T,). Ademais, são 


tabelados os valores conseguidos de indicadores de desempenho e esforços de controle. 


Perturbação de +10% Ca, e +1%T, 


O comportamento transiente da concentração é mostrado pelo gráfico seguinte. 
Nele, são apresentadas as variações do sistema com as malhas acopladas (SA), com as 
malhas desacopladas (SD) e com as malhas parcialmente desacopladas; D1 se refere ao 
sistema desacoplado parcialmente, de modo a atenuar a interação da malha de controle de 
concentração na malha referente ao controle de temperatura; D2 também representa um 
sistema desacoplado parcialmente, o qual visa a não interferência da malha de controle da 


temperatura na malha de controle da concentração. 
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Figura 7- A-Sistema em malha fechada correspondente à variação da concentração ao longo do tempo 
para a perturbação de +10% Ca e +1%T,; B- Variação da vazão ao longo do tempo para a perturbação 
de +10% Ca e +1%T, 


Fonte: acervo do autor 
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Ao se elevar a concentração inicial, a conversão aumenta, visto que é maior a 
probabilidade de colisões efetivas com o aumento no número de partículas reagentes 
(Brady; Humiston, 1986). Para o sistema ser ajustado ao set-point, a vazão é alterada de 
modo a diminuir (aumentar o tempo espacial, elemento representativo do tempo em que 
as moléculas levam para percorrer o elemento espacial do reator) e, desta forma, elevar a 
conversão para que a concentração de A volte ao valor estabelecido inicialmente. 

NA Figura 7 B, observa-se que há um aumento significativo da vazão para o 
sistema com as malhas desacopladas e com desacoplamento parcial, D1. Este fato implica 
diretamente a geração de um ajuste ruim, o que pode ser comprovado pela análise da 
resposta para a variável controlada (ambos os controladores citados apresentaram maiores 
valores de off-set ao se avaliar a diferença entre o set-point e o valor de concentração no 
intervalo de tempo analisado). 

Na Figura 8 A, é mostrado o comportamento transiente da temperatura no interior 
do reator. O gráfico denotativo à retirada de calor no decorrer do tempo está alocado na 
Figura 8 B. 


3554 
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Figura 8- A- Sistema em malha fechada correspondente à variação da temperatura ao longo do tempo 
para a perturbação de +10% Ca, e +1%T,; B- Variação do calor ao longo do tempo para a perturbação 
de +10% Ca e +1%T, 


Fonte: acervo do autor 


Ao se elevar a temperatura, a conversão sofre um aumento pelo fato decorrente do 
acréscimo de energia cinética das moléculas, acarretando uma maior facilidade em romper 
a barreira reacional imposta pela energia de ativação (Brady; Humiston, 1986). Para efetuar 
um controle de temperatura, a quantidade de calor retirada deve aumentar (fato que foi 
observado no gráfico acima para a maioria dos controladores). Todavia, para o sistema 
totalmente desacoplado, houve um declínio considerável a partir de 1500s, de tal modo em 
que a temperatura não é controlada de forma acurada. O fato descrito, não esperado, é 
decorrente de uma ineficiência do controle da temperatura somente pela manipulação da 


quantidade de calor. 
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Na Tabela 5 são apresentados os valores dos indicadores de desempenho para a 
concentração (ITAE Ca) e temperatura (ITAE T). Além disso, estão alocados os esforços de 
controle para a vazão (Esf. T) e sistema de resfriamento (Esf. Q). 


SA SD D1 D2 
ITAE Ca 1,06.10'º 4,1.1010 1,31.10'º 1,69.10'º 
ITAE T 7,89.107 6,84.10º 4,0.108 2,56.108 
Esf. F 7,6714101? 1,69.10" 9,20.10!2 5,91.10!2 
Esf. Q 8,/9.10º 6,17.10º 5,63.10º 1,03.10º 


Tabela 5- Indicadores de desempenho e esforços de controle 


Fonte: acervo do autor 


Com base na tabela anterior, constata-se o melhor desempenho, no sistema 
acoplado, visto que este apresenta os menores valores de ITAE, parâmetro associado ao 
erro. Nota-se que os esforços inerentes à vazão apresentaram valores pequenos, isto é, 
não há um efeito pronunciado no equipamento transportador de fluido ao se efetuar o ajuste 
da vazão. Quanto ao calor, o menor esforço de controle é conseguido pelo controlador 
parcialmente desacoplado D1, para o qual é requisitado um menor gasto com refrigeração 


para a execução do ajuste que se mostrou robusto. 


Perturbação de -10% Ca e 1%T, 


Para a referente perturbação foram gerados os gráficos relativos do comportamento 
ao longo do tempo para as variáveis controladas e manipuladas. Na Figura 9 A, é mostrado 
o perfil dinâmico da concentração. No presente caso, ao se reduzir a concentração inicial, 
a conversão é reduzida e, por conseguinte, a vazão necessita ser aumentada, objetivando 
a elevação da concentração do reagente ocasionada pela redução da conversão reacional. 
A Figura 9 B ilustra o comportamento transiente da vazão. 
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Figura 9- A- Sistema em malha fechada correspondente à variação da concentração ao longo do 
tempo para a perturbação de -10% Ca, e -1%T,; B- Variação da vazão ao longo do tempo para a 
perturbação de -10% Ca, e -1%T, 


Fonte: acervo do autor 


Mais uma vez, nota-se uma incongruência para o controlador com as duas malhas 
desacopladas, a vazão diminui de forma abrupta em torno de 1500s e, no tempo aproximado 
de 6200s, fica estagnada em um valor extremamente baixo, o que não permite a execução 
de um controle significativo. 

De forma posterior, foram gerados os gráficos referentes ao controle de temperatura 
e do calor extraído do sistema. 
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Figura 10- A- Sistema em malha fechada correspondente à variação da temperatura ao longo do tempo 
para a perturbação de -10% Ca, e -1%T,; B- Variação do calor ao longo do tempo para a perturbação 
de -10% Ca, e -1%T, 


Fonte: acervo do autor. 


Observa-se que a variável de estado temperatura, apresentou valores destoantes 
ao set-point, para todos os controladores, em especial, para SD. Com o decréscimo em 
To deve-se reduzir a quantidade de calor retirado para que a temperatura da meio volte 
a subir. Com base na última afirmação, o comportamento do calor para SD é absurdo, 
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visto que a retirada de energia aumenta de maneira significativa, acarretando um controle 
ruim da temperatura. Pode-se afirmar, então, para o sistema em estudo são necessárias, 
obrigatoriamente, duas variáveis manipuladas para o ajuste acurado das variáveis 


controladas. 


SA SD D1 D2 
ITAE Ca 2,46.10º 6,9.1010 1,8.10º 1,24.101º 
ITAE T 1,75.108 1,02.10º 2,07.108 2,31.108 
Esf. F 4,56.10!2 1,83.10" 5,22.10!2 6,01.1072 
Esf. Q 6,88.10º 6,70.10º 5,35.10º 9,89.10º 


Tabela 6- Indicadores de desempenho e esforços de controle 


Fonte: acervo do autor 


Mediante a observação da Tabela 6, verifica-se que o controlador para o sistema 
totalmente desacoplado, apresentou o pior desempenho, analisando os valores de ITAE. 
Esse fato era esperado, visto que o ajuste ao set-point, confirmado pela avaliação dos 
gráficos, foi extremamente defasado. Para D1, foi obtido o menor esforço de controle para 
o sistema de resfriamento. 


5 | CONCLUSÃO 


Com o aporte teórico fornecido pela literatura foi possível a aplicação da estratégia 
de controle multivariável através do desacoplamento entre as malhas de controle para 
a concentração e temperatura de um Reator CSTR, envolvendo uma reação química 
de grande relevância nos processos industriais, retirada de um problema prático, para o 
qual engenheiros de processos são solicitados, com o objetivo de melhorar o processo 
químico descrito. Como era esperado, obtiveram-se respostas diferentes para cada tipo de 
controlador, mediante as flutuações aplicadas às variáveis de entrada. 

Para o esforço da manipulação do calor retirado, os menores valores foram 
encontrados para o caso de desacoplamento parcial D1. Consequentemente, este é o mais 
indicado visando à economia de fluido refrigerante para a aplicação no processo analisado. 

Uma menção a ser feita é que fica a cargo de projetistas a análise melhor definida 
de quais são as perturbações mais comuns nesse tipo de reator e a seleção do melhor 
controlador com base em economia de recursos financeiros e adequações a condições 
estabelecidas. 
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RESUMO: Os profissionais que trabalham em 
regimes de turno, como na mineração, têm o 
seu ciclo circadiano bastante alterado, podendo 
ser conduzidos a um estado de fadiga. Nessa 
condição há uma maior propensão à tomada de 
decisões erradas ou tardias, que comumente 
são causas de acidentes. Nesse contexto, este 
trabalho visa verificar a efetividade da aplicação 
do sistema de gerenciamento de fadiga 
OptAlertê, que mede a sonolência de operadores 
de caminhão fora de estrada, em uma mina a 
céu aberto. Três fases de teste, sob distintas 
condições de tratativa de fadiga, durante o 
período de seis semanas cada foram realizadas 
com vinte e quatro operadores. Os resultados 
mostraram que a implantação de um sistema de 
gerenciamento capaz de quantificar a sonolência 
de operadores e estabelecer procedimentos e 
tratativas específicas é efetiva no aumento do 
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nível de alerta dos operadores, podendo assim 
ser útil para redução do risco de acidentes. 
PALAVRAS-CHAVE: Gerenciamento de fadiga; 
sonolência; caminhão fora de estrada; risco de 
acidente. 


FATIGUE MANAGEMENT OF TRUCK 
OPERATORS WHEN OFF WORK : CASE 
STUDY IN AN OUTSIDE MINE 


ABSTRACT: Shift workers, as in mining, have 
their circadian cycle greatly altered, and may 
be led to a state of fatigue. In this condition, 
there is a greater propensity to make wrong or 
late decisions, which are commonly causes 
of accidents. In this context, this work aims to 
verify the effectiveness of the application of the 
OptAlertê fatigue management system, which 
measures the drowsiness of off-road truck 
operators, in an open pit mine. Three test phases, 
under different fatigue management conditions, 
during a period of six weeks each were performed 
with twenty-four operators. The results have 
shown that the implementation of a management 
system capable of quantifying the sleepiness of 
operators and establishing specific procedures 
and treatments is effective in increasing the 
operator's alertness level, thus being able to be 
useful for reducing the risk of accidents. 
KEYWORDS: Fatigue management; drowsiness, 
off-road truck, risk of accident. 


INTRODUÇÃO 


A fadiga representa um grande desafio 
no seguimento da mineração. Dependendo 
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da ocupação específica e da configuração operacional, os trabalhos do setor podem ter 
uma quantidade razoável de tarefas intensas, monótonas e repetitivas. Combinado com 
as cumprimento da jornada de trabalho e horários de turno, a aparição de fadiga nos 
colaboradores é algo esperado, como descrito por Bauerle, Dugdale e Poplin [1] 

É prática comum que operações mineiras adotem regimes de trabalho em turnos, de 
maneira a atender as demandas de mercado, sempre em acordo com as regulamentações 
trabalhistas locais. Os regimes de trabalho adquirem diversas formas, desde a adoção de 
sistema fly-in/fly-out aos regimes de turno em diferentes configurações. Boivin e Boudreau 
[2] estudaram os impactos na saúde de trabalhadores e os distúrbios do sono associados 
ao trabalho em regimes de turno. Os autores concluíram que o trabalho sob estes regimes 
aumenta o risco de erros e acidentes no ambiente de trabalho. 

A fadiga é um fenômeno que interessa várias áreas de estudos da sociedade em 
geral. Por isso, não se obtêm uma unanimidade de definição em seu conceito. Para este 
estudo se adotará a definição de Everitt [3], que diz que a fadiga pode ser inferida por 
meio de várias características observáveis e monitoráveis. Segundo Shen, Barbera e 
Shapiro [4], a fadiga se manifesta, principalmente, por meio de quatro fatores: psicológico 
(cansaço, falta de motivação, ações induzidas pelo estresse), fisiológico (perda de força e 
resistência, consumo de energia), cognitivo (tempo de reação retardado, esquecimento) e 
comportamental (fechamento da pálpebra ou aceno de cabeça, fala mais lenta, diminuição 
da produtividade). 

Um estudo publicado em 2007 pela Caterpillar Global Mining [5] estimou que até 65% 
dos acidentes de transporte de materiais em minerações a céu aberto foram relacionados à 
fadiga dos operadores. Pelders e Nelson [6] investigaram os fatores chave que contribuem 
para a fadiga em operadores em quatro minas de ouro e em uma mina de platina na África 
do Sul. Qualitativamente, a fadiga pôde ser atribuída ao regime de trabalho de turno, bem 
como às questões cotidianas, como problemas pessoais de cada indivíduo. 

Johns et al. [8] propuseram uma escala (JDS ou Johns Drowsiness Scale) derivada 
de uma análise de regressão múltipla stepwise backward, validada estatisticamente, .para 
medição de sonolência através de uma combinação ponderada de variáveis oculares por 
oculografia de refletância de infravermelho. Essa escala foi proposta através de um estudo 
que consistiu na análise de falhas de desempenho em dois diferentes testes de tempos de 
reação (tempo de reação e tempo de decisão) e também em testes de condução simulada 
de veículos. Trinta e um voluntários saudáveis participaram dos dois testes de tempo 
de reação por quinze minutos, com e sem privação de sono (27 - 33 horas). Em ambos 
cenários, a sonolência deles foi medida pela escala JDS (0 — 10) a cada minuto. Como 
resultado, o experimento concluiu que após a privação de sono o tempo de reação cresceu 
e os erros de omissão (falha em responder ao estímulo em até 2 segundos) ocorreram 
com mais frequência nos dois tipos de teste. Saídas de pista também aumentaram 
consideravelmente de acordo com o aumento da pontuação da escala e o risco por 
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minuto cresceu progressivamente. Portanto, o risco da falha em estado de sonolência 
pode ser amplamente associado às altas pontuações na escala JDS, refletindo mudanças 
nas características de movimentação de pálpebra durante as piscadas. Foram ainda 
identificadas três faixas de interesse na escala JDS, a saber: o risco de falha começava 
a aumentar apenas levemente à medida que a pontuação ultrapassa o valor de 3 (risco 
baixo); para pontuações de 4,5 até 5 (risco médio) o risco aumentava substancialmente 
e acima de 5 (risco alto) o risco de falha de desempenho era tão alto que o indivíduo era 
instruído que estava sonolento demais para continuar a dirigir. 

Em um experimento conduzido por Aidman et al. [8] foi realizado o monitoramento 
de condutores de veículos leves utilizando o sistema de monitoramento de fadiga OptAleri 
em rodovias, que usa a escala JDS. Esse grupo foi supervisionado por alguns dias com o 
sistema funcionando em dois cenários: no primeiro cenário o sistema apenas monitorava e 
não gerava nenhum alerta de sonolência ao condutor. No segundo cenário, os condutores 
eram comunicados sempre que o sistema alertava e gerava alarme de sonolência. Os 
resultados obtidos foram expressivos e verificados através das variações do JDS. Deste 
modo, verificou-se que no segundo cenário as pontuações de JDS foram menores do que 
no primeiro. Consequentemente, verificou-se que o sistema corroborou para um maior nível 
de alertar dos condutores. 

Nesse contexto, este trabalho visa verificar a efetividade da aplicação do sistema 
de gerenciamento de fadiga OptAlert no aumento do nível de alerta dos operadores de 
caminhão fora de estrada em uma mina a céu aberto, contribuindo, consequentemente, na 


redução do risco de acidentes relacionados à fadiga. 


METODOLOGIA 


Foram selecionados aleatoriamente vinte e quatro operadores de caminhões fora de 
estrada, que trabalham em regime de turno, para avaliação do sistema de monitoramento 
de fadiga. A mina em estudo opera vinte e quatro horas por dia, sete dias por semana, com 
os turnos configurados da seguinte maneira: 


* Noturno: 01:00 às 07:00h; 
. Diurno: 07:00 às 16:00h; 
* | Tarde: 16:00 às 01:00; 


Durante o experimento os operadores trabalharam em escala fixa em seis caminhões, 
no quais foram instalados os kits do sistema OptAlertê. Os operadores foram avisados do 
início do experimento e treinados no funcionamento do sistema. 

O período de avaliação foi de cinco meses, superando a ordem de onze mil horas 
monitoradas. O experimento foi dividido em três fases de seis semanas cada, conforme 
Figura 1. 
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Fase 2 Fase 3 


Central de 


Intervenções na tela monitoramento 





Figura 1: Fases do experimento. 


* Fase 1: os operadores foram orientados a operar os equipamentos utilizando 
os óculos OptAlert?), mas sem receber nenhum tipo de aviso ou interferência 
sobre o estado de sonolência medido. A tela do dispositivo instalado na cabine 
do equipamento foi desabilitada, não informando o operador sobre os níveis 
de sonolência apresentados. Essa fase, chamada de “modo cego”, teve como 
objetivo estabelecer o padrão de alerta e sonolência da mina, de maneira a co- 
nhecer os níveis de sonolência da operação previamente à implementação do 
sistema, possibilitando análises e comparações futuras. 


* | Fase2:os operadores passaram a receber, em tempo real, alertas visuais e so- 
noros através da tela do dispositivo instalada dentro da cabine do equipamento 
quando atingiam níveis de sonolência médios e altos. Nessa fase objetivou-se 
estabelecer o grau de eficiência do sistema, de maneira gradual, assim como 
a mudança comportamental dos operadores quando sujeitos a intervenções do 
sistema 


* Fase 3: nessa fase, além dos alertas sonoros e visuais na cabine dos equipa- 
mentos presentes na fase 2, houve também a intervenção da equipe de moni- 
toramento do sistema de gerenciamento de fadiga. A equipe de monitoramento 
foi responsável por monitorar, vinte e quatro horas por dia, o sistema de gestão 
de fadiga, acionando o procedimento estabelecido sempre que alarmes médios 
ou altos fossem recebidos por um operador. O procedimento preconiza que, 
caso o operador receba de um a quatro alarmes médios, deve-se paralisar o 
equipamento, realizar a ginástica laboral e/ou ingerir um lanche caso tenha se 
passado muito tempo desde a última ingestão de glicose. No caso de mais de 
quatro alarmes médios ou um alarme de nível alto, o operador deve parar o 
equipamento em local seguro e ser substituído naquele turno, como ilustrada 
na Figura 2. 
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Figura 2: Procedimento operacional para alarmes de fadiga. 


Vale ressaltar que na fase 3, foi designado um profissional habilitado para monitorar 
os resultados advindos do sistema e realizar as intervenções, quando necessário. 

Durante todo o período do experimento os dados foram coletados pelo sistema 
OptAlert?, possibilitando análises sobre o nível de risco da mina antes e depois da 
implantação do sistema e das intervenções feitas no caso de alertas de risco. 


DISCUSSÃO DE RESULTADOS 


Os dados foram analisados de maneira a verificar a efetividade da aplicação do 
sistema e das intervenções na redução do nível de risco dos operadores da operação. 
Inicialmente, analisou-se a pontuação média, medida em JDS, das fases do experimento. 
A Tabela 1 sumariza as principais estatísticas descritivas por fase, também ilustradas pelos 
boxplots da Figura 3, em que é possível observar a redução de cerca de 22% (1,70 para 
1,32) na pontuação média entre as fases 1 e 3. 


Estatísticas Descritivas Fase 1 Fase 2 Fase 3 
Média 1,70 1,47 1,32 
Erro padrão 0,00301 0,00279 0,00282 
Mediana 1,50 1,30 1,20 
Desvio padrão 1,28 1,13 1,03 
Variância da amostra 1,04 1,28 1,06 
Mínimo 0,00 0,00 0,00 
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Máximo 7,90 7,00 6,90 
n 180695 164136 132392 


Tabela 1: Estatísticas da pontuação média nas três fases do experimento. 


9,0 
8,0 
7,0 
6,0 
5,0 
4,0 
3,0 
2,0 
1,0 
0,0 


0000 00 00 0 


Pontuação Média (JDS) 








E Fasei BE Fase2 DD Fase3 


Figura 3: Boxplot da pontuação média por fase do experimento. 


A Figura 4 avalia a evolução da pontuação média por hora do dia. Não se pode 
deixar de notar a gradual evolução do nível de sonolência ao longo dos turnos, em que o 
turno noturno apresenta um nível de sonolência médio maior que os outros turnos e que o 
turno diurno é o que apresenta o menor nível. Os picos de sonolência acontecem no final do 
turno da tarde (entre 22h e Oh) e próximos ao final do turno noturno (entre 04h e 06h). Esse 
comportamento sazonal ao longo do dia se repetiu independente da fase do experimento. 
Além disso, há clara diferença na pontuação média entre as 3 fases do experimento, em 
que a intervenção direta da central de monitoramento, na fase 3, resultou em menores 
níveis de sonolência por hora. 
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Figura 4: Pontuação média (JDS) por hora do dia em cada fase do experimento. 


Outra análise se refere à quantidade de alarmes de risco gerados por hora do dia. 
Nesse sentido, a frequência média de alarmes por hora foi reduzida em cerca de 86% entre 
as fases 1 e 3, conforme Tabela 2. Na Figura 5 é possível observar que, além do evidente 
deslocamento de médias entre as fases do experimento, há redução na variabilidade dos 


resultados. 

Estatísticas Descritivas Fase 1 Fase 2 Fase 3 
Média 0,67 0,26 0,10 
Erro padrão 0,05077 0,02613 0,02328 
Mediana 0,70 0,25 0,07 
Desvio padrão 0,25 0,13 0,11 
Variância da amostra 0,06 0,02 0,01 
Mínimo 0,19 0,07 0,01 
Máximo 1,09 0,50 0,53 

n 24 24 24 


Tabela 2: Estatísticas da frequência de alarmes por hora nas três fases do experimento. 
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Figura 5: Frequência de alarmes por hora nas três fases do experimento. 


A redução da frequência de alarmes por hora se deu tanto para alarmes de risco 
de nível médio quanto para os de nível alto, bem como nos diferentes turnos e horas do 
dia, como pode ser comparado na Figura 6, Figura 7 e Figura 8. Nota-se que, mesmo nas 
diferentes fases do experimento, a frequência de alarmes por hora aumenta ao longo dos 
turnos, se mantendo mais concentrada no final do turno da tarde e do turno noturno, em 
conformidade com o resultado analisado de média de pontuação por hora. Tais resultados 
corroboram com o ciclo circadiano humano, uma vez que os operadores do experimento 


estavam mais cansados ao final dos turnos que se passam no período noturno. 
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Figura 6: Frequência de alertas por hora do dia - Fase 1. 
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Figura 7: Frequência de alertas por hora do dia - Fase 2. 
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Figura 8: Frequência de alertas por hora do dia - Fase 83. 


As análises mostram que a simples intervenção de alarmes sonoros e visuais na 
cabine dos equipamentos aumentaram o nível de alerta dos operadores e reduziram o nível 
de risco operacional, com uma redução de cerca de 61% na frequência de alarmes entre 
as fases 1 e 2. Já a intervenção direta de uma equipe responsável pelo monitoramento e 
gestão da fadiga dos operadores acarretou na uma redução ainda mais acentuada no perfil 
de risco da operação, visto que a frequência de alarmes entre as fases 1 e 3 foi reduzida 


em cerca de 85% e em 62% entre as fases 1 e 2. 


CONCLUSÃO 


Este estudo apresentou uma metodologia para verificar a efetividade de um sistema 
de gerenciamento de fadiga em aumentar o nível de alerta dos operadores de caminhão 
fora de estrada em uma mina a céu aberto e assim contribuir na redução de acidentes. 
Três etapas de teste, bem definidas, a saber: etapa “modo cego” (uso de óculos, porém 
sem alertas de sonolência na tela ou intervenções); etapa com alertas na tela e etapa 
com intervenção da central de monitoramento se fizeram necessárias para observar 
o aumento a evolução do nível de alerta dos operadores. Com base nos resultados e 


análises apresentados, concluiu-se que a implantação de um sistema de monitoramento e 
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gestão de fadiga capaz de determinar de maneira quantitativa a sonolência de operadores 
e estabelecer procedimentos e tratativas específicos é de grande valia e efetiva no que 
tange à gestão de riscos e pode ser útil na prevenção de acidentes em operações a céu 
aberto. Ademais, verificou-se que a intervenção através de uma central de monitoramento 
responsável pela gestão do sistema é de suma importância para o aumento do nível de 
alerta dos operadores e consequente redução do risco relacionado à fadiga. 


REFERÊNCIAS 


[1] Bauerle, T., Dugdale, Z., Poplin, G. (2018). Mineworker fatigue: A review of what we know and future 
decisions, 70 (3), pp. 33-40 


[2] Boivin, D., Boudreau, P. (2014). Impacts on shift work on sleep and circadian rhythms.Payhologir 
Biologie, vol 62 (5), pp.292-301. 


[3] Everitt, B. An introduction to latent variable models. Chapman and Hall; London: 1984.107 pp. 


[4] Shen J, Barbera J, Shapiro C. Distinguishing sleepiness and fatigue: Focus on definition and 
measurement. Sleep Medicine Reviews. 2006;10:63-76. [PubMed] 


[5] Catterpillar Global Mining, Viewpoint. Operator Fatigue Detection Technology Review Caterpillar, 
Peoria, IL (2008). 


[6] Pelders e Nelson (2019) Contributors to Fatigue of MineWorkers in the South African Gold and 
platinum sector. 10:188-195. 


[7] Johns, M., Chapman, R., Crowley, K., Tucker, A. (2008) A new method for assessing the risks of 
drowsiness while driving. Somnologie - Schlafforschung und Schlafmedizin. 12, 66-74(2008) 


[8] Aidman,E., Chadunow,C., Johnson, K., Reece, J. Real-time driver drowsiness feedback improves 
driver alertness and self-reported driving performance. Accident Analysis & Prevention. 81: 8-13. 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 10 EE 


CAPÍTULO 11 


CHARACTERISTICS OF PB* DOPED CSI MATRIX 
UNDER GAMMA AND NEUTRON EXCITATIONS 


Data de aceite: 01/11/2021 
Data de submissão: 08/10/2021 


Maria da Conceição Costa Pereira 

Instituto de Pesquisas Energéticas e Nucleares 
(IPEN / CNEN - SP) 

São Paulo, SP — Brazil 
http://lattes.cnpg.br/0977986453736808 


Tufic Madi Filho 

Instituto de Pesquisas Energéticas e Nucleares 
(IPEN / CNEN - SP) 

São Paulo, SP — Brazil 

Universidade Paulista — UNIP 

Campus Tatuapé 

São Paulo - Brazil 
http:/lattes.cnpqg.br/3251192790149021 


José Roberto Berretta 

Instituto de Pesquisas Energéticas e Nucleares 
(IPEN / CNEN - SP) 

São Paulo, SP — Brazil 
http:/lattes.cnpqg.br/4444753560161862 


Lucas Faustino Tomaz 

Instituto de Pesquisas Energéticas e Nucleares 
(IPEN / CNEN - SP) 

São Paulo, SP — Brazil 
http:/lattes.cnpq.br19015111335731083 


ABSTRACT: In recent years, there has been an 
increasing interest in finding new fast scintillating 
material or improve the characteristics of known 
scintillatorsfor the demand of high energy physics, 
industrial and nuclear medical applications. 
Divalent lead ions Pb2* built in some crystal 
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structures are efficient emission centers and 
their applications in scintillators wereand aresstill 
the reason of an intensive study of emission 
properties of different compounds containing 
these ions. In this context, the crystals of Pb2* 
doped Csl matrix were grown by the vertical 
Bridgman techniqueandsubjected to annealing in 
vacuum of 10º mbar and constanttemperatureof 
350ºC, for 24 hours, and then they were 
employed. To evaluate the response of the 
Csl:Pb scintillator crystal to gamma radiation, 
radioactive sources of !ºCs (662 keV), “Co 
(1173 keV and 1333 keV), 2Na (511 keV and 
1275 keV) and '“Ba (355 keV) wereused. The 
operating voltage of the photomultiplier was 
2700 V for the detection of gamma rays and 
the accumulation time in the counting process 
was 600 s. The scintillator response to neutron 
radiation from a radioactive source of AmBe with 
energy range of 1 to 12 MeV was available. The 
activity of the AmBe source was 1 Ci Am. The 
emission ratewas 2.2x 10º neutrons/ second. The 
operating voltage of the photomultiplier tube was 
1300 V. The accumulation time in the counting 
process was 600 s. With the results obtained, it 
may be observed that the crystals are sensitive to 
these radiations. 

KEYWORDS:Scintillator, crystal growth, gamma 
radiation, neutron radiation. 


CARACTERÍSTICAS DA MATRIZ Cs 
DOPADA COM Pb?* SOB EXCITAÇÕES 
GAMA E NÉUTRONS 


RESUMO: Nos últimos anos, tem havido 
um interesse crescente em encontrar novos 


Capítulo 11 


materiais de cintilação rápida ou melhorar as características de cintiladores conhecidos para 
atender a demanda de física de alta energia, aplicações médicas, industriais e nucleares. 
Os íons divalentes de chumbo Pb? inseridos em algumas estruturas cristalinas são centros 
emissores eficientes e suas aplicações em cintiladores foram e ainda são motivo de um 
estudo intensivo das propriedades de emissão de diferentes compostos contendo esses íons. 
Nesse contexto, os cristais de matriz Csl dopada com Pb2+ foram crescidos pela técnica de 
Bridgman vertical e submetidos ao tratamento térmico a vácuo de 10º mbar e temperatura 
constante de 350 º C, por 24 horas. Para avaliar a resposta do cristal cintilador Csl:Pb à 
radiação gama, foram utilizadas fontes radioativas de 'º/Cs (662 keV), “Co (1173 keV e 
1333 keV), 2Na (511 keV e 1275 keV) e “Ba (355 keV). A tensão de operação do tubo 
fotomultiplicador foi de 2700 V para a detecção dos raios gama e o tempo de acumulação no 
processo de contagem foi de 600 s. A resposta do cintilador à radiação de nêutrons de uma 
fonte radioativa de AmBe com intervalo de energia de 1 a 12 MeV foi estudada. A atividade 
da fonte AmBe foi de 1 Ci Am. A taxa de emissão foi de 2,2 x 10º nêutrons/segundo. A tensão 
de operação do tubo fotomultiplicador foi de 1300 V. O tempo de acumulação no processo 
de contagem foi de 600 s. Com os resultados obtidos, pode-se observar que os cristais são 
sensíveis a essas radiações. 

PALAVRAS-CHAVE: Cintilador, crescimento de cristal, radiação gama, radiação de nêutrons. 


11 INTRODUCTION 


The scintillation method is widely used for detection of various types of ionizing 
radiation. A scintillator consists of material able to convert energy of absorber ionizing 
radiation into a flash light. The flash light, usually visible or ultraviolet, emitted by a 
scintillator, is called scintillation. Any scintillator should be transparent to the wavelength of 
its own emission. The use of scintillation to detect radiation is a century old. The discovery 
of scintillator materials may be divided into three periods. The first period included the 
earliest scintillators: CaWO, first used in the year following Rontgen's discovery for X-rays; 
uranyl salts used by Becquerel to discover radioactivity and by Rutherford to study alpha 
particle scattering. The second period began with Hofstadter's development of the thallium 
activated Nal. The third period, the past two decades, has witnessed a veritable renaissance 
in research and development of scintillator materials, prompted in large part by the need 
for scintillators for precision calorimetry in high energy physics and for high light output 
scintillators for medical imaging, geophysical exploration, tomography, astrophysics and 
numerous other scientific and industrial applications. [7,16] 

There is still considerable interest in finding inorganic materials with perfect 
scintillation properties or improving the traditional scintillators, including álcali halides. In 
order to enhance the characteristics of a known scintillator or to obtain a new crystal, it 
is necessary to solve many problems, such as crystalline structure, chemical binding or 
impurity concentration, which is more advantageous for high light output of a crystal. [17] 

Scintillators based on cesium iodide are the leading ones among the materials 
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available for radiation detectors. Cesium iodide based scintillators have high luminescence 
yield accompanied by high energy resolution and radiation stability. Cesium iodide has 
elements of ionic radius, namely: the cesium ion (Cs*) radius is 1.652 Á, electronegativity 
0.7 eV and the iodine ion (I') radius is 2.256 Á; electronegativity 2.5 eV the equilibrium value 
lattice separation is 4.512 Á. This facilitates the incorporation of a large number of cationic 
impurities. Various impurities have been included into alkali halides and their effects on 
electrical, magnetic and optical properties have been investigated. 

The Csl matrix is poorly hygroscopic, its constituents have a high atomic number, 
are easy to handle and have relatively low cost.[15] Babin et al. [2] mention the interest in 
Csl:Pb crystals as possible materials for scintillators, since their luminescence is practically 
not studied. Some data reported by different authors are quite contradictory. This indicates 
a complexity of this system and a strong dependence on the optical characteristics of the 
crystals that may be related to the method used for growth. According to Keszthelyi et al. 
[6], the scintillation characteristics of the crystals are also influenced by the growth method 
used. 

Divalent lead ions Pb? in some crystal structures are efficient emission centers and 
their application in scintillators was and is still the reason of an intensive study of emission 
properties of different compounds containing these ions. The study of pure and doped 
cesium iodide has been the theme of search in our laboratory. In continuation, a systematic 
investigation has been carried out in Pb?* doped cesium iodide crystals under gamma and 
neutron excitations and the results are discussed in this paper. 


21 MATERIALS AND METHODS 


Lead doped Csl crystals were grown using the vertical Bridgman technique. 
[3,12,13,9,14] The grown crystals were subjected to annealing. In this procedure, vacuum 
of 10º mbar and constant temperature of 350º C, for 24 hours, were employed. The crystals 
were cut into cylindrical form and their faces were polished with ethylene glycol p.A. 
(C,H,O,). The side surfaces were left unpolished to enhance internal reflection. 

The luminescence decay time of each crystal was determined under excitation of 
662 keV gamma rays from a !º Cs source. The scintillation photons were detected using a 
bi-alkaline photomultiplier tube (model RCA 8575, 21 pins). The conformation of the signal 
from the photomultiplier tube was obtained using the RC (resistor-capacitor) circuit having a 
time constant of approximately 5 ns and fed directly to the oscilloscope (Philips PM3295A), 
resulting in the measurement of the decay of the scintillation pulses. In Fig. 1 the electronic 


system used for luminescence decay time measurements is schematically show. 
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Figure 1: Diagram of the electronic system used for luminescence decay time measurements. 


In the gamma ray response study, the crystals were cut with dimensions of 20 mm 
diameter and 20 mm high, polished and coupled direcily to the bi-alkaline photomultiplier 
tube (model RCA 8575, 21-pin) using silicone grease (Dow Corning) of 0.5 McStokes 
viscosity as optical interface. This ensures a uniform refractive index across the contact 
surface between crystal and photomultiplier tube. The sides of the crystal which were not 
in contact with the photosensor were covered with several layers of polytetrafluoroethylene 
(PTFE) tape to assure good reflection of light. 

The radioactive sources were positioned in the center of the upper face of the crystal. 
Conventinal ORTEC eletronic modules constitute the electronics for the treatment of the 
signals from the photomultiplier tube. Gamma radiation sources with energies in the range 
of 335 keV to 1333 keV were used. The operating voltage of the photomultiplier was 2700 V 
and the accumulation time in the counting process was 600 s. In Fig. 2, the diagram of the 
measuring system is showed. 
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Figure 2: Schematic representation of the electronics associated with the Csl:Pb scintillator crystal 
coupled to the photomultiplier tube. 
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The AmBe source with energy range of 1 to 12 MeV was used to study the neutron 


radiation response. The activity of the AmBe source was 1Ci Am. The emission rate was 2.6 


x 10º neutrons/second. The operating voltage of the photomultiplier tube was 1700 V; the 


accumulation time in the counting process was 1800 s. The scintillator used was cut with 


dimensions of 20 mm diameter and 10 mm high. 


The detection efficiency of the scintillator crystal was measured with the AmBe source 


positioned at a distance of 100 mm from the photomultiplier tube. A foil of cadmium (Cd) 


was placed around the detector and the photomultiplier tube in order to avoid the scattered 


neutron contribution, as showed in Fig.3. Paraffin shielding to protect the worker, was used. 
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Figure 3: Exploded schematic representation of the array used for measurements. 


31 RESULTS AND DISCUSSION 
In Fig. 5, the luminescence decay curves of the Csl:Pb crystals are presented. For 


comparison purposes, the luminescence decay constant was determined for the pure Csl 


crystal, as shown in Fig. 4. 
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Figure 4: Luminescence decay time spectrum of pure Csl crystal. 
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Figure 5: Luminescence decay time spectra of Csl:Pb crystals at various concentrations: (a) 5x10“ M, 
(b) 10% M, (c) 1,5x10% M and (d) 102 M. 


Fig. 5 and Fig.4 show the luminescence decay curves for the Csl:Pb and pure Csl 
crystals. As it may be observed in these curves, there were no significant differences for 
the various concentrations of lead in the crystals evaluated. The incorporation of the Pb 
into the Csl matrix, did not, practically, change the luminescence decay time compared to 
that obtained for pure Csl, whose values are between 13 ns and 17 ns. These results are 
very promising, i. g. to obtain crystals having shorter decay time comparable to pure Csl 
scintillator crystal. The crystal of pure Csl produced in this work, for comparison purposes 
showed a Iluminescence decay constant of 12 ns in accordance with the literature [5,1,10]. 

The luminescence decay time parameter is one aspect of crystal quality, especially 
for applications requiring high radiation rates or high energy radiations. The luminescence 
decay time results found in this work are encouraging since, with the introduction of dopant 
ions such as TI* or Na* in the Csl matrix, it results in increased luminescence decay time 
for 1 s and 630 ns respectively. The Csl:TI crystal is one of the most efficient scintillators in 
gamma ray detection per unit volume, but the long decay time, in the order of 1s, is one ofits 
disadvantages, increasing the probability of occurrence of overlapping signals, when more 
intense radioactive sources are used. The decay time in the order of 1s, practically makes it 


unfeasible rapid temporal studies, such as those of the decay structure performed in nuclear 
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physics. The Csl:Na crystal presents good light production, but the luminescence decay 
time is still high for specific applications and its main disadvantage is the hygroscopy [7,16]. 

Fig. 6 to 9 show the gamma radiation spectroscopy results of !*”Cs (662 keV), “Co 
(1173 keV and 1333 keV), 2Na (511 keV and 1275 keV) and '“Ba (355 keV) obtained with 
the Csl:Pb crystals. It may be observed that the crystals are sensitive to these radiations; 


however, the photopeak of these radionuclides have not been identified. 
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Figure 6: Spectra obtained for the radiation of "Cs, “Co, Na and '““Ba sources with Csl:Pb 5x10“M 
crystal. 
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Figure 7: Spectra obtained for the radiation of "Cs, “Co, Na and '“Ba sources with Csl:Pb 10ºM 
crystal. 
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Figure 8: Spectra obtained for the radiation of "Cs, “Co, 2Na and '““Ba sources with Csl:Pb 1.5x10º%M 
crystal. 
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Figure 9: Spectra obtained for the radiation of "Cs, “Co, “Na and '““Ba sources with Csl:Pb 102M 
crystal. 


From the curves obtained for the Csl:Pb crystals, using gamma radiation from the 
sources of 3 Cs, “Co, 2Na and '“Ba (Figs. 6-9) it was not possible to extract the resolution 
of these detectors. The low resolution capacity of the crystals of Csl:Pb may be attributed to 
the self- absorption of the Iluminescence in the crystals themselves, due to the presence of 
Pb in their structures. However, for applications where the energy resolution is not a relevant 
item, the Csl:Pb crystals may be used when a radiation detector with good mechanical 
resistance, relatively low cost of production and different sizes and geometries is desired. 

In electromagnetic calorimeters, the need for large-scale high-energy particle 
detectors has led to international collaboration to develop, optimize and produce scintillators 
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for this purpose. The Compact Muon Solenoid (CMS) experiment in the Large Hadron 
Collider (LHC) accelerator at CERN uses approximately 80,000 PbWO, (lead tungstate) 
crystals, totalling approximately 100 tonnes [11]. Another project, ALICE, also in the 
accelerator (LHC) will use 18,000 PLWO, crystals, about 13 tons [4]. The crystal scintillator 
PbWO, presents light output of 250 photons / MeV and decay time of 5 to 15 ns [4,8]. 

Fig. 10 shows of measurements performed using the Csl:Pb crystals with 


concentrations 102 M, 10ºM and 5x104 M and the neutron source. 
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Figure 10: Results of the measurements using Csl:Pb crystal and neutron source. 


Csl:Pb crystals, when excited with radioactive source of AmBe with energy range 
of 1 to 12 MeV, showed sensitivity to fast neutron, and the highest counting efficiency was 


obtained in the concentration of 102 M. 


41 CONCLUSION 


The Csl: Pb crystals presented short decay time values, demonstrating the viability 
of their use for measurements of high dose and high energy rates. These results have 
shown to be very promising to obtain crystals having short decay time, when compared to 
the pure Csl scintillator crystal. 

The crystals of Csl:Pb showed sensitivity to gamma radiation detecting the energies 
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of the radioactive sources used, although it was not possible to identify the corresponding 
photopeak. 

The addition of the lead ions (Pb*) to the Csl matrix resulted in crystals with promising 
results, when excited with neutron radiation. The crystal showed to be sensitive to fast and 
thermal neutron. 

In the dopant concentration range (Pb) studied, the best neutron detection efficiency 
was obtained for the molar fraction 107. 

The results showed the feasibility of using the crystals developed in this work for 
application in high energy physics. In our work environment, they may be applied as area 
of a neutron detector and used in the personal dosimeter for the workers of the IEA-R1 and 
IPEN/MB-01 nuclear research reactors. 
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RESUMO: O método de Taguchi/projeto robusto 
foi desenvolvido pelo engenheiro e estatístico 
Genichi Taguchi (1924 - 2012) sendo denominado 
por ele de controle off-line. Este método, segundo 
o estatístico, foi idealizado como uma nova 
abordagem da qualidade com ênfase no produto 
e no processo. A qualidade é medida pelo desvio 
que uma característica funcional apresenta em 
relação ao valor base esperado da mesma. Os 
fatores responsáveis por esses desvios são 
chamados de “ruídos”(temperatura, poeira, 
vibração, umidade, etc) que acabam afetando a 
qualidade do produto. A função perda, proposta 
por Taguchi, modela estas perdas causadas 
pelos ruídos. Afim de ter um produto de qualidade 
devemos buscar otimizar a função perda do 
mesmo utilizando técnicas estatísticas para 
poder ter um resultado com a menor perda por 
ruídos possível. Utilizando o método de Taguchi 
buscamos obter um produto robusto, ou seja, 
um produto que não altere suas condições de 
funcionamento ou trabalho por uma grande gama 
de variáveis que venham a produzir ruídos. Além 
da otimização da função perda o método engloba 
um conjunto de técnicas propostas para atingir o 
objetivo, por exemplo, os experimentos fatoriais 
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com matrizes ortogonais. Por ser um método 
muito prático, sendo esta uma das preocupações 
no momento de sua idealização, ficou bem 
reconhecido por empresas que começaram a 
implantá-lo em suas diretrizes. O projeto feito 
neste artigo foi desenvolvido utilizando o método 
de Taguchi, especificamente os experimentos 
fatoriais com matrizes ortogonais. O objetivo 
do trabalho é procurar a melhor combinação de 
ingredientes buscando o melhor sabor possível 
para a bebida junto aos clientes. 
PALAVRAS-CHAVE:  Taguchi, 
Matriz ortogonal, Produtividade. 


Otimização, 


OPTIMIZING THE FLAVOR OF A 
BEVERAGE USING THE TAGUCHI 
METHOD 


ABSTRACT: The Taguchi method / robust 
design was developed by the engineer and 
statistician Genichi Taguchi (1924 - 2012) being 
called by him of control offline. This method, 
according to the statistician, was conceived as 
a new approach to quality with an emphasis on 
product and process. The quality is measured 
by the deviation that a functional characteristic 
presents in relation to the expected base value 
of the same. The factors responsible for these 
deviations are called “noises” (temperature, dust, 
vibration, humidity, etc.) that end up affecting the 
quality of the product. The loss function, proposed 
by Taguchi, models these losses caused by 
noise. In order to have a quality product we must 
seek to optimize the loss function of the same 
using statistical techniques to be able to have a 
result with the least noise loss possible. Using 
the Taguchi method we seek to obtain a robust 
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product, that is, a product that does not alter its working or working conditions by a wide range 
of variables that will produce noises. In addition to the optimization of the loss function, the 
method encompasses a set of techniques proposed to achieve the objective, for example, 
the orthogonal matrix factorial experiments. Being a very practical method, being one of the 
preoccupations at the moment of its idealization, was well recognized by companies that 
began to implant it in its directives. The design made in this article was developed using the 
Taguchi method, specifically the orthogonal matrix factorial experiments. The objective of the 
work is to look for the best combination of ingredients seeking the best possible flavor for the 
beverage with the customers. 

KEYWORDS: Taguchi, Optimization, Orthogonal Matrix, Productivity. 


11 INTRODUÇÃO 


Atualmente, a indústria está procurando maneiras de reduzir custos, despesas 
e maximizar a produtividade, consequentemente aumentando os lucros. A aplicação de 
métodos estatísticos nos processos de produção e prestação de serviços tem se tornado 
um fator importante e diferenciador entre as organizações. 

O estudo estatístico utilizando a variância (ANOVA) de determinado fenômeno 
permite o controle da variabilidade, ou seja, um padrão de comportamento pode ser 
definido, diminuindo os ruídos, que são os principais responsáveis pela falta de precisão e 
acurácia nos processos dentro de uma indústria. 

O planejamento é embasado em análises de dados, que em sua grande maioria, são 
realizadas de maneira estatística. 

O método de Taguchi pode ser usado para determinar a combinação ótima de 
fatores e interações que influenciam o comportamento das variáveis de resposta em um 
determinado processo. Em geral, esse método requer menos esforço de amostragem e 
menores custos de teste, sem afetar significativamente as conclusões tiradas (ROSA et 
al., 2009). 

Este artigo trata da aplicação do método Taguchi como ferramenta estatística para a 
análise e investigação dos principais fatores e interações que influenciam no melhor sabor 


de suco dentre as combinações propostas pelo estudo para uma pesquisa de mercado. 


21 FUNDAMENTAÇÃO TEÓRICA 


A pesquisa de mercado é uma etapa que auxilia a indústria no que tange os anseios 
dos clientes. Deve-se saber o que vender, para isso, saber o que o cliente quer comprar 
se torna um quesito chave. Para responder a esta pergunta nada mais prático que ouvir do 
próprio interessado no produto, ou seja, o consumidor, sobre o que necessita, o que deseja 
adquirir. 

Lima et al. (2011) explicam que planejar experimentos é a forma mais adequada 


para estudar a complexidade dos vários fatores do processo e suas interações, a fim 
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de aumentar a probabilidade de resolução de problemas por meio de análise estatística 
(ANOVA). Esse método, que se acredita ter melhorado a qualidade e a produtividade, tem 
sido cada vez mais aplicado nos últimos anos na indústria brasileira (SILVA e SILVA, 2008). 

A obtenção de dados, nesse caso, geralmente é realizada por meio de questionários 
ou entrevistas com os consumidores em potencial. Para a análise desses dados, a 
ferramenta taguchi/projeto robusto, que foi idealizada pelo estudioso japonês Dr. Genichi 
Taguchi, se torna um artifício muito prático, didático e de fácil interpretação, tendo em 
vista que se baseia na melhora das características dos processos de um produto através 
do controle e ajuste dos chamados fatores controláveis. Estes minimizam a variação 
do produto final em comparação com o resultado pretendido que é estabelecido no 
começo da análise. Ajustando os fatores importantes para o nível ótimo, a fabricação dos 
produtos pode ser realizada de forma que se adeque a uma maior gama de situações, tal 
característica é observada nos ditos produtos robustos. Os fatores incontroláveis, como 
condições ambientais, variação dimensional, dentre vários outros são os responsáveis pela 
variabilidade do processo, portanto, quanto mais robusto um produto se apresenta, menos 
se observa variação no seu objetivo final. 

Nos anos 1950 e 1960 o Dr. Genichi Taguchi desenvolveu um design robusto e 
aplicou essa abordagem em muitas áreas (SOUZA et al., 2011). 

Primeiramente devemos passar pelas etapas básicas para a aplicação da 
metodologia Taguchi. Esta se divide em quatro: 


l. Identificação dos fatores 


Nesta etapa será realizada a identificação dos fatores (ruído e fatores do ambiente 
e processo de fabricação) e os parâmetros dos produtos. 

Para estes fatores são analisados as possíveis influências e iterações com os outros. 

Deve-se tomar cuidado para não deixar de fora algum parâmetro muito importante 
e significativo na hora da escolha dentre os vários parâmetros que o problema venha a ter. 
Caso isso aconteça ocorrerá uma distorção da função perda. 


Il. Planejamento e condução dos experimentos 


Este é o momento do planejamento de coleta dos dados experimentais. A partir 
destes dados podemos construir a relação sinal/ruído. 

Para realizar o planejamento deve-se iniciar pela escolha do tipo de planejamento, 
ou seja, pela escolha da matriz ortogonal que melhor se aplica ao problema. 

A escolha das matrizes dependem principalmente do número de fatores e da 
quantidade de corridas (ou seja, de casos de experimentos) que poderiam ser realizados 
conforme a disponibilidade de tempo e custo. 

Em seguida são especificados valores para os diferentes níveis dos parâmetros. 

Com estes dados basta escolher de forma aleatória as corridas e programar a 
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realização dos ensaios. 


II. Análise dos níveis Ótimos dos parâmetros 


É feita a otimização dos parâmetros do produto a partir da relação com o sinal/ruído. 

Isto significa obter um modelo estatístico desta relação com os dados coletados no 
experimento e aplicar, neste modelo, técnicas de otimização para encontrar os valores dos 
parâmetros Ótimos dos produtos. 

Ao final desta etapa tem-se um conjunto de valores de parâmetros (ou características) 
do produto que tornam seu desempenho robusto e estável em relação as características 


ambientais e às variações do processo. 


IV. Validação dos resultados 


Como os níveis ótimos dos parâmetros obtidos anteriormente são fruto de um modelo 
estatístico, e, portanto, uma aproximação da realidade, deve-se realizar uma etapa de 
validação dos resultados encontrados, ou seja, verificação dos níveis ótimos especificados 
para os parâmetros. 

Isto é feito conduzindo um experimento com um protótipo cujos parâmetros são 
ajustados conforme os valores ótimos obtidos na fase anterior. 

Os resultados deste experimento devem coincidir com àqueles encontrados por 
meio do modelo, dentro, é claro, da devida margem de segurança. 

Caso isto ocorra significa que o modelo obtido é confiável e, portanto, pode-se 
aprovar estes parâmetros como especificações para o projeto. 

O arranjo ortogonal é uma matriz com os dados, parâmetros e níveis dos 
experimentos realizados. Existem diferentes tipos de dimensões para as tabelas e cabe 
ao projetista analisar e utilizar dentro do seu número de parâmetro e níveis do problema 
investigado a dimensão para o seu caso, explica Taguchi (1990). Para a construção de um 
arranjo ortogonal é preciso que cada fator seja ortogonal ao outro, ou seja, tenha os seus 
níveis aparecendo o mesmo número de vezes. 

Taguchi explica que o arranjo ortogonal está diretamente ligado ao grau de liberdade, 
que é calculado através da subtração do número de fatores do processo estudado menos 
um. Através da contagem do grau de liberdade é concebido o arranjo ortogonal (ROSS, 
1995). 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 12 RE 


Engenharias: Da genialidade à profissão e seu desenvolvimento 2 Capítulo 12 








Es | g En 
: no 7 
E y ! E É p 
- E | 5s E 1 





Tabela 1 - Matriz ortogonal para o método experimental de Taguchi. 


Fonte: adaptado Montgomery (2004) 


A matriz varia conforme o número de fatores e níveis que serão utilizados. Tais 
fatores são determinados após consulta na literatura ou por protocolos já existentes 
(HUANG; HUNG; YANG, 2016). 

Em uma análise experimental Taguchi sugere a análise da resposta média de cada 
linha no arranjo mais interno e também a análise da variação utilizando uma escolha 
apropriada de relação Sinal/Ruído (S/N). Essas relações S/N derivam da função-perda 
quadrática, tendo três tipos padrão para análise, segundo Souza et. al. (2011): 


a. Nominal é melhor 


Idealmente, para este tipo de função de perda, a resposta deve ter um valor de 
trabalho igual ao valor base especificado. Por exemplo, na fabricação de peças de máquinas, 
o cumprimento dimensional estrito é exigido de acordo com as regras de montagem. Caso 
contrário, as peças não funcionarão juntas, pois geralmente são feitas de um material muito 
duro e durável. Portanto, o valor nominal necessário para o modelo é ótimo. 


A equação que rege este modelo é uma função quadrática descrita logo a seguir: 


L(y)=k(y-m? 
Onde: 


- Léa perda associada a uma partícula característica de performance, no caso 
Oy, 
*- mé o valor objetivo; 
*  Kéo parâmetro de perdas; 
* yéaresposita de cada experimento. 
A relação sinal/ruído para esta situação é designada por: 


) 


+ 
12 





SN, = 10log 





Onde: 
* Yéa média de cada experimento; 


*- Séo desvio padrão das respostas para determinada combinação de níveis de 
parâmetros. 


b. Menor é melhor 


Como exemplo para este caso podemos citar o consumo de combustível de um 


carro ou o tempo que o computador demora pra iniciar. É regida peça equação: 


A relação sinal/ruído para esta situação é definida pela equação: 


SN, =-10log | E pa ] 


Hal 
Onde: 
* Yéaresposta de cada experimento; 


*. né o número de respostas de cada experimento. 


c. Maior melhor 
Neste caso o valor ideal seria infinito e neste ponto a perda seria zero. Porém, isso 
não é possível, portanto, através deste método podemos otimizar para um maior valor 


positivo possível. A equação que rege este método é: 





Onde: 
*  Kkéo parâmetro de perdas; 
* yéaresposta de cada experimento. 


A relação Sinal/Ruído para esta situação é: 





SN, =—10l08] — 5 — 


31 METODOLOGIA 


Este projeto tem como objetivo através da combinação de alguns ingredientes 
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determinar qual o melhor sabor para uma bebida, levando em consideração dosagens e 


tipos de frutas. Para isto estabelecemos os parâmetros e níveis a seguir: 


Parâmetro Descrição | Nível! | Nível II 


A Suco de laranja 125 ml 


Abacaxi 20 g 
Maçã 120 9 
Água 40 ml 

60 g 





MIO O W 


Tabela 2 - Parâmetros e níveis do experimento 


Fonte: Dados colhidos mediante pesquisa do autor. 


A matriz ortogonal escolhida para este projeto foi a L8 com 7 observações para cada 


experimento: 


OI" SO AB ÇGOIN ma 
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Tabela 3 - Tabela L8 dos experimentos 


Fonte: Dados colhidos mediante pesquisa do autor. 


A tabela a seguir foi obtida pelos experimentos feitos. Para este procedimento 
foram analisadas as respostas da seguinte forma, desagradável (1), regular(3), bom(5) 
e foi calculado a média dos experimentos. Cada combinação de sabor foi provada por 
duas pessoas. Para manter o critério de avaliação constante, cada um provou as oito 


combinações possíveis. Para o primeiro indivíduo: 
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1 5 1 1 1 5 5 1 2,11 
2 5 3 3 1 3 1 3 2,71 
3 3 5 1 3 1 5 5 3,29 
4 5 3 3 1 3 1 1 2,43 
5 5 5 1 3 3 1 5 3,29 
6 5 5 5 5 5 3 3 4,43 
7 5 5 3 3 1 3 3 3,29 
Tabela 4 - Experimentos para o primeiro indivíduo 
Fonte: Dados colhidos mediante pesquisa do autor. 
Para o segundo indivíduo: 

LO A B| c| D| E Fl c|(ivy 
1 5 3 1 1 5 5 1 3 
2 5 1 3 1 3 1 3 2,42 
3 3 3 1 3 1 5 5 3 
4 5 5 3 1 3 1 1 2,11 
5 5 5 1 & 3 1 5 3,29 
6 5 5 5 5 5 3 3 4,43 
7 5 5 3 3 3 3 3 3,57 





Tabela 5 - Experimentos para o segundo indivíduo 


Fonte: Dados colhidos mediante pesquisa do autor. 


Portanto, a tabela com os experimentos e a relação sinal ruído (maior melhor) é: 


Nº S/N(dB) 
9,07 
8,14 
9,92 
8,16 
10,34 
12,93 
10,68 
8,16 


Oo" o qr GO INa 





Tabela 6 - Experimentos e resultados 


Fonte: Dados colhidos mediante pesquisa do autor. 
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A influência de cada parâmetro é dada pela média das respostas para o nível auto 
do parâmetro menos a média das respostas para o nível baixo do parâmetro. Portanto, 
calculando a influencia de cada parâmetro: 


Parâmetro A: 
E(A) = [(3,29 + 4,43 + 3,43 + 2,57) - (2,865 + 2,56 + 3,14 + 2,57)]/4 
E (A) = 0,65 


Parâmetro B: 
E(B) = [(2,85 + 2,56 + 3,29 + 4,43) - (3,14 + 2,57 + 3,43 + 2,57)]/4 
E(B) = 0,355 


Parâmetro C: 
E(C) = [(3,14+ 2,57 + 3,29 + 4,43) - (2,85 + 2,56 + 3,43 + 2,57)]/4 
E(C) = 0,505 


Parâmetro D: 
E(D) = [(2,85 + 3,14 + 3,29 + 3,43) - (2,56 + 2,57 + 4,43 + 2,57)]/4 
E(D) = 0,145 


Parâmetro E: 
E(E) = [(2,56 + 2,57 + 3,29 + 3,43) - (2,865 + 3,14 + 4,43 + 2,57)]/4 
E(E) = -0,285 


Podemos então analisar a influência de cada parâmetro e destacar que três deles 


tem um impacto grande positivamente no resultado, calculando a interação entre eles: 


Média 
2,85 
2,06 
3,14 
2,51 
3,29 
4,43 
3,43 
2,91 





Tabela 7 - Influência dos parâmetros 


Fonte: Dados colhidos mediante pesquisa do autor. 
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Interação AB: 
E(AB) = [(2,95 + 2,57 +3,29+ 2,57) - (2,56 + 3,14 + 4,43 + 3,43)]/4 


E(AB) = -0,57 
Interação AC: 
E(AC) = [(2,85 + 3,14 + 4,43 + 2,57) - (2,56 + 2,57 + 3,29 + 3,43)]/4 
E(AC) = 0,285 


Interação BC: 
E(BC) = [(2,85 + 2,56 + 3,43 + 2,57) - (3,14 + 2,57 + 3,29 + 4,43)]/4 
E(BC) = -0,505 


Interação (ABC): 
E(ABC) = [(2,56 + 3,14 + 3,29 + 2,57) - (2,85 + 2,57 + 4,43 + 3,43)]/4 
E(ABC) = -0,43 


41 RESULTADOS E DISCUSSÃO 


Como estamos buscando maximizar a qualidade do sabor do suco vamos analisar 
quais parâmetros mais contribuíram positivamente para um bom sabor da bebida. 
|. O parâmetro “A”, suco de laranja, contribui com 0,65 unidades positivamente. 


|. O parâmetro “B”, abacaxi, contribui com 0,355 unidades positivamente. 
II. O parâmetro “C”, maçã, contribui com 0,505 unidades positivamente. 
IV. O parâmetro “D”, água, contribui com 0,145 unidades positivamente. 
V. O parâmetro “E”, açúcar, contribui com 0,285 unidades negativamente. 


VI. A interação “AB”, suco de laranja e abacaxi, contribuiu com 0,57 unidades 
negativamente. 


VII. Ainteração “AC”, suco de laranja e maçã, contribuiu com 0,285 positivamente. 


VII. A interação “BC”, abacaxi e maçã, contribuiu com 0,505 unidades 
negativamente. 


IX. A interação “ABC”, suco de laranja, abacaxi e maçã, contribuiu com 
0,43 unidades negativamente. experimento 6 foi o que mais contribuiu 
positivamente para o sabor. Ele contribuiu em média com 4,43 unidades. 
Como exposto, os parâmetros “A”, “B” e “C” foram os que mais contribuíram 
positivamente para a bebida com uma porcentagem de 77,83% da contribuição total dos 
parâmetros. O açúcar, parâmetro “E”, contribuiu negativamente, ou seja, está diminuindo o 
grau de excelência da bebida. 


Podemos ainda definir o nível mais importante para o sabor de cada parâmetro: 
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Parâmetros/níveis A B C E 
Nível nz 13,13 11,41 12,71 12,99 
Nível II 13,72 11,71 13,43 12,18 11,85 











Tabela 8 - Contribuição de cada nível 


Fonte: Dados colhidos mediante pesquisa do autor. 


Sabendo a contribuição de cada parâmetro e de cada nível, dentre os experimentos 


podemos analisar qual combinação de sabor tem um maior impacto positivo no palato dos 


clientes, ou seja, aquela combinação que está próxima do valor bom (5). O experimento 6 


é O que mais se adequa a exigência, abaixo podemos ver a sua combinação que ficou com 
4,43 unidades positivas, ou seja, 88,6% do valor julgado bom (5). A pior combinação foi a 


do experimento 2 com uma contribuição de 2,56 unidades positivas, ou seja, ficou abaixo 


da condição regular (3) e com uma porcentagem em relação ao valor considerado bom (5) 


de 51,2%. 


Parâmetro 
A - Nível 2 
B - Nível 1 
C - Nível 2 
D - Nível 1 
E - Nível 1 


Tabela 9 - Melhor combinação de parâmetros 


Descrição 


125 ml de suco de laranja 


40 g de abacaxi 


120 g de maçã 
40 ml de água 


20 g de açúcar 


Fonte: Dados colhidos mediante pesquisa do autor. 


Podemos analisar a resposta e comparar com o gráfico da contribuição dos 


parâmetros para as médias: 
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Gráfico de efeitos principais para as médias 
Suco de laranja Abacaxi Maça Água Açúcar 
34 
33 


3,2 


Médias 


3,1 
3,0 
Za 
28 


2 


Figura 1 - Efeito dos parâmetros nas médias. 


Fonte: Dados colhidos mediante pesquisa do autor. 


Efeitos principais dos parâmetros para a relação sinal ruído, segundo Phadke (1989): 


Gráfico de Efeitos Principais para Razões S/N 


Suco de laranja Abacaxi Maça Agua Açúcar 
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Sinal-ruído: Maior é melhor 


Figura 2 - Efeitos principais dos parâmetros em relação ao S/N 


Fonte: Dados colhidos mediante pesquisa do autor. 
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5 | CONCLUSÕES 


O método de Taguchi é uma ferramenta muito eficiente e prática para a otimização de 
produtos. O processo experimental fatorial com a matriz ortogonal permite analisar por meio 
de ensaios o comportamento e a resposta do produto a condições impostas pelo projetista. 

Neste trabalho podemos perceber que esta análise pode ser realizada por grandes 
empresas como multinacionais, por exemplo, pois além de aumentar a qualidade do 
produto se aproxima do cliente e de seus gostos, com isto, também se torna uma análise 
de mercado, tendo em vista que quanto maior a aceitação da sua bebida pelo palato dos 
clientes, mais vendas você terá e consequentemente mais lucro. 

O processo maximiza a produtividade pois busca eliminar ou atenuar os ruídos que 
causam a variabilidade e a falta de padrão dos produtos em sua fase final. A padronização 
de um serviço ou produto demonstra a perícia naquilo que está sendo trabalhado por parte 
da empresa, contribuindo para uma melhor experiência do cliente e alavancando o nome 


da organização no mercado de trabalho, evidenciando a qualidade nos serviços prestados. 
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RESUMO: O Ciclismo é a modalidade esportiva 
que acopla o atleta e uma bicicleta, envolvendo 
movimentos sincronizados entre uma série 
de articulações que proporcionam o esforço 
constante da região pélvica-lombar e membros 
inferiores, sendo uma corrida constante contra o 
tempo. O Ciclismo permite ao indivíduo a prática 
esportiva de forma individual ou em grupos. 
No Município de Porto Velho — RO é comum 
observar diversos grupos de Ciclistas ou grupos 
de “Pedal”, realizando percursos ou trilhas 
diversas através de toda extensão territorial da 
capital do estado de Rondônia. Em virtude dos 
fatos constatados pensou-se em desenvolver 
uma plataforma que pudesse proporcionar ao 
esportista uma experiência nova e repleta de 
interatividade, criando uma área onde pudesse 
desenvolver relacionamentos, buscar informação 
e que lhe gerasse praticidade, para uma atividade 
física mais segura. O aplicativo MUNDOBIKE 
tem por finalidade interligar todos os ciclistas de 
um bairro, município, estado ou até mesmo de 
um País. Tomando por base que a sociedade 
atual está cada vez mais integrada a plataformas 
digitais para a sua vivência do dia a dia, seja por 
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meio de redes sociais ou aplicativos que facilitem 
as suas tomadas de decisão, o MUNDOBIKE 
tem a meta de gerar conexão entre as pessoas, 
em uma sociedade que está se mostrando mais 
digitalizada a cada dia que passa. 
PALAVRAS-CHAVE: Aplicativo; 
Interligação; Virtual; MundoBike. 


Ciclismo; 


LINKING CYCLING 


ABSTRACT: Cycling is a sport that combines 
the athlete and a bicycle, involving synchronized 
movements between a series of joints that 
provide the constant effort of the pelvic-lumbar 
region and lower limbs, being a constant race 
against time. Cycling allows individuals to 
practice sporis individually or in groups. In the 
Municipality of Porto Velho — RO, it is common 
to observe several groups of Cyclists or groups 
of “Pedal”, performing different routes or trails 
through the entire territorial extension of the 
capital of the state of Rondonia. Due to the 
facts presented, it was thought to develop a 
platform that could provide the athlete with a new 
experience full of interactivity, creating an area 
where they could develop relationships, seek 
information and generate practicality, for a safer 
physical activity. The MUNDOBIKE app aims to 
connect all cyclists in a neighborhood, city, state 
or even a country. Based on the fact that today's 
society is increasingly integrated with digital 
platforms for their daily lives, either through social 
networks or apps that facilitate decision-making, 
MUNDOBIKE aims to generate connection 
between people, in a society that is showing itself 
to be more digitalized every day. 
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KEYWORDS: Apps; Cycling; Interconnection; Virtual; MundoBike. 


11 INTRODUÇÃO 


Com o advento da Internet, e a constante globalização, o surgimento de aplicativos 
que possam trazer conforto e praticidade pra vida do ser humano é cada vez mais comum. 
Em uma sociedade que é mais capitalista, tecnológica e materialista a cada dia que passa, 
constata-se que o ser humano busca ferramentas que lhe proporcione mais agilidade na 
hora de executar uma determinada tarefa, facilitando-lhe, o mínimo que seja, a sua vida. 

Segundo a própria definição do dicionário, aplicativo é “um programa de computador 
concebido para processar dados eletronicamente, facilitando e reduzindo o tempo de 
execução de uma tarefa pelo usuário”. Sendo assim, levando em consideração a crescente 
demanda social por um estilo de vida mais saudável, tendo cada vez mais pessoas optando 
por uma maneira mais “fit” de levar a vida, viu-se a possibilidade de se desenvolver um 
aplicativo focado em um nicho especifico e que gerasse essa praticidade e interatividade 
aos seus usuários. 

A partir de estudo feito entre ciclistas locados no município de Porto Velho — RO 
constatou-se que haveria espaço no mercado para um aplicativo que gerasse maior 
interligação entre ciclistas de uma determinada área proporcionando uma maior visibilidade 
a modalidade e aproveitando a facilidade de comunicação e integração das plataformas 
digitais pensou-se em criar um aplicativo que promovesse e desenvolvesse uma nova 
experiência dentro do esporte, ao mesmo tempo em que incentivasse a saúde pessoal e o 
bem estar social como um todo. 


21 O APLICATIVO 


Primeiramente o aplicativo irá nivelar a categoria do Ciclista, conforme a sua 
experiência com a modalidade, para que o usuário possa estar se identificando e progredindo 
conforme o uso da ferramenta. O aplicativo será denominado de “*MUNDOBIKE” e oferecerá 
ao seu usuário uma série de ferramentas que lhe permita desenvolver o seu conhecimento 
do esporte, dando-lhe dicas de como progredir e até mesmo lhe oferecendo uma rede de 
outros usuários, permitindo assim a interligação entre todos os praticantes do esporte no 
seu Estado. 

O MUNDOBIKE usará uma interface estilo mapa GPS, assim como outros aplicativos 
de localização (Waze, Google Maps, dentre outros) onde estarão ilustrados percursos e 
distâncias percorridos pelo ciclista, irá sinalizador grupos de ciclistas próximos ao usuário, 
dará sugestão de metas e desafios a serem executados pelo usuário e proporcionar a 
experiência mais completa possível. A princípio, o aplicativo irá dividir os usuários em 04 
categorias: Iniciante, Intermediário, Experiente e Profissional oferecendo diversas metas 
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ou desafios para que esteja incentivando e desafiando a evolução constante do ciclista 
sendo que conforme ele for completando essas tarefas ocorrerá a progressão de nível 
dentro do próprio aplicativo, assim mantendo a sua curiosidade pelo esporte e a promoção 


de uma vida mais saudável e menos sedentária. 


2.1 Iniciante 


Como o próprio nome já diz, aqui será enquadrado aquele indivíduo que está tendo 
um primeiro contato com o esporte ou que tenha pouca vivência dentro do Ciclismo. O 
aplicativo irá dar todas as dicas para tornar a experiência a mais descontraída possível, sem 
ter aquela pressão por resultado de categorias mais avançadas, de preferência apontando 
percursos com menor comprimento e de curta duração. O aplicativo agiria como um guia 


para o iniciante dando dicas de segurança e zelando pelo bem estar físico do usuário. 


2.2 Intermediário 


Neste estágio de evolução, o ciclista já possui um nível de conhecimento adequado 
do aplicativo e assim estará sujeito a percursos mais longos e com trechos que exigem uma 
afinidade maior do ciclista com a sua bicicleta. 


2.3 Experiente 


Em se tratando de ciclistas experientes, o aplicativo estará direcionando um usuário 
com mais vivência, e assim apontará grupos que percorrem trajetos mais longos ou trajetos 


que não sejam muito explorados dentro do aplicativo. 


2.4 Profissional 


O ciclista profissional está totalmente habituado a rotina de pedalar e praticamente 
transformou-a em um estilo de vida. O aplicativo exigirá metas de alta dificuldade do usuário 
e estabelecerá grupos mais restritos que tenham regularidade maior dentro do aplicativo 
com distâncias percorridas de ampla exigência e duração. O foco será proporcionar com 


maior dificuldade técnica e exigente ao usuário. 


31 RESULTADO E DISCUSSÕES 


Em se tratando da aplicabilidade da ferramenta, desenvolveu-se um estudo 
conforme necessidades de ciclistas quanto à viabilidade e interesse em uma possível 
plataforma digital que lhes permitisse encontrar grupos que dividissem o mesmo interesse 
que eles, sem que houvesse qualquer tipo de prejuízo a sua privacidade digital. Apesar 
da existência de ferramentas semelhantes ao aplicativo MUNDOBIKE, a forma com que o 
usuário utilizará o aplicativo e, sendo provocado por esse, para que atinja novos níveis de 
desempenho, é único e intrigará a curiosidade do público em um real interesse pelo uso 


da ferramenta. 
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41 CONCLUSÃO 


Assim, conforme apresentação, projeção e futura demonstração, se comprovará 
a real viabilidade comercial e econômica no desenvolvimento do aplicativo MUNDOBIKE 
para proporcionar uma experiência virtual mais completa ao esportista, com a finalidade 
de desenvolver o Ciclismo como prática esportiva benéfica a saúde e assim gerando uma 
comunidade mais ativa primeiramente no Município de Porto Velho — RO, futuramente 


visando uma maior visibilidade estadual e atingindo patamar nacional. 
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RESUMO: O presente trabalho analisou o 
resultado da avaliação dos cursos de Engenharia 
da Universidade Presbiteriana Mackenzie (UPM), 
sob a Ótica de seus egressos, no período de 
2012 a 2017, visando utilizar os resultados como 
ferramenta de gestão do ensino. A avaliação de 
um curso de graduação em Engenharia, sob a 
Ótica dos seus egressos, é um dos indicadores 
da qualidade do curso, pois estes têm uma visão 
ampla do processo de ensino e capacidade de 
verificar a contribuição do curso na sua formação 
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UPM 


pessoal e profissional, permitindo verificar se as 
metas de formação propostas estão sendo bem 
recebidas pelo mercado de trabalho globalizado. 
Para tanto, devem ser resgatados os conceitos 
relacionados a avaliação educacional e 
avaliação da educação superior no Brasil, prática 
desenvolvida por meio do Sistema Nacional de 
Avaliação da Educação Superior - SINAES, 
criado pela Lei nº 10.861, de 14 de abril de 
2004. Os resultados obtidos, fornecem subsídios 
para a crítica, a reflexão e o diálogo sobre o 
projeto pedagógico dos cursos oferecidos, 
apresentando tanto os seus diferenciais, que 
devem ser preservados e reforçados, quanto as 
suas fragilidades que indicam os aspectos que 
subsidiarão discussões que possam conduzir a 
melhorias dos cursos, devendo ser considerados 
como ferramentas de controle estratégico e 
operacional. 

PALAVRAS-CHAVE: Avaliação de curso. 
Percepção do Egresso. Qualidade dos Cursos 
de Engenharia. Formação. Mercado de Trabalho. 


PERCEPTION OF GRADUATE ABOUT 
THE QUALITY OF UPM ENGINEERING 
COURSES 


ABSTRACT: The present work analyzed the 
evaluation of the Engineering courses of the 
Mackenzie Presbyterian University, from the 
point of view of its graduates, in the period from 
2012 to 2017, aiming to use the results as a 
teaching management tool. The evaluation of 
an undergraduate course in Engineering, from 
the perspective of its graduates, is one of the 
indicators of the quality of the course. These 
have a vision of the teaching process and ability 
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to verify the contribution in their personal and professional formation, making it possible to 
verify whether the proposed training goals are being well received by the globalized labor 
market. Therefore, should be rescued the concepts related to the educational evaluation and 
evaluation of higher education in Brazil, a practice developed through the National System for 
the Evaluation of Higher Education - SINAES, created by Law No. 10,861 of April 14, 2004. 
The results obtained provide support for criticism, reflection and dialogue on the pedagogical 
project of the offered courses, presenting both their differentials, which must be preserved 
and reinforced, as well as their fragilities that indicate the aspects that will support discussions 
that may lead course improvementis, and should be considered as a strategic and operational 
control tool. 

KEYWORDS: Course evaluation. Perception of egress Quality of engineering courses. 
Formation. Job market. 


11 INTRODUÇÃO 
De acordo com SINAES (2004), 


Os principais objetivos da avaliação envolvem melhorar o mérito e o valor das 
instituições, áreas, cursos e programas, nas dimensões de ensino, pesquisa, 
extensão, gestão e formação; melhorar a qualidade da educação superior 
e orientar a expansão da oferta, além de promover a responsabilidade 
social das IES, respeitando a identidade institucional e a autonomia de cada 
organização. 

A avaliação é uma prática que possibilita a transparência das realizações 
institucionais, tanto para a comunidade acadêmica como para a sociedade, fornecendo 
subsídios para a análise de seus processos, possibilitando uma visão da problemática 
educacional implantada. Deve ser um processo sistemático e contínuo que tem por 
finalidade oferecer meios para o aprimoramento do currículo, apresentando, entre outros 
aspectos a associação entre avaliação e o processo de tomada de decisão (VIANNA, 2000). 

Como os processos avaliativos são instrumentos de gestão universitária capazes 
de indicar os melhores caminhos e formas de melhorias para um ensino de qualidade, a 
avaliação pelos egressos é uma das formas de se verificar a opinião deste grupo que já 
está exercendo suas atividades profissionais, possibilitando uma visão das transformações 
que ocorre com o egresso devido à influência exercida pelo currículo cursado, no seu 
desempenho como critério de avaliação de curso (VASCONCELOS, 2012). 

A Avaliação Institucional na UPM já está consolidada e vem sendo realizada desde 
1990, e de forma sistematizada desde 2001, objetivando a promoção da melhoria do ensino- 
aprendizagem e da avaliação global da eficiência da Universidade que se efetiva pela 
análise dos resultados obtidos (CASTANHEIRA; MASSON; MIRANDA, 2006; MIRANDA 
et al, 2011). 

O objetivo maior dos Cursos de Engenharia da UPM é proporcionar ao aluno uma 


excelente formação integral e, seus egressos, podem fornecer subsídios para avaliar o 
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alcance desse macro objetivo, sendo que o Projeto Didático Pedagógico (PPC), segundo 
Brito Cruz (2002), deve criar condições para capacitar seus egressos a criar e aplicar 
conhecimentos, promover a pesquisa cooperativa com empresas através do apoio em 
pesquisa e desenvolvimento (P&D), criar condições para a sustentabilidade em sua área 
específica, ser empreendedor e gestor. Muitos desses engenheiros deverão atuar no 
desenvolvimento de alta tecnologia em uma relação intima entre a tecnologia e a ciência. 

Assim, torna-se importante avaliar a percepção dos egressos para verificar se os 
objetivos estão sendo atingidos (VASCONCELOS; PEREIRA, 2015). 


21 AVALIAÇÃO REALIZADA PELOS EGRESSOS 


A avaliação é um tema que requer um detalhamento especial, pela sua própria 
natureza e principalmente pela função especial que cumpre dentro do currículo, que se 
fundamenta no princípio de que a aprendizagem não é alcançada espontaneamente pelo 
tecnicismo, mas, requer um processo acumulativo de assimilação reflexiva, de modo que 
o educando, a partir da reflexão sobre suas experiências e percepções iniciais, observe, 
reelabore e sistematize seu conhecimento acerca do objeto em estudo (DAVINI, 2008). 

Objetivando conhecer e acompanhar o percurso dos egressos, quanto ao mercado 
de trabalho, à satisfação atual com o curso realizado e à imagem da UPM, a Comissão 
Própria de Avaliação-CPA desenvolveu um instrumento de investigação para ser aplicado 
periodicamente. 

A importância da pesquisa com os egressos tem ganhado relevância na medida 
em que se percebe que o desempenho desse grupo pode ser considerado como um dos 
indicadores da qualidade dos cursos oferecidos pela Escola de Engenharia (Civil, Elétrica, 
de Materiais, Mecânica e de Produção), bem como do seu compromisso social, pois um 
dos objetivos da Universidade é a preparação de seu aluno para o mundo do trabalho. 

Como parte do instrumento, o egresso é convidado a manifestar-se sobre a 
contribuição da Universidade na sua formação profissional e no seu desenvolvimento 
cultural. A pesquisa anterior teve uma abrangência de um período de 10 anos (egressos 
desde 2005/2º semestre a 2015/1º semestre), e a pesquisa atual, realizada em 2017, 
teve uma abrangência de 5 anos, correspondendo aos egressos formados no período de 
2012/1ºsemestre a 2016/2º semestre. 

O método de pesquisa utilizado foi o mail survey, que se caracteriza pelo 
questionamento direto das pessoas cuja opinião ou comportamento se deseja conhecer. 
Neste caso os levantamentos são utilizados quando o pesquisador deseja responder 
acerca da distribuição de uma variável ou de maneira como ocorrem em as relações entre 
características de pessoas ou grupos, de forma como ocorrem em situações naturais, 
sendo esta a estratégia mais apropriada para a análise de fatos e descrições (MARTINS; 
THEOPHILO, 2009). 
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Além disso, este tipo de pesquisa facilita a obtenção de um maior número de 
respondentes (UMBACH, 2004). Neste caso, os respondentes são os egressos de cursos 
de graduação da Escola de Engenharia da UPM. Os resultados estão divulgados em escala 
de O a 1 (onde 1 é o patamar máximo de qualidade — aprovação máxima do indicador; e O 
é o patamar mínimo de qualidade — reprovação absoluta do indicador). 

Os indicadores apresentados de O a 1 representam as porcentagens de repostas 
positivas, retratando uma indicação de qualidade e o nível de satisfação do egresso sobre 
o tipo de serviço prestado pela universidade e a percepção da sociedade sobre a Escola 
de Engenharia da UPM. 

O questionário utilizado é composto de 14 questões fechadas e 1 questão aberta 


(sugestões), de acordo com a Tabela 1. 


Tipos Número de Questões 
Sobre Profissão e Mercado de Trabalho 7 
Sobre Profissão e o Mackenzie 3 
Satisfação com o Mackenzie 4 


Tabela 1 : Tipos de Questão 


Considerando o instrumento, o objetivo desse momento avaliativo foi conhecer 
a opinião do egresso sobre a qualidade da preparação oferecida pela UPM para o 
enfrentamento dos desafios da sociedade em geral e do complexo mundo coorporativo. 


31 RESULTADOS OBTIDOS 


O processo avaliativo foi realizado em dois meses, abrangendo o envio do 
instrumento aos egressos e o fechamento da coleta dos dados. Dos 3273 alunos convidados 
a responder o instrumento, 658 responderam, ou seja, obteve-se uma adesão de 20,13%. 

Para um intervalo de confiança de 95% da pesquisa, ou seja, se todos os indivíduos 
tivessem a mesma probabilidade de responder ao instrumento, a margem de erro estimada, 
para um índice de confiança de 95%, seria de 3,4%, mas neste caso deve-se levar em 
conta os erros estatísticos devido a dois fatores: não foram enviados e-mail para todos os 
formandos (somente para aqueles com os e-mails ativos), e nem todos responderam. 

Comparando o número de participantes da pesquisa atual (20,13%) com a realizada 
anteriormente (20,67%) observou-se que o percentual de participantes foi praticamente 
igual (decréscimo de 0,54%). 

O instrumento é composto por dois indicadores sensoriais que auxiliam na 
construção de uma base de dados para as análises que serão realizadas. A esses dois 
indicadores seguem três referentes ao ingresso no mundo do trabalho, dessa forma os 
cinco primeiros indicadores são de natureza objetiva e não abordam a percepção do 
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egresso. A partir do sexto indicador busca-se a opinião do egresso sobre a qualidade 
percebida sobre a formação recebida assim como o modo como a sociedade, em especial 
o mercado de trabalho, percebe a UPM. Os últimos indicadores referem-se a situação atual 
do egresso no mundo corporativo. Esses indicadores podem ser associados a um critério 
de empregabilidade do egresso. 

As Tabelas 2 a 4 apresentam os resultados das respostas válidas sobre as 
informações gerais. 


Sexo Percentual (%) 





Tabela 2: Sexo do Egresso 


Os resultados mostram que ainda o percentual feminino é menor que o masculino 
na área das engenharias. 


Idade Percentual (%) 
<25 anos 241 36,63 
26 a 29 anos 282 42,86 
mais de 50 anos E 0,00 


Tabela 3: Idade atual do Egresso 


Comparando os resultados obtidos na pesquisa atual com a realizada anteriormente, 
observou-se que houve um aumento no percentual de egressos com faixa etária menor que 25 
anos (aumento de 8,98%), bem como na faixa etária entre 26 a 29 anos (aumento de 3,80%). 


Estado Percentual (%) 
SP 624 94,83 
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Exterior 


Sem Resposta 





Tabela 4: Estado onde o Egresso reside atualmente. 


Comparando os resultados obtidos na pesquisa atual com a realizada anteriormente, 
observou-se que houve um aumento no percentual de egressos com faixa etária menor que 
25 anos (aumento de 8,98%), bem como na faixa etária entre 26 a 29 anos (aumento de 
3,80%). 

Predominantemente os alunos residem na região sudeste (em São Paulo a maioria) 
e apenas 1,67% (11 alunos) estão residindo no exterior. Comparando os resultados obtidos 
na pesquisa atual com a anterior, houve um aumento no percentual de egressos residindo 
no exterior (aumento de 83%), e um decréscimo de no percentual dos egressos que residem 
em São Paulo. 


A Tabela 5 apresenta o tempo entre a formatura e o início da atividade profissional. 


Tempo Quantidade de respostas Percentual (%) 
Já estava exercendo atividade 





profissional antes da formatura 435 66,11 
Menos de 1 ano 18,09 
De 1a2 anos 10,18 
De2a3 anos 2,89 
De 3a4 anos o 9 1,37 
Mais de 4 anos 0,76 
Sem Resposta 0,61 





Tabela 5: Tempo entre a formatura e o início da atividade profissional. 


Embora a variação do tempo entre a formatura e o início da atividade profissional 
tenha sido praticamente mantida nesta pesquisa, observou-se que houve um discreto 
aumento no tempo para o egresso ser absorvido pelo mercado de trabalho, quando se 
compara os resultados obtidos na pesquisa anterior. 


A Tabela 6 apresenta o tipo de empresa em que o egresso trabalha. 
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Tipo de Empresa Percentual (%) 





Empresa privada 58,06 
Profissional liberal ou autônomo 10,48 
Proprietário de empresa ou firma individual 4,71 


Outros 26,175 





Tabela 6: Natureza do emprego atual 


Comparando os resultados obtidos na pesquisa atual com a realizada anteriormente, 
observou-se que houve um decréscimo no percentual dos egressos sendo absorvidos por 
empresas privadas (-2,21%) e que sejam proprietários de empresas ou firma individual 
(-1,44%), e um aumento em atividades como profissional liberal ou autônomo (+1,12%), 
o que aponta para um aumento da absorção do egresso em atividades profissionais 
empreendedoras. 


A Tabela 7 apresenta a renda média mensal dos egressos. 


Renda Percentual (%) 
(R$) respostas 
Até 5.000 227 34,50 
De 5.001 até 7.000 293 44,53 
De 7.001 até 9.000 58 8,81 


Acima de 9.500 53 8,05 
Sem Resposta 21 4,10 





Tabela 7: Renda Média Mensal (bruta) em reais. 


Comparando os resultados obtidos na pesquisa atual com a realizada anteriormente, 
observa-se que houve um aumento do percentual de egressos com salários menores que 
R$5000,00 (9,31%), embora na pesquisa anterior o universo da mesma tenha atingido 
egressos com maior tempo exercendo a atividade profissional. 


3.1 Indicadores institucionais: avaliações dos anos de 2015 e 2017 


A Tabela 8 apresenta os resultados obtidos para o grau de satisfação do Egresso em 
relação à Instituição nas pesquisas realizadas nos anos de 2015 e 2017 (questões 4 a 13). 
O indicador varia de O (nota mínima) a 1 (nota máxima). 
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Perguntas Avaliação 2015 | Avaliação 2017 | Variação 


(2ºS/2006- (1ºS/2012- 
1ºS/2015) 2ºS/2017) 


4. Qual o nível de satisfação com a sua ; 
profissão? (Alto) 0,58 os 0,02 





5. Você se sentiu preparado para o mundo do 
trabalho quando se formou? (Sim , muito) 0,78 0,74 +0,01 


6. O conjunto de componentes curriculares 
contribuiu para o seu desempenho profissional? -0,03 
(Muito) 


7.0 curso como um todo colaborou para o seu 0,73 0,77 +0,04 
desenvolvimento cultural e social? (Sim , muito) 


8. Qual o conceito que você atribui aos 0,89 0,88 001 
professores do curso que você fez”? 


9. Você mantem algum contato atual com a +0.05 
Universidade Presbiteriana Mackenzie? (Sim) , 
10.Você escolheria a UPM novamente para +0.04 
realizar um curso de pós-graduação? (Sim) : 
11. Como você avalia a imagem do Mackenzie ? e -0,02 
12. Qual o seu nível de satisfação com a sua -0.05 
situação profissional atual no aspecto financeiro? ú 
13. Qual o seu nível de satisfação com a sua -0.04 
situação atual, no aspecto social? 0,8 0,7 Ê 


Tabela 8: Grau de satisfação do Egresso em relação à Instituição. 


Embora a média dos indicadores seja maior que 70%, observa-se que há indicadores 
que exigem atenção. 

Os indicadores que exigem maior atenção são aqueles correspondentes as 
questões 4 e 9 que apresentaram índices de satisfação menor ou igual a 60%. 

Comparando-se os indicadores institucionais obtidos no ano de 2016 em relação 
ao ano de 2015, observa-se que dos 10 indicadores avaliados houve um decréscimo em 7 
indicadores. 

Salienta-se que entre os indicadores que sofreram acréscimo foram: questão 5 
(Você se sentiu preparado para o mundo do trabalho quando se formou”), questão 7 (O 
curso como um todo colaborou para o seu desenvolvimento cultural e social?), questão 
9 (Você mantém algum contato atual com a Universidade Presbiteriana Mackenzie?) e 
questão 10(Você escolheria o Mackenzie novamente para realizar um curso de pós- 
graduação?), com aumentos de 0,01; 0,04; 0,05 e 0,04%, respectivamente. 

A queda média nos indicadores foi de -0,01. 

A maior queda foi atribuída à questão 12 (Qual o seu nível de satisfação com a sua 
situação profissional atual no aspecto financeiro?), cujo decréscimo foi de -0,05. 

Em relação ao nível de satisfação profissional e à percepção quanto à formação 
obtida na universidade para o mundo do trabalho, os indicadores obtidos foram abaixo de 
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0,60, significando que grande parte dos egressos avaliou esse indicador como negativo, 
sendo que estes apresentaram uma queda em relação à avaliação anterior. 

Os resultados continuam apontando que se deve dar uma atenção especial para 
uma efetiva preparação para o mundo do trabalho, sendo que este aspecto precisa ser 
fortemente considerado na revisão dos PPCs, bem como na modernização das novas 
matrizes curriculares, levando-se em conta as novas tendências e exigências profissionais 
atuais. 

Os resultados obtidos mostram que 71% dos alunos avaliam de forma satisfatória 
a formação recebida, sendo que tanto o atendimento como as estruturas curriculares são 
consideradas fatores que atenderam às necessidades do aluno. 

A satisfação do egresso com a Instituição, bem como a sua qualidade são atestadas 
pelos índices obtidos para a intenção de complementação da formação na UPM. 


3.2 Índice de satisfação geral 


Neste item são abordados dois parâmetros em que se agregam as questões em 
médias aritméticas para dois conjuntos de séries: satisfação geral com trabalho (questões: 
4,12 e 13) e satisfação geral com o Curso/UPM (questões: 5,6,7,8,10, 11). A Tabela 9 e o 
Gráfico da Figura 1 apresentam os resultados obtidos para os índices de satisfação geral 


por curso. 


Satisfação geral com Satisfação geral com o 
Curso trabalho (questões: 4, 12 Curso/UPM (questões: 
e 13) 5,6,7,8,10, 11) 


Eng. Materiais 069 | 0,88 
Eng. Mecânica 0,83 
Eng. Produção 0,75 


Tabela 9: Índices de satisfação geral. 
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Eng. Eng. Eng. 
Materiais Mecânica Produção 


Eng. Civil Eng. Elétrica 


E Satisfação geral com trabalho (questões: 4, 12 e 13) 


= Satisfação geral com Mackenzie (questões: 5,6,7,8.10, 11) 


Figura 1 — Gráfico dos Índices de satisfação geral. 


Por meio dos resultados pode-se observar ue na Escola de Engenharia da UPM, o 
curso que apresenta o maior índice de satisfação com o trabalho é o Curso de Engenharia 
de Materiais (0,69) e o que apresenta o menor índice de satisfação é o Curso de Engenharia 
Elétrica (0,55). 

Quanto ao índice de satisfação com o Mackenzie o curso que apresenta o maior 
índice de satisfação é o Curso de Engenharia de Materiais (0,88) e o que apresenta o 
menor índice de satisfação é o Curso de Engenharia Elétrica (0,67). 


3.3 Grau de empregabilidade 


A Tabela 10 e o Gráfico da Figura 2 apresentam os resultados obtidos para o grau 
de empregabilidade. 


Curso Sim, na área | Sim, mas fora da 

de formação | área de formação 
Engenharia Civil 52,22% 18,89% 21,18% 1,11% 
Engenharia Elétrica 43,90% 21,95% 31,/2% 2,44% 


Sem resposta 

Engenharia de Materiais 44,74% 39,47% 15,79% (0) 
Engenharia Mecânica 48,84% 34,88% 16,28% 0 

Engenharia de Produção 50,00% 41,91% 6,06% 2,03% 


Média EE 49,95% 21,85% 21,27% 0,93% 


Tabela 10: Grau de empregabilidade (exercício de alguma atividade profissional). 
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Eng. Civil Eng. Elétrica Ene. Eng. Eng. 
Materiais Mecânica Produção 


m Satisfação geral com trabalho (questões: 4, 12 e 13) 


m Satisfação geral com Mackenzie (questões: 5,6,7,8,10, 11) 


Figura 2 — Gráfico do Grau de empregabilidade (exercício de alguma atividade profissional). 


Os resultados obtidos mostram que cerca de 50% dos alunos trabalham na área 
de sua formação e 77,80% se encontram empregados. O curso que apresenta o maior 
grau de empregabilidade é o curso de Engenharia de Produção (91,91%) e o com menor 
empregabilidade é o de Engenharia Elétrica (65,85%). 

Os egressos dos cursos de engenharia são os que vem sofrendo maior impacto 
com a crise econômica pela qual o país vem passando, pois o mercado de trabalho para 
este profissional foi um dos mais afetados. Neste contexto, a ênfase que vem sendo 
implementada em toda a universidade para a educação empreendedora, torna-se um 
importante diferencial para a inserção do egresso da UPM no mercado de trabalho. 


41 CONSIDERAÇÕES FINAIS 


Os processos avaliativos são fundamentais para a melhoria dos cursos de Engenharia, 
pois por meio das informações obtidas nesses levantamentos, os projetos pedagógicos são 
atualizados, incorporando as solicitações globalizadas exigidas pelo mercado de trabalho, 
com enfoque nas observações dos alunos e dentro do contexto institucional. 

Pela importância socioeconômica da engenharia, o seu ensino deve ser coerente com 
as aspirações da sociedade onde se insere. Para isso, deverá atualizar-se constantemente 
criando condições para o pleno desenvolvimento técnico e humanístico de seus alunos, que 
devem estar preparados para enfrentar situações diversas, conscientes de suas expertises. 

Desta forma a avaliação é a Única forma de assegurar o acompanhamento do ensino 
e da aprendizagem, da gestão, da pesquisa e das atividades extracurriculares preservando 
assim a identidade institucional e sua capacidade de se adaptar às exigências do mercado 
sem perder de vista a formação integral do aluno. 
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As avaliações realizadas com os egressos da Escola de Engenharia da UPM, em 
2017, utilizando a metodologia mail survey contemplaram um número significativo de 
alunos (cerca de 20%). O ideal seria pelo menos 40%, o que é extremamente difícil de se 
conseguir depois de algum tempo, pois os dados são dinâmicos. 

Os resultados obtidos mostraram que a satisfação positiva do egresso com a 
Instituição, bem como a sua qualidade são atestadas pelos índices obtidos para a intenção 
de complementação da formação. 

Os projetos didático-pedagógicos dos cursos de engenharia da UPM têm sido 
atualizados periodicamente considerando os resultados destas avaliações e de outros 
estudos e análises que envolvem dados levantados como o perfil do aluno que se deseja 
formar e com as ações necessárias para cumprir os objetivos, deixando explicito o perfil do 
profissional a ser formado. 

Para atingir estes objetivos, dentre as ações implementadas, citam-se: 
reformulação e flexibilização curricular com revisão da estrutura acadêmica buscando o 
aumento da qualidade; mudança de horário de aulas, de forma a envolver mais o aluno 
como a implantação de atividades integradoras ao longo de todo o curso, incentivando 
o protagonismo estudantil, visando despertar maior interesse pela sua área escolhida; 
maior incentivo ao empreendedorismo, estímulo à inovação, bem como a participação em 
programas de iniciação científica, além da implantação de práticas pedagógicas modernas 
objetivando o decréscimo nas taxas de evasão. 

Destaca-se também o incentivo à internacionalização, por meio dos convênios com 
Universidades Estrangeiras, que permitem aos alunos vivenciar novas realidades, novos 
métodos de ensino, proporcionando uma formação mais abrangente, além de propiciar a 


dupla titulação que possibilita o futuro exercício profissional fora do país. 
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